Title: OVERCOMING CATASTROPHIC INTERFERENCE USING CONCEPTOR-AIDED BACKPROPAGATION

Abstract: Catastrophic interference has been a major roadblock in the research of continual learning. Here we propose a variant of the back-propagation algorithm, "conceptor-aided backprop" (CAB), in which gradients are shielded by conceptors against degradation of previously learned tasks. Conceptors have their origin in reservoir computing, where they have been previously shown to overcome catastrophic forgetting. CAB extends these results to deep feedforward networks. On the disjoint and permuted MNIST tasks, CAB outperforms two other methods for coping with catastrophic interference that have recently been proposed.

1. INTRODUCTION:

Agents with general artificial intelligence are supposed to learn and perform well on multiple tasks. Continual learning refers to the scenarios where a machine learning system can retain previously acquired skills while learning new ones. However, when trained on a sequence of tasks, neural networks usually forget about previous tasks after their weights are adjusted for a new task. This notorious problem known as catastrophic interference (CI) (McCloskey & Cohen, 1989; Ratcliff, 1990; French, 1999; Kumaran et al., 2016) poses a serious challenge towards continual learning.

Many approaches have been proposed to overcome or mitigate the problem of CI in the last three decades (Hinton & Plaut, 1987; French, 1991; Ans & Rousset, 1997; French, 1997; Srivastava et al., 2014) . Especially recently, an avalanche of new methods in the deep learning field has brought about dramatic improvements in continual learning in neural networks. Kirkpatrick et al. (2017) introduced a regularization-based method called elastic weight consolidation (EWC), which uses the posterior distribution of parameters for the old tasks as a prior for the new task. They approximated the posterior by a Gaussian distribution with the parameters for old tasks as the mean and the inverse diagonal of the Fisher information matrix as the variance. Lee et al. (2017) introduced two incremental moment matching (IMM) methods called mean-IMM and mode-IMM. Mean-IMM approximates the distribution of parameters for both old and new tasks by a Gaussian distribution, which is estimated by minimizing its KL-divergence from the mixture of two Gaussian posteriors, one for the old task and the other one for the new task. Mode-IMM estimates the mode of this mixture of two Gaussians and uses it as the optimal parameters for both tasks.

In the field of Reservoir Computing (Jaeger, 2001; Maass et al., 2002) , an effective solution to CI using conceptors was proposed by Jaeger (2014) to incrementally train a recurrent neural network to generate spatial-temporal signals. Conceptors are a general-purpose neuro-computational mechanism that can be used in a diversity of neural information processing tasks including temporal pattern classification, one-shot learning, human motion pattern generation, de-noising and signal separation (Jaeger, 2017) . In this paper, we adopt and extend the method introduced in Jaeger ( 2014) and propose a conceptor-aided backpropagation (CAB) algorithm to train feed-forward networks. For each layer of a network, CAB computes a conceptor to characterize the linear subspace spanned by the neural activations in that layer that have appeared in already learned tasks. When the network is trained on a new task, CAB uses the conceptor to adjust the gradients given by backpropagation so that the linear transformation restricted to the characterized subspace will be preserved after the gradient descent procedure. Experiment results of two benchmark tests showed highly competitive performance of CAB.

The rest of this paper is structured as follows. Section 2 introduces conceptors and their application to incremental learning by ridge regression. Section 3 extends the method to stochastic gradient descent and describes the CAB algorithm. Section 4 compares its performance on the permuted and disjoint MNIST tasks to recent methods that address the same problem. Finally we conclude our paper in Section 5.

2. INCREMENTAL RIDGE REGRESSION BY CONCEPTORS:

This section reviews the basics of conceptor theory and its application to incrementally training linear readouts of recurrent neural networks as used in reservoir computing. A comprehensive treatment can be found in (Jaeger, 2014) . In brief, a matrix conceptor C for some vector-valued random variable x ∈ R N is defined as a linear transformation that minimizes the following loss function.

2.1. CONCEPTORS:

E x [||x -Cx|| 2 ] + α -2 ||C|| 2 fro (1)

where α is a control parameter called aperture and || • || fro is the Frobenius norm. This optimization problem has a closed-form solution

C = R(R + α -2 I) -1

(2)

where

R = E x [xx ]

is the N × N correlation matrix of x, and I is the N × N identity matrix. This result given in (2) can be understood by studying the singular value decomposition (SVD) of C. If R = U ΣU is the SVD of R, then the SVD of C is given as U SU , where the singular values s i of C can be written in terms of the singular values σ i of R:

s i = σ i /(σ i + α -2 ) ∈ [0, 1).

In intuitive terms, C is a soft projection matrix on the linear subspace where the samples of x lie. For a vector y in this subspace, C acts like the identity: Cy ≈ y, and when some noise orthogonal to the subspace is added to y, C de-noises: C(y + ) ≈ y. Figure 1 shows the ellipsoids corresponding to three sets of R 3 points. We define the quota Q(C) of a conceptor to be the mean singular values:

Q(C) := 1 N N i=1 s i .

Intuitively, the quota measures the fraction of the total dimensions of the entire vector space that is claimed by C.

Moreover, logic operations that satisfy most laws of Boolean logic can be defined on matrix conceptors as the following:

EQUATION

C i ∨ C j :=(R i + R j )(R i + R j + α -2 I) -1 (4) C i ∧ C j :=¬(¬C i ∨ ¬C j ) (5)

where ¬C softly projects onto a linear subspace that can be roughly understood as the orthogonal complement of the subspace characterized by C. C i ∨ C j is the conceptor computed from the union of the two sets of sample points from which C i and C j are computed. It describes a space that is approximately the sum of linear subspaces characterized by C i and C j , respectively. The definition of C i ∧ C j reflects de Morgan's law. Figure 2 illustrates the geometry of these operations.

Figure 2 : Geometry of Boolean operations on 2-dimensional conceptors. The OR (resp. AND) operation gives a conceptor whose ellipsoid approximately is the smallest (largest) ellipsoid enclosing (contained in) the argument conceptor's ellipsoids.

2.2. INCREMENTAL RIDGE REGRESSION:

This subsection explains how conceptors can be applied to master continual learning in a simple linear model trained on a supervised task by ridge regression. The training is done sequentially on multiple input-to-output mapping tasks. This simplified scenario illustrates the working principle of continual learning with conceptors and will later be used repeatedly as a sub-procedure in the CAB algorithm for training multilayer feed-forward networks.

Consider a sequence of m incoming tasks indexed by j. We denote the training dataset for the j-th task by {(x j 1 , y j 1 ), • • • , (x j n , y j n )}, where x j i ∈ R N are input vectors and y j i ∈ R M their corresponding target outputs. Whenever the training dataset for a new task is available, the incremental learning method will compute a matrix conceptor C j for the input variable of the new task using Equation 2and update the linear model, resulting in a sequence of linear models W 1 , . . . W m such that W j solves not only the j-th task but also all previous tasks: for k ≤ j, y k ≈ W j x k . The conceptor C j is a soft projection matrix onto the linear subspace spanned by input patterns from the j-th task. Then,

A j-1 = C 1 ∨ • • • ∨ C j-1

characterizes the memory space already claimed by the tasks 1, . . . , j -1 and F j = ¬A j-1 , the orthogonal complement of A j -1, represents the memory space still free for the j-th task. Here "memory space" refers to the linear space of input vectors. In detail, this method proceeds in the following way:

• Initialization (no task trained yet):

W 0 = 0 M ×N , A 0 = 0 N ×N .

• Incremental task learning: For tasks j = 1, . . . , m do:

1. Store the input vectors from the j-th training dataset of size n into a N × n sized input collection matrix X j , and store the output vectors into a M ×n sized output collection matrix Y j . 2. Compute the conceptor for this task by C j = R j (R j + α -2 I) -1 , where R j = 1 n X j X j 3. Train an increment matrix W j inc (to be added to W j-1 , yielding W j ), with the crucial aid of a helper conceptor F j : (a) F j := ¬A j-1 (comment: this conceptor characterizes the "still disposable" memory space for the j-th task), (b) T := Y j -(W j-1 X j ) (comment: this matrix consists of target values for a linear regression to compute W j inc ), (c) S := F j X j (comment: this matrix consists of input arguments for the linear regression), (d) W j inc = ((SS /n + λ -2 I) -1 ST /n) (comment: carry out the regression, regularized by λ -2 ), 4. Update W j : W j = W j-1 + W j inc . 5. Update A : A j = A j-1 ∨ C j (comment: this is possible due to the associativity of the ∨ operation on conceptors)

The weight increment W j inc does not interfere much with the previously learned weights W j-1 because the regularization in step 3(d) constrains the row space of W j inc to be only the linear subspace spanned by input arguments defined in 3(c), which are inside the kernel of W j-1 due to the projection by F j . Intuitively speaking, when learning a new task, this algorithm exploits only the components of input vectors in the still unused space (kernel of W j-1 , characterized by F j ) to compensate errors for the new task and leaves the directions in the already used memory space (row space of W j-1 , characterized by A j-1 ) intact.

3. CONCEPTOR-AIDED SGD AND BACK-PROP:

In this section, we first derive a stochastic gradient descent version of the algorithm described in the previous section, then present the procedure of CAB.

3.1. SGD:

In the algorithm introduced in the previous section, W j inc is computed by ridge regression, which offers a closed-form solution to minimize the following cost function

EQUATION

)

where t = y j -W j-1 x j , s = F j x j . One can also minimize this cost function by stochastic gradient descent (SGD), which starts from an initial guess of W j inc and repeatedly performs the following update

EQUATION

)

where η is the learning rate and the gradient is given by:

EQUATION

)

Substituting t by y j -W j-1 x j and s by F j x j = (I -A j-1 )x j in (8), we get

EQUATION

Due to the regularization term in the cost function, as the optimization goes on, eventually W inc will null the input components that are not inside the linear subspace characterized by F j , hence W j inc A j-1 x j will converge to 0 as the algorithm proceeds. In addition, since W j = W j-1 + W j inc , (10) can be simplified to

EQUATION

)

Adding W j-1 to both sides of (7), we obtain the update rule for W j :

EQUATION

)

where e := W j x j -y j . In practice, at every iteration, the expected value can be approximated by a mini-batch of size n B , indexed by i B :

Ê[es ] = 1 n B L i B =0 (W j x j i B -y j i B )(F j x j i B ) = 1 n B L i B =0 (W j x j i B -y j i B )x j i B F j (13)

where the transpose for F j can be dropped since it is symmetric.

If we only train the j-th task without considering the previous tasks, the update rule given by normal SGD is

W j ← W j -2ηE[ex j ] + 2ηλ -2 W j (14)

Comparing this to the update rule in (12), we notice two modifications when a conceptor is adopted to avoid CI: first, the gradient of weights are calculated using the conceptor-projected input vector s = F j x j instead of the original input vector x j ; second, regularization is done on the weight increment W j inc rather than the final weight W j . These two modifications lead to our design of the conceptor-aided algorithm for training multilayer feed-forward networks.

3.2. BACKPROP:

The basic idea of CAB is to guide the gradients of the loss function on every linear component of the network by a matrix conceptor computed from previous tasks during error back-propagation (Rumelhart et al., 1986) , repeatedly applying the conceptor-aided SGD technique introduced in the previous section in every layer.

Consider a feed-forward network with L + 1 layers, indexed by l = 0, . . . L, such that the 0-th and the L-th layers are the input and output layers respectively. W (l) represents the linear connections between the (l -1)-th and the l-th layer, where we refer to the former as the pre-synaptic layer with respect to W (l) , and to the latter as the post-synaptic layer. We denote by N (l) the size of the l-th layer (excluding the bias unit) and A (l) j a conceptor characterizing the memory space in the l-th layer used up by the first j tasks. Let σ(•) be the activation function of the nonlinear neurons and θ all the parameters of the network to be trained. Then the incremental training method with CAB proceeds as follows:

• Initialization (no task trained yet): ∀l = 0, . . . , L -1, A (l) 0 := 0 (N (l) +1)×(N (l) +1) , and randomly initialize W (l+1) 0 to be a matrix of size N (l+1) × (N (l) + 1).

• Incremental task learning: For j = 1, . . . , m do:

1. ∀l = 0, . . . , L -1, F (l) j = ¬A (l) (j-1) . (This conceptor characterizes the still disposable vector space in layer l for learning task j)

2. Update the network parameters θ (j-1) obtained after training the first j -1 tasks to θ j by stochastic gradient descent, where the gradients are computed by CAB instead of the classical backprop. Algorithms 1 and 2 detail the forward and backward pass of CAB, respectively. Different from classical backprop, the gradients are guided by a matrix conceptor F (l) j , such that in each layer only the activity in the still disposable memory space will contribute to the gradient. Note that the conceptors remain the same until convergence of the network for task j.

3. After training on the j-th task, run the forward procedure again on a batch of n B input vectors, indexed by i B , taken from the j-th training dataset, to collect activations h

(l) i B j

of each layer into a N (l) × n B sized matrix H (l) j , and set the correlation matrix

R (l) j = 1 n B H (l) j (H (l) j ) .

4. Compute a conceptor on the l-th layer for the j-th pattern by C (l) j = R (l) j (R (l) j + α -2 I N (l) ×N (l) ) -1 , ∀l = 0, . . . , L -1. Finding an optimal aperture can be done by a cross-validation search 1 .

5. Update the conceptor for already used space in every layer: A (l) j = A (l) j ∨ C (l) j , ∀l = 0, . . . , L -1.

Algorithm 1 The forward procedure of conceptor-aided backprop, adapted from the traditional backprop. Input vectors are passed through a feed-forward network to compute the cost function. L(ŷ j , y j ) denotes the loss for the j-th task, to which a regularizer Ω(θ j inc ) = Ω(θ j -θ j-1 ) = ||θ j -θ j-1 || 2 fro is added to obtain the total cost J , where θ contains all the weights (biases are considered as weights connected to the bias units). The increment of parameters rather than the parameters themselves are regularized, similar to the conceptor-aided SGD. Require: Network depth, l Require: W (l) j , l ∈ {1, . . . , L}, the weight matrices of the network Require: x j , one input vector of the j-th task Require: y j , the target output for x j 1: h (0) = x j 2: for l = 1, . . . L do -1) , 1] , include the bias unit 4: l) 5:

3: b (l) = [h (l

a (l) = W (l) j b (

h (l) = σ(a (l) ) 6: end for 7: ŷj = h (l) 8: J = L(ŷ j , y j ) + λΩ(θ j inc )

Algorithm 2 The backward procedure of conceptor-aided backprop for the j-th task, adapted from the traditional backprop. The gradient g of the loss function L on the activations a (l) represents the error for the linear transformation W (l) j between the (l -1)-th and the l-th layers. In the standard backprop algorithm, the gradient of L on W (l) j is computed as an outer product of the post-synaptic errors g and the pre-synaptic activities h (l-1) . This resembles the computation of the gradient in the linear SGD algorithm, which motivates us to apply conceptors in a similar fashion as in the conceptor-aided SGD. Specifically, we project the gradient ∇ W (l) j L by the matrix conceptor F (l-1) j that indicates the free memory space on the pre-synaptic layer.

1:

g ← ∇ ŷ J = ∇ ŷ L(ŷ, y) 2: for l = L, L -1, . . . , 1 do 3:

Convert the gradient on the layer's output into a gradient on the pre-nonlinearity activation ( denotes element-wise multiplication):

g ← ∇ a (l) J = g σ (a (l) ) 4:

Compute the gradient of weights, project it by F (l-1) j , and add it to the regularization term on the increment:

∇ W (l) j J =g(F (l-1) j b (l-1) ) + λ∇ W (l) j Ω(θ j inc ) = gb (l-1) F (l-1) j + 2λW (l) inc j =gb (l-1) F (l-1) j + 2λ(W (l) j -W (l) j-1 ) 5:

Propagate the gradients w.r.t. the next lower-level hidden layers activations:

g ← ∇ h (l-1) J = W (l) j g 6: end for To test the performance of CAB, we evaluated it on the permuted MNIST experiment (Srivastava et al., 2013; Goodfellow et al., 2014; Kirkpatrick et al., 2017; Lee et al., 2017) , where a sequence of pattern recognition tasks are created from the MNIST dataset (LeCun et al., 1998) . For each task, a random permutation of input image pixels is generated and applied to all images in MNIST to obtain a new shuffled dataset, equally difficult to recognize as the original one, the objective of each task is to recognize these images with shuffled pixels.

For a proof-of-concept demonstration, we trained a simple but sufficient feed-forward network with [784-100-10] of neurons to classify 10 permuted MNIST datasets. The network has logistic sigmoid neurons in both hidden and output layers, and is trained with mean squared error as the cost function.

Vanilla SGD was used in all experiments to optimize the cost function. Learning rate and aperture were set to 0.1 and 4, respectively. For comparison, we also tested EWC on the same task with the same network architecture, based on the implementation by Seff (2017) . The parameters chosen for the EWC algorithm were 0.01 for the learning rate and 15 for the weight of the Fisher penalty term. Figure 3 shows the performance of CAB on this task, the average testing accuracy is 95.2% after learning all 10 tasks sequentially. Although a fair amount of effort was spent on searching for optimal parameters for EWC, the accuracies shown here might still not reflect its best performance. However, the same experiment with EWC was also conducted in Kemker et al. (2017) , where the authors reimplemented EWC on a network with higher capacity (2 hidden layers and 400 ReLU neurons per layer) and the resulting average accuracy after learning 10 tasks sequentially was shown to be around 93%.

Since all tasks are generated by permuting the same dataset, the portion of the input space occupied by each of them should have the same size. However, as more tasks are learned, the chance that the space of a new task will overlap with the already used input space increases. Figure 4 shows the singular value spectra and quota of the input and hidden layer conceptors every time after a new task is learned. As the incremental learning proceeds, it becomes less likely for a new task to be in the free space. For example, the second task increases the quota of the input layer memory space by 0.1, whereas the 10th task increases it by only 0.03. However, CAB still manages to make the network learn new tasks based on their input components in the non-overlapping space.

(a) Singular value spectra of conceptors A (0) j on the input layer.

(b) Singular value spectra of conceptors A (1) j on the hidden layer.

Figure 4 : The development of singular value spectra of conceptors for "used-up" space on the input layer and hidden layer during incremental learning of 10 permuted MNIST tasks. Quota of these conceptors are displayed in the legends.

4.2. DISJOINT MNIST EXPERIMENT:

We then applied CAB to categorize the disjoint MNIST datasets into 10 classes (Srivastava et al., 2013; Lee et al., 2017) . In this experiment, the original MNIST dataset is divided into two disjoint datasets with the first one consisting of data for the first five digits (0 to 4), and the second one of the remaining five digits (5 to 9). This task requires a network to learn these two datasets one after the other, then examines its performance of classifying the entire MNIST testing images into 10 classes. The current state-of-the-art accuracy on this task, averaged over 10 learning trials, is 94.12(±0.27)%, achieved by Lee et al. (2017) using IMM. They also tested EWC on the same task and the average accuracy was 52.72(±1.36)%.

To test our method, we trained a feed-forward network with [784-800-10] neurons. Logistic sigmoid nonlinearities were used in both hidden and output layers, and the network was trained with vanilla SGD to minimize mean squared errors. The aperture α = 9 was used for all conceptors on all layers, learning rate η and regularization coefficient λ were chosen to be 0.1 and 0.005 respectively. The accuracy of CAB on this task, measured by repeating the experiment 10 times, is 94.91(±0.30)%.

It is worth mentioning that the network used by Lee et al. (2017) for testing IMM and EWC had [784-800-800-10] rectified linear units (ReLU), so CAB achieved better performance with fewer layers and neurons.

4.3. COMPUTATIONAL COST:

If a conceptor is computed by ridge regression, the time complexity is O(nN 2 + N 3 ) when the design matrix is dense, where n is the number of samples and N the number of features. In terms of wall time measures, the time taken to compute a conceptor from the entire MNIST training set (in this case, n = 55000 images and N = 784 pixels, corresponding to the input layer in our networks) is 0.42 seconds of standard notebook CPU time on average. Although we did not implement it in these experiments, incremental online adaptation of conceptors by gradient descent is also possible in principle and would come at a cost of O(N 2 ) per update.

5. CONCLUSION:

In this work, we first reviewed the conceptor-based incremental ridge regression algorithm, introduced in section 3.11 of Jaeger (2014) for memory management in recurrent neural networks. Then we derived its stochastic gradient descent version for optimizing the same objective. Finally we designed a conceptor-aided backprop algorithm by applying a conceptor to every linear layer of a feed-forward network. This method uses conceptors to guide gradients of parameters during the backpropagation procedure. As a result, learning a new task interferes only minimally with previously learned tasks, and the amount of already used network capacity can be monitored via the singular value spectra and quota of conceptors.

In Jaeger (2014) , different scenarios for continual learning are investigated in a reservoir computing setting. Two extreme cases are obtained when (i) the involved learning tasks are entirely unrelated to each other, versus (ii) all tasks come from the same parametric family of learning tasks. The two cases differ conspicuously with regards to the geometry of involved conceptors, and with regards to opportunities to re-use previously acquired functionality in subsequent learning episodes. The permuted MNIST task is an example of (i) while the disjoint MNIST task rather is of type (ii).

Conceptors provide an analytical tool to discuss the "family relatedness" and enabling/disabling conditions for continual learning in geometrical terms. Ongoing and future research is devoted to a comprehensive mathematical analysis of these phenomena which in our view lie at the heart of understanding continual learning.

:

Jaeger (2014) proposes a number of methods for analytical aperture optimization. It remains for future work to determine how these methods transfer to our situation.

Figures and Tables:

Figure 1: 3D point clouds (black dots) and their corresponding conceptors, represented by ellipsoids whose axes are the singular vectors of conceptors and the lengths of these axes match the singular values of conceptors. Each edge of the plot boxes range from -1 to +1 admitted by neural dynamics with a tanh nonlinearity; conceptor ellipsiods lie inside the unit sphere.

Figure 3: Average performance across already learned permuted MNIST tasks using CAB or EWC

Peer review (iclr2018):
 This paper introduces a method for learning new tasks, without interfering previous tasks, using conceptors. This method originates from linear algebra, where a the network tries to algebraically infer the main subspace where previous tasks were learned, and make the network learn the new task in a new sub-space which is "unused" until the present task in hand. The paper starts with describing the method and giving some context for the method and previous methods that deal with the same problem. In Section 2 the authors review conceptors. This method is algebraic method closely related to spanning sub spaces and SVD. The main advantage of using conceptors is their trait of boolean logics: i.e., their ability to be added and multiplied naturally. In section 3 the authors elaborate on reviewed ocnceptors method and show how to adapt this algorithm to SGD with back-propagation. The authors provide a version with batch SGD as well. In Section 4, the authors show their method on permuted MNIST. They compare the method to EWC with the same architecture. They show that their method more efficiently suffers on permuted MNIST from less degradation. Also, they compared the method to EWC and IMM on disjoint MNIST and again got the best performance. In general, unlike what the authors suggest, I do not believe this method is how biological agents perform their tasks in real life. Nevertheless, the authors show that their method indeed reduce the interference generated by a new task on the old learned tasks. I think that this work might interest the community since such methods might be part of the tools that practitioners have in order to cope with learning new tasks without destroying the previous ones. What is missing is the following: I think that without any additional effort, a network can learn a new task in parallel to other task, or some other techniques may be used which are not bound to any algebraic methods. Therefore, my only concern is that in this comparison the work bounded to very specific group of methods, and the question of what is the best method for continual learning remained open.