## Toward Understanding Catastrophic Forgetting in Continual Learning

Cuong V. Nguyen , Alessandro Achille , Michael Lam , Tal Hassner , † † † ‡∗ Vijay Mahadevan , Stefano Soatto † †

## † Amazon Web Services

{ nguycuo,aachille,michlam,vmahad,soattos @amazon.com }

## ‡ Facebook Inc.

{ talhassner @gmail.com }

## Abstract

We study the relationship between catastrophic forgetting and properties of task sequences. In particular, given a sequence of tasks, we would like to understand which properties of this sequence influence the error rates of continual learning algorithms trained on the sequence. To this end, we propose a new procedure that makes use of recent developments in task space modeling as well as correlation analysis to specify and analyze the properties we are interested in. As an application, we apply our procedure to study two properties of a task sequence: (1) total complexity and (2) sequential heterogeneity. We show that error rates are strongly and positively correlated to a task sequence's total complexity for some state-of-the-art algorithms. We also show that, surprisingly, the error rates have no or even negative correlations in some cases to sequential heterogeneity. Our findings suggest directions for improving continual learning benchmarks and methods.

## 1. Introduction

Continual learning (or life-long learning) [32, 37, 44] is the ability of a machine learning model to continuously learn from a stream of data, which could possibly be non-iid or come from different but related tasks. A continual learning system is required to adapt its current model to the new tasks or datasets without revisiting the previous data. Such a system should be able to positively transfer its current knowledge (summarized in its model) to the new tasks using as few data as possible, to avoid catastrophically forgetting the old tasks, and to transfer back its knowledge from new tasks to old tasks in order to improve overall performance.

In recent years, interest in continual learning has risen [1, 21, 26, 27, 29, 39, 40, 42, 51], especially from the deep

∗ Work done at Amazon.

learning research community, due to its potential to reduce training time and training set sizes (e.g., by continuously adapting from previous models), both of which are critical to the training of modern deep neural networks. Solving continual learning is also an essential step toward artificial general intelligence as it allows machines to continuously adapt to changes in the environment with minimal human intervention, a process analogous to human learning.

However, continual learning by deep models has proven to be very challenging due to catastrophic forgetting , a long known problem of training deep neural networks [4, 5, 14, 15, 28, 31, 34]. Catastrophic forgetting refers to the tendency of a model to forget all its previously learned tasks if not trained properly on a new task, e.g., when fine-tuning on the new task for a long time without proper regularization to the previous model parameters.

Recent work attempted to tackle this problem either by better training algorithms [21, 25, 33, 51], structure sharing [36, 41, 49], episodic memory [9, 27, 29], machinegenerated pseudo-data [18, 26, 42], or a combination of these approaches [29, 39]. Benchmarks to compare these methods typically constructed a sequence of tasks and then measured the algorithms' performance when transferring from one task to another. Two popular examples of these benchmarks are the permuted MNIST [15] and split MNIST [51].

In this paper, we seek to understand catastrophic forgetting at a more fundamental level. Specifically, we investigate the following question:

Given a sequence of tasks, which properties of the tasks influence the hardness of the entire sequence?

We measure task sequence hardness by the final error rate of a model trained sequentially on the tasks in the sequence.

An answer to this question is useful for continual learning research in several ways. First, it helps us estimate the hardness of a benchmark based on its individual tasks, thereby potentially assisting the development of new and bet-

ter benchmarks for continual learning. Additionally, knowing the hardness of a task sequence allows us to estimate a priori the cost and limits of running continual learning algorithms on it. Crucially, by gaining a better understanding of catastrophic forgetting at a more fundamental level, we gain more insights to develop better methods to mitigate it.

This work is the first attempt to answer the above question. We propose a new and general procedure that can be applied to study the relationship between catastrophic forgetting and properties of task sequences. Our procedure makes use of recent developments in task space modeling methods, such as the Task2Vec framework [2], to specify the interested properties. Then, we apply correlation analysis to study the relationship between the specified properties and the actual measures of catastrophic forgetting.

As an application, we use our procedure to analyze two properties of a task sequencetotal complexity and sequential heterogeneity -and design experiments to study their correlations with the sequence's actual hardness. We refer to total complexity as the total hardness of individual tasks in the sequence, while sequential heterogeneity measures the total dissimilarity between pairs of consecutive tasks.

We show how these two properties are estimated using the Task2Vec framework [2], which maps datasets (or equivalently, tasks) to vectors on a vector space. We choose these two properties for our analysis because of their intuitive relationships to the hardness of task sequences: since continual learning algorithms attempt to transfer knowledge from one task to another, both the hardness of each individual task and the dissimilarity between them should play a role in determining the effectiveness of the transfer.

The findings from our analysis are summarized below.

- · Total complexity has a strong correlation with the task sequence hardness measured by the actual error rate.
- · Sequential heterogeneity has little or no correlation with the task sequence hardness. When factoring out the task complexity, we even find negative correlations in some cases.

The first finding, although expected, emphasizes that we should take into account the complexity of each task when designing new algorithms or benchmarks, which is currently lacking in continual learning research. Besides, the research community is currently somewhat divided on the issue whether task similarity helps or hurts continual learning performance. Some authors showed that task similarity helps improve performance in the context of transfer learning [2, 3, 35], while some others conjectured that task dissimilarity could help improve continual learning performance [13]. Our second finding gives evidence that supports the latter view.

Deeper analysis into these phenomena suggests that (a) the task sequence hardness also depends on the ability to backward transfer (i.e., learning a new task helps a previous task) and (b) continual learning algorithms should be customized for specific task pairs to improve their effectiveness. We give detailed analysis and discussions in Sec. 7.

## 2. Continual learning algorithms and existing benchmarks

We overview modern continual learning algorithms and existing benchmarks used to evaluate them. For more comprehensive reviews of continual learning, we refer to Chen et al. [11] and Parisi et al. [30]

## 2.1. Continual learning algorithms

The simplest and most common approaches to continual learning use weight regularization to prevent catastrophic forgetting. Weight regularization adds a regularizer to the likelihood during training to pull the new weights toward the previous weights. It has been improved and applied to continual learning of deep networks in the elastic weight consolidation (EWC) algorithm [21], where the regularizer is scaled by the diagonal of the Fisher information matrix computed from the previous task. Since the diagonal Fisher information approximates the average Hessian of the likelihoods, EWC is closely related to Laplace propagation [17, 43], where Laplace's approximation is applied after each task to compute the regularizers. Besides Fisher information, the path integral of the gradient vector field along the parameter optimization trajectory can also be used for the regularizer, as in the synaptic intelligence (SI) approach [51].

Another form of regularization naturally arises by using Bayesian methods. For instance, variational continual learning (VCL) [29, 45] applied a sequence of variational approximations to the true posterior and used the current approximate posterior as prior for the new task. The KullbackLeibler term in the variational lower bound of VCL naturally regularizes the approximate posterior toward the prior. Improved training procedures have also been developed for this type of approximate Bayesian continual learning through the use of natural gradients [10, 47], fixed-point updates [52], and local approximation [6]. More expressive classes of variational distributions were also considered, including channel factorized Gaussian [22], multiplicative normalizing flow [22], or structured Laplace approximations [33].

The above methods can be complemented by an episodic memory, sometimes called a coreset , which stores a subset of previous data. Several algorithms have been developed for utilizing coresets, including gradient episodic memory (GEM) [27], averaged GEM [9], coreset VCL [29], and Stein gradients coreset [10].

Other algorithmic ideas to prevent catastrophic forgetting include moment matching [25], learning without forgetting [26], and deep generative replay [13, 19, 42]. Structure

sharing [36, 39] is also another promising direction that can be combined with the above algorithmic solutions to improve continual learning.

## 2.2. Existing benchmarks

The most common benchmarks for continual learning use MNIST [24] as the base dataset and construct various task sequences for continual learning. For example, permuted MNIST [15] applies a fixed random permutation on the pixels of MNIST input images for each task, creating a sequence of tasks that keep the original labels but have different input structures. Split MNIST [51], on the other hand, considers five consecutive binary classification tasks based on MNIST: 0/1, 2/3, . . . , 8/9. Another variant is rotated MNIST [27], where the digits are rotated by a fixed angle between 0 and 180 degrees in each task. Similar constructions can also be applied to the not-MNIST set [7], the fashion MNIST set [48], or the CIFAR set [23] such as in the split not-MNIST [29] and split CIFAR benchmarks [51].

Other continual learning benchmarks include ones typically used for reinforcement learning. For instance, Kirkpatrick et al. [21] tested the performance of EWC when learning to play Atari games. Schwarz et al. [38] proposed a new benchmark for continual learning based on the StarCraft II video game, where an agent must master a sequence of skills without forgetting the previously acquired skills.

## 3. Analysis of catastrophic forgetting

Recent developments in task space modeling, such as Task2Vec [2] and Taskonomy [50], provide excellent tools to specify and analyze relationships between different tasks from data. In this paper, we propose a novel and general procedure that utilizes these tools to study catastrophic forgetting. Our procedure is conceptually simple and can be summarized in the following steps:

- 1. Specify the properties of a task sequence that we are interested in and estimate these properties using a suitable task space modeling methods.
- 2. Estimate actual measures of catastrophic forgetting from real experiments. In our case, we measure catastrophic forgetting by the task sequence hardness, defined as the final error rate of a model trained sequentially on the sequence.
- 3. Use correlation analysis to study the correlations between the estimated properties in Step 1 and the actual measures in Step 2.

This procedure can be used even in other cases, such as transfer or multi-task learning, to study properties of new algorithms. For the rest of this paper, we demonstrate its use for analyzing two properties of task sequences and their effects on continual learning algorithms.

## 4. Total complexity and sequential heterogeneity of task sequences

We define two properties that we would like to investigate: the total complexity and sequential heterogeneity of a task sequence, and detail the methodology used to estimate these quantities from data. We start by first introducing the Task2Vec framework [2], the main tool that we employ to quantify the above properties.

## 4.1. Preliminaries: Task2Vec

Task2Vec [2] is a recently developed framework for embedding visual classification tasks as vectors in a real vector space. The embeddings have many desirable properties that allow reasoning about the semantic and taxonomic relations between different visual tasks. This is one of several recent attempts to provide tools for understanding the structure of task space. Other related efforts that can be used as alternatives to Task2Vec include, e.g., [12, 46, 50].

Given a labeled classification dataset, D = ( { x , y i i ) } N i =1 , Task2Vec works as follows. First, a network pre-trained on a large dataset (e.g., ImageNet), called the probe network, is applied to all the images x i in the dataset to extract the features from the last hidden layer (i.e., the value vectors returned by this layer). Using these features as new inputs and the labels y i , we then train the classification layer for the task. After the training, we compute the Fisher information matrix for the feature extractor parameters. Since the Fisher information matrix is very large for deep networks, in practice we usually approximate it by (1) using only the diagonal entries and (2) averaging the Fisher information of all weights in the same filter. This results in a vector representation with size equal to the number of filters in the probe network. In this paper, we will use a ResNet [16] probe network that only has convolutional layers.

Task2Vec embeddings have many properties that can be used to study the relationships between tasks. We discuss two properties that are most relevant to our work. The first of these properties is that the norms of the embeddings encode the difficulty of the tasks. This property can be explained intuitively by noticing that easy examples (those that the model is very confident about) have less contributions to the Fisher information while uncertain examples (those that are near the decision boundary) have more contributions. Hence, if the task is difficult, the model would be uncertain on many examples leading to a large embedding.

The second property that we are interested in is that Task2Vec embeddings can encode the similarity between tasks. Achille et al. [2] empirically showed this effect on the iNaturalist dataset [53], where the distances between Task2Vec embeddings strongly agree with the distances between natural taxonomical orders, hinting that the dissimilarity between tasks can be approximated from the distance

between them in the embedding space. The embeddings were also shown to be useful for model selection between different domains and tasks.

## 4.2. Total complexity

We now discuss the notions of total complexity and sequential heterogeneity of task sequences, and how we can estimate them from Task2Vec embeddings. We note that these definitions only capture specific aspects of sequence complexity and heterogeneity; however, they are enough to serve the purpose of our paper. In future work, we will consider more sophisticated definitions of sequence complexity and heterogeneity.

We define the total complexity of a task sequence as the sum of the complexities of its individual tasks. Formally, let T = ( t 1 , t 2 , . . . , t k ) be a sequence of k distinct tasks and C t ( ) be a function measuring the complexity of a task t . The total complexity of the task sequence T is:

$$C ( T ) = \sum _ { i = 1 } ^ { k } C ( t _ { i } ). \quad \ \ \ ( 1 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

We slightly abuse notation by using the same function C ( ) · for the complexity of both sequences and tasks.

For simplicity, we only consider sequences of distinct tasks where data for each task are only observed once. The scenario where data for one task may be observed many times requires different definitions of total complexity and sequential heterogeneity. We will leave this extension to future work.

A simple way to estimate the complexity C t ( ) of a task t is to measure the error rate of a model trained for this task. However, this method often gives unreliable estimates since it depends on various factors such as the choice of model and the training algorithm.

In this work, we propose to estimate C t ( ) from the Task2Vec embedding of task t . Specifically, we adopt the suggestion from Achille et al. [2] to measure the complexity of task t by its distance to the trivial task (i.e., the task embedded at the origin for standard Fisher embedding) in the embedding space. That is,

$$C ( t ) = d ( e _ { t }, e _ { 0 } ),$$

where e t and e 0 are the embeddings of task t and the trivial task respectively, and d ( · , · ) is a symmetric distance between two tasks in the embedding space. Following Achille et al. [2], we choose d ( · , · ) to be the normalized cosine distance:

$$d ( e _ { 1 }, e _ { 2 } ) = \cos \left ( \frac { e _ { 1 } } { e _ { 1 } + e _ { 2 } }, \frac { e _ { 2 } } { e _ { 1 } + e _ { 2 } } \right ), \quad \ ( 3 ) \quad \ \ i s n$$

where e 1 and e 2 are two task embeddings and the division is element-wise. This distance was shown to be well correlated with natural distances between tasks [2].

The total complexity in Eq. (1) depends on the sequence length. We can also consider the total complexity per task, C T /k ( ) , which does not depend on sequence length. In our analysis, however, we will only consider sequences of the same length. Hence, our results are not affected whether total complexity or total complexity per task is used.

We note that our total complexity measure is very crude and only captures some aspects of task sequence complexity. However, as we will show in Sec. 6.2, our measure is positively correlated with catastrophic forgetting and thus can be used to explain catastrophic forgetting. A possible future research direction would be to design better measures of task sequence complexity that can better explain catastrophic forgetting (i.e., by giving better correlation scores).

## 4.3. Sequential heterogeneity

We define the sequential heterogeneity of a task sequence as the sum of the dissimilarities between all pairs of consecutive tasks in the sequence. Formally, for a task sequence T = ( t 1 , t 2 , . . . , t k ) of distinct tasks, its sequential heterogeneity is:

$$\cdot \\ F ( T ) = \sum _ { i = 1 } ^ { k - 1 } F ( t _ { i }, t _ { i + 1 } ),$$

where F t, t ( ′ ) is a function measuring the dissimilarity between tasks t and t ′ . Note that we also use the same notation F ( ) · for sequential heterogeneity and task dissimilarity here, but its interpretation should be clear from the context.

The dissimilarity F t, t ( ′ ) can be naively estimated by applying transfer learning algorithms and measuring how well we can transfer between the two tasks. However, this would give a tautological measure of dissimilarity that is affected by both the model choice and the choice of the transfer learning algorithm.

To avoid this problem, we also propose to estimate F t, t ( ′ ) from the Task2Vec embedding. For our purpose, it is clear that we can use the distance d ( · , · ) of Eq. (3) as an estimate for F ( · , · ) . That is,

$$\overset { \kappa } { e } F ( t, t ^ { \prime } ) = d ( e _ { t }, e _ { t ^ { \prime } } ) = \cos \left ( \frac { e _ { t } } { e _ { t } + e _ { t ^ { \prime } } }, \frac { e _ { t ^ { \prime } } } { e _ { t } + e _ { t ^ { \prime } } } \right ).$$

The sequential heterogeneity in Eq. (4) only considers pairs of consecutive tasks, under the assumption that catastrophic forgetting is mostly influenced by the dissimilarity between these task pairs. In general, we can define other measures of heterogeneity, such as the total dissimilarity between all pairs of tasks. We will leave these extensions to future work.

Our choice of using Task2Vec to estimate C T ( ) and F T ( ) is more compatible with the multi-head models for continual learning [29, 51], which we will use in our experiments. In multi-head models, a separate last layer (the SoftMax layer) is trained for each different task and the other weights

are shared among tasks. This setting is consistent with the way Task2Vec is constructed in many cases. For instance, if we have two binary classification tasks whose labels are reversed, they would be considered similar by Task2Vec and are indeed very easy to transfer from one to another in the multi-head setting, by changing only the head.

## 5. Correlation analysis

Having defined total complexity and sequential heterogeneity, we now discuss how we can study their relationships to the hardness of a task sequence. Given a task sequence T = ( t 1 , t 2 , . . . , t k ) , we measure its actual hardness with respect to a continual learning algorithm A by the final error rate obtained after running A on the tasks t 1 , t 2 , . . . , t k sequentially. That is, the hardness of T with respect to A is:

$$H _ { A } ( T ) = \text{err} _ { A } ( T ). \text{ \quad \ \ } ( 6 ) \text{ \quad } \text{n}$$

In this paper, we choose final error rate as the measure of actual hardness as it is an important metric commonly used to evaluate continual learning algorithms. In future work, we will explore other metrics such as the forgetting rate [8].

To analyze the relationships between the hardness and total complexity or sequential heterogeneity, we employ correlation analysis as the main statistical tool. In particular, we sample M task sequences T , T 1 2 , . . . , T M and compute their hardness measures ( H A ( T i )) M i =1 as well as their total complexity ( C T ( i )) M i =1 and sequential heterogeneity ( F T ( i )) M i =1 measures. From these measures, we compute the Pearson correlation coefficients between hardness and total complexity measures or between hardness and sequential heterogeneity measures. These coefficients tell us how correlated these quantities are.

Formally, the Pearson correlation coefficient between two corresponding sets of measures X = ( X i ) M i =1 and Y = ( Y i ) M i =1 is defined as:

$$r _ { X Y } = \frac { \sum _ { i = 1 } ^ { M } ( X _ { i } - \bar { X } ) ( Y _ { i } - \bar { Y } ) } { \sqrt { \sum _ { i = 1 } ^ { M } ( X _ { i } - \bar { X } ) ^ { 2 } } \sqrt { \sum _ { i = 1 } ^ { M } ( Y _ { i } - \bar { Y } ) ^ { 2 } } }, \quad ( 7 ) \quad \text{$learn$} \\ \text{$the I$} \quad \text{and}$$

where ¯ X and ¯ Y are the means of X and Y respectively. In addition to the correlation coefficients, we can also compute the p-values, which tell us how statistically significant these correlations are.

When computing the correlations between the hardness measures H A ( T i ) and the total complexities C T ( i ) , it is often a good idea to constrain the task sequences T i to have the same length. The reason for this normalization is that longer sequences tend to have larger complexities, thus the correlation may be biased by the sequence lengths rather than reflecting the complexity of individual tasks.

Similarly, when computing the correlations between H A ( T i ) and the sequential heterogeneity F T ( i ) , it is also a good idea to constrain the total complexity of the task sequences T i to be the same, so that the individual tasks' complexities would not affect the correlations. This can be achieve by using the same set of individual tasks for all the sequences (i.e., the sequences are permutations of each other). We call the sequential heterogeneity obtained from this method the normalized sequential heterogeneity .

## 6. Experiments

We next describe the settings of our experiments and discuss our results. More detailed discussions on the implications of these results to continual learning and catastrophic forgetting research are provided in Sec. 7.

## 6.1. Settings

Datasets and task construction. We conduct our experiments on two datasets: MNIST and CIFAR-10, which are the most common datasets used to evaluate continual learning algorithms. For each of these sets, we construct a more general split version as follows. First, we consider all pairs of different labels as a unit binary classification task, resulting in a total of 45 unit tasks. From these unit tasks, we then create 120 task sequences of length five by randomly drawing, for each sequence, five unit tasks without replacement.

We also construct 120 split task sequences which are permutations of a fixed task set containing five random unit tasks to compute the normalized sequential heterogeneity. For each unit task, we train its Task2Vec embedding using a ResNet18 [16] probe network pre-trained on a combined dataset containing both MNIST and CIFAR-10.

Algorithms and network architectures. We choose two recent continual learning algorithms to analyze in our experiments: synaptic intelligence (SI) [51] and variational continual learning (VCL) [29]. For the experiments on MNIST, we also consider the coreset version of VCL (coreset VCL). These algorithms are among the state-of-the-art continual learning algorithms on the considered datasets, with SI representing the regularization-based methods, VCL representing the Bayesian methods, and coreset VCL combining Bayesian and rehearsal methods.

On CIFAR-10, we run SI with the same network architecture as those considered in [51]: a CNN with four convolutional layers, followed by two dense layers with dropout. Since VCL was not originally developed with convolutional layers, we flatten the input images and train with a fully connected network containing four hidden layers, each of which has 256 hidden units. On MNIST, we run both SI and VCL with a fully connected network containing two hidden layers, each of which has 256 hidden units. We denote this setting by MNIST256 2 .

Since MNIST is a relatively easy dataset, we may not observe meaningful results if all the errors obtained from different sequences are low and not very different. Thus,

to make the dataset harder for the learning algorithms, we also consider smaller network architectures. In particular, we consider fully connected networks with a single hidden layer, containing either 50 hidden units (for MNIST50 ) or 20 hidden units (for MNIST20 ). Following [29, 51], we also use the multi-head version of the models where a separate last layer (the SoftMax layer) is trained for each different task and the other weights are shared among tasks. For coreset VCL, we use random coresets with sizes 40, 40, 20 for MNIST256 2 , MNIST50 and MNIST20 respectively.

Optimizer settings. For both SI and VCL, we set the regularization strength parameter to the default value λ = 1 . In all of our experiments, the models are trained using Adam optimizer [20] with learning rate 0.001. Similar to [51], we set the batch size to be 256 in CIFAR-10 and 64 in MNIST settings for SI. We run this algorithm for 60, 10, 10, 5 epochs per task on CIFAR-10, MNIST256 2 , MNIST50 and MNIST20 respectively.

For VCL and coreset VCL, we set the batch size to be the training set size [29] and run the algorithms for 50, 120, 50, 20 epochs per task on CIFAR-10, MNIST256 2 , MNIST50 and MNIST20 respectively. For all algorithms, we run each setting ten times using different random seeds and average their errors to get the final error rates.

## 6.2. Results

Tables 1(a-c) show the correlation coefficients and their p-values obtained from our experiments for the total complexity, sequential heterogeneity, and normalized sequential heterogeneity, respectively. We also show the scatter plots of the errors versus these quantities, together with the linear regression fits for the CIFAR-10 dataset in Fig. 1. All plots in the experiments, including those for the MNIST dataset, are provided in Fig. 4, 5, and 6.

Table 1(a) and Fig. 1(a) show strong positive correlations between error rate and total complexity for both SI and VCL in the CIFAR-10 setting, with a correlation coefficient of 0.86 for the former algorithm and 0.69 for the latter. These correlations are both statistically significant with p-values less than 0.01. On the MNIST256 2 settings, SI and coreset VCL have weak positive correlations with total complexity, where the algorithms have correlation coefficients of 0.24 and 0.28, both with p-values less than 0.01, respectively.

When we reduce the capacity of the network and make the problem relatively harder (i.e., in the MNIST50 and MNIST20 settings), we observe stronger correlations for all three algorithms. With the smallest network (in MNIST20 ), all the algorithms have statistically significant positive correlation with total complexity.

In terms of sequential heterogeneity, Table 1(b) and Fig. 1(b) show that it has a weak positive correlation with error rate in the CIFAR-10 setting. In particular, SI and VCL have correlation coefficients of 0.30 and 0.21 (both statistically significant), respectively. Interestingly, we find no significant correlation between error rate and sequential heterogeneity in all the MNIST settings, which suggests that heterogeneity may not be a significant factor determining the performance of continual learning algorithms on this dataset.

Since the complexity of each individual task in a sequence may influence the heterogeneity between the tasks (e.g., an easy task may be more similar to another easy task than to a hard task), the complexity may indirectly affect the results in Table 1(b). To avoid this problem, we also look at the normalized sequential heterogeneity in Table 1(c) and Fig 1(c), where the set of tasks is fixed and thus task complexity has been factored out.

Surprisingly, Table 1(c) reports some negative correlations between error rate and sequential heterogeneity. For example, the correlation coefficient for SI on CIFAR-10 is -0.25 with a p-value less than 0.01, while there is no significant correlation for this algorithm on the MNIST dataset. VCL, on the other hand, has negative correlations with coefficients -0.20 and -0.21, respectively on MNIST50 and MNIST20 , with p-values less than 0.05. Coreset VCL also has negative correlation between its error rate and sequential heterogeneity on MNIST50 , with coefficient -0.26 and pvalue less than 0.01. These unexpected results suggest that in some cases, dissimilarity between tasks may even help continual learning algorithms, a fact contrary to the common assumption that the performance of continual learning algorithms would degrade if the tasks they need to solve are very different [3, 35].

## 7. Discussions

On total complexity. The strong positive correlations between error rate and total complexity found in our analysis show that task complexity is an important factor in determining the effectiveness of continual learning algorithms. However, this factor is usually not taken into consideration when designing new algorithms or benchmarks. We suggest that task complexity is explicitly considered to improve algorithm and benchmark design. For example, different transfer methods can be used depending on whether one transfers from an easy task to a hard one or vice versa, rather than using a single transfer technique across all task complexities, as currently done in the literature. Similarly, when designing new benchmarks for continual learning, it is also useful to provide different complexity structures to test the effectiveness of continual learning algorithms on a broader range of scenarios and difficulty levels.

To illustrate the usefulness of comparing on various benchmarks, we construct two split MNIST sequences, one of which has high total complexity while the other has low total complexity. The sequences are constructed by starting with the binary classification task 0/1 and greedily adding tasks that have the highest (or lowest) complexity C t ( ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-516895e9-05da-4f1f-b06e-c5242249a27b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-9896bd2d-7524-48e6-9e78-f4b535d31b63.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-3e3eb63e-aab8-4857-997c-b39ad5c345bb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-e0eb69f5-af1a-450f-bc0c-4e120fa435ac.jpg)

(b) Sequential heterogeneity

(c) Normalized sequential heterogeneity

Figure 1. Error vs. (a) total complexity, (b) sequential heterogeneity and (c) normalized sequential heterogeneity on CIFAR-10 , together with the linear regression fits and 95% confidence intervals. Green (red) color indicates statistically significant positive (negative) correlations. Black color indicates negligible correlations.

|     | Variable                            | Algorithm          | MNIST- 256 2                                                  | MNIST- 50                                                      | MNIST- 20                                                     | CIFAR-10                                    |
|-----|-------------------------------------|--------------------|---------------------------------------------------------------|----------------------------------------------------------------|---------------------------------------------------------------|---------------------------------------------|
| (a) | Total Complexity                    | SI VCL Coreset VCL | 0.24 ( p < 0 . 01 ) 0.05 ( p = 0 . 59 ) 0.28 ( p < 0 . 01 )   | 0.22 ( p < 0 . 05 ) 0.17 ( p = 0 . 07 ) 0.41 ( p < 0 . 01 )    | 0.36 ( p < 0 . 01 ) 0.21 ( p < 0 . 05 ) 0.37 ( p < 0 . 01 )   | 0.86 ( p < 0 . 01 ) 0.69 ( p < 0 . 01 ) -   |
| (b) | Sequential Heterogeneity            | SI VCL Coreset VCL | -0.01 ( p = 0 . 86 ) 0.04 ( p = 0 . 69 ) 0.09 ( p = 0 . 31 )  | 0.05 ( p = 0 . 55 ) 0.01 ( p = 0 . 88 ) 0.12 ( p = 0 . 18 )    | 0.07 ( p = 0 . 48 ) 0.05 ( p = 0 . 58 ) 0.18 ( p = 0 . 05 )   | 0.30 ( p < 0 . 01 ) 0.21 ( p < 0 . 05 ) -   |
| (c) | Normalized Sequential Heterogeneity | SI VCL Coreset VCL | -0.07 ( p = 0 . 43 ) 0.03 ( p = 0 . 76 ) -0.08 ( p = 0 . 37 ) | -0.04 ( p = 0 . 65 ) -0.20 ( p < 0 . 05 ) -0.26 ( p < 0 . 01 ) | 0.05 ( p = 0 . 58 ) -0.21 ( p < 0 . 05 ) -0.16 ( p = 0 . 07 ) | -0.25 ( p < 0 . 01 ) -0.17 ( p = 0 . 06 ) - |

Table 1. Correlation coefficients (p-values) between error rate and (a) total complexity, (b) sequential heterogeneity, and (c) normalized sequential heterogeneity of three state-of-the-art continual learning algorithms (SI, VCL, coreset VCL) on four different tests conducted with the CIFAR-10 and MNIST datasets. Results with statistical significance ( p &lt; 0 05 . ) are shown in bold.

Fig. 3 shows these sequences and the error rates of VCL, coreset VCL and SI when evaluated on them. We also show the error rates of the algorithms on the standard split MNIST sequence for comparison. From the figure, if we only compare on the standard sequence, we may conclude that coreset VCL and SI have the same performance. However, if we consider the other two sequences, we can see that SI is in fact slightly better than coreset VCL. This small experiment suggests that we should use various benchmarks, ideally with different levels of complexity, for better comparison of continual learning algorithms.

It is also worth noting that although the correlation between error rate and task complexity seems trivial, we are still not very clear which definition of task sequence complexity would be best to explain catastrophic forgetting (i.e., to give the best correlations). In this paper, we propose the first measure for this purpose, the total complexity.

Onsequential heterogeneity. The weak or negative correlations between error rate and sequential heterogeneity found in our analysis show an interesting contradiction to our intuition on the relationship between catastrophic forgetting and task dissimilarity. We emphasize that in our context, the weak and negative correlations are not a negative result, but actually a positive result. In fact, some previous work showed that task similarity helps improve performance in the context of transfer learning [2, 3, 35], while some others claimed that task dissimilarity could help continual learning [13] although their discussion was more related to the permuted MNIST setting. Our finding gives evidence that supports the latter view in the split MNIST and split CIFAR-10 settings.

To identify possible causes of this phenomenon, we carefully analyze the changes in error rates of VCL and SI on CIFAR-10 and observe some issues that may cause the neg-

Figure 2. Details of the error rates of VCL and SI on two typical task sequences from CIFAR-10. Each column shows the errors on a particular task when subsequent tasks are continuously observed. Sequence 1 contains the binary tasks 2/9, 0/4, 3/9, 4/8, 1/2 with sequential heterogeneity 0.091, while sequence 2 contains the tasks 1/2, 2/9, 3/9, 0/4, 4/8 with sequential heterogeneity 0.068 (the labels are encoded to 0, 1, . . . , 9 as usually done for this dataset). For both algorithms, the final average errors (the last points in the right-most plots) on sequence 2 are higher than those on sequence 1, despite sequence 1's higher sequential heterogeneity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-c2c5467e-d481-4fb0-bf3e-b721c46c68b0.jpg)

Figure 3. Average error rates of VCL, coreset VCL and SI on 3 task sequences from MNIST with different complexity levels. The high complexity sequence contains the binary tasks 0/1, 2/5, 3/5, 2/3, 2/6 with total complexity 0.48, while the low complexity sequence contains the tasks 0/1, 1/8, 1/3, 1/5, 7/8 with total complexity 0.35. The standard sequence contains the common split 0/1, 2/3, 4/5, 6/7, 8/9 with total complexity 0.41.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-5bb4ec3e-5600-4bc0-9a8f-c00a9a2821c3.jpg)

but the algorithm tends not to transfer very well on the second sequence. This inability to transfer leads to higher error rates on tasks 3, 4, and 5 even when the algorithm learns them for the first time. One possible cause of this problem could be that a fixed regularization strength λ = 1 is used for all tasks, making the algorithm unable to adapt to new tasks well. This explanation suggests that we should customize the algorithm (e.g., by tuning the λ values or the optimizer) for effectively transferring between different pairs of tasks in the sequence.

ative correlations. For illustration, we show in Fig. 2 the detailed error rates of these algorithms on two typical task sequences where the final average error rates do not conform with the sequential heterogeneity. Both of these sequences have the same total complexity, with the first sequence having higher sequential heterogeneity.

From the changes in error rates of VCL in Fig. 2, we observe that for the first sequence, learning a new task would cause forgetting of its immediate predecessor task but could also help a task learned before that. For instance, learning task 3 and task 5 increases the errors on task 2 and task 4 respectively, but helps reduce errors on task 1 (i.e., backward transferring to task 1). This observation suggests that the dissimilarities between only consecutive tasks may not be enough to explain catastrophic forgetting, and thus we should take into account the dissimilarities between a task and all the previously learned tasks.

From the error rates of SI in Fig. 2, we observe a different situation. In this case, catastrophic forgetting is not severe,

Future directions. The analysis offered by our paper provides a general and novel methodology to study the relationship between catastrophic forgetting and properties of task sequences. Although the two measures considered in our paper, total complexity and sequential heterogeneity, can explain some aspects of catastrophic forgetting, the correlations in Table 1 are not very strong (i.e., their coefficients are not near 1 or -1). Thus, they can still be improved to provide better explanations for the phenomenon. Besides these two measures, we can also design other measures for properties such as intransigence [8].

## 8. Conclusion

This paper developed a new analysis for studying relationships between catastrophic forgetting and properties of task sequences. An application of our analysis to two simple properties suggested that task complexity should be considered when designing new continual learning algorithms or benchmarks, and continual learning algorithms should be customized for specific transfers. Our analysis can be extended to study other relationships between algorithms and task structures such as the effectiveness of transfer or multi-task learning with respect to properties of tasks.

Figure 4. Total complexity vs. average error , together with the linear regression fit and 95% confidence interval, for each algorithm and test in Table 1(a). Green color indicates statistically significant positive correlations. Black color indicates negligible correlations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-89e5efdf-3d4e-4285-80d4-7679ddf108fe.jpg)

Figure 5. Sequential heterogeneity vs. average error , together with the linear regression fit and 95% confidence interval, for each algorithm and test in Table 1(b). Green color indicates statistically significant positive correlations. Black color indicates negligible correlations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-227a2fe4-9ae9-4a0d-a5bf-c42f2a70d82c.jpg)

Figure 6. Normalized sequential heterogeneity vs. average error , together with the linear regression fit and 95% confidence interval, for each algorithm and test in Table 1(c). Red color indicates statistically significant negative correlations. Black color indicates negligible correlations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01091v1-figure-9d16d42f-90b7-4d20-ae4d-c45be42b1b2c.jpg)

## References

- [1] A. Achille, T. Eccles, L. Matthey, C. Burgess, N. Watters, A. Lerchner, and I. Higgins. Life-long disentangled representation learning with cross-domain latent homologies. In Advances in Neural Information Processing Systems , pages 9895-9905, 2018. 1
- [2] A. Achille, M. Lam, R. Tewari, A. Ravichandran, S. Maji, C. Fowlkes, S. Soatto, and P. Perona. Task2Vec: Task embedding for meta-learning. arXiv:1902.03545 , 2019. 2, 3, 4, 7
- [3] H. B. Ammar, E. Eaton, M. E. Taylor, D. C. Mocanu, K. Driessens, G. Weiss, and K. Tuyls. An automated measure of MDP similarity for transfer in reinforcement learning. In AAAI Conference on Artificial Intelligence Workshops , 2014. 2, 6, 7
- [4] B. Ans and S. Rousset. Avoiding catastrophic forgetting by coupling two reverberating neural networks. Comptes Rendus de l'Acad´ emie des Sciences-Series III-Sciences de la Vie , 320(12):989-997, 1997. 1
- [5] B. Ans and S. Rousset. Neural networks with a self-refreshing memory: Knowledge transfer in sequential learning tasks without catastrophic forgetting. Connection Science , 12(1):119, 2000. 1
- [6] T. D. Bui, C. V. Nguyen, S. Swaroop, and R. E. Turner. Partitioned variational inference: A unified framework encompassing federated and continual learning. arXiv:1811.11206 , 2018. 2
- [7] Y. Bulatov. notMNIST dataset. http://yaroslavvb. blogspot.com/2011/09/notmnist-dataset. html , 2011. Accessed: 2019-01-10. 3
- [8] A. Chaudhry, P. K. Dokania, T. Ajanthan, and P. H. Torr. Riemannian walk for incremental learning: Understanding forgetting and intransigence. In European Conference on Computer Vision , 2018. 5, 8
- [9] A. Chaudhry, M. Ranzato, M. Rohrbach, and M. Elhoseiny. Efficient lifelong learning with A-GEM. In International Conference on Learning Representations , 2019. 1, 2
- [10] Y. Chen, T. Diethe, and N. Lawrence. Facilitating Bayesian continual learning by natural gradients and Stein gradients. In Continual Learning Workshop @ NeurIPS , 2018. 2
- [11] Z. Chen and B. Liu. Lifelong Machine Learning . Morgan &amp; Claypool Publishers, 2016. 2
- [12] H. Edwards and A. Storkey. Towards a neural statistician. In International Conference on Learning Representations , 2017. 3
- [13] S. Farquhar and Y. Gal. Towards robust evaluations of continual learning. arXiv:1805.09733 , 2018. 2, 7
- [14] R. M. French. Catastrophic forgetting in connectionist networks. Trends in cognitive sciences , 3(4):128-135, 1999. 1
- [15] I. J. Goodfellow, M. Mirza, D. Xiao, A. Courville, and Y. Bengio. An empirical investigation of catastrophic forgetting in gradient-based neural networks. In International Conference on Learning Representations , 2014. 1, 3

- [16] K. He, X. Zhang, S. Ren, and J. Sun. Deep residual learning for image recognition. In IEEE Conference on Computer Vision and Pattern Recognition , pages 770-778, 2016. 3, 5
- [17] F. Husz´ ar. Note on the quadratic penalties in elastic weight consolidation. Proceedings of the National Academy of Sciences , 115(11):E2496-E2497, 2018. 2
- [18] D. Isele and A. Cosgun. Selective experience replay for lifelong learning. In AAAI Conference on Artificial Intelligence , 2018. 1
- [19] N. Kamra, U. Gupta, and Y. Liu. Deep generative dual memory network for continual learning. arXiv:1710.10368 , 2017. 2
- [20] D. P. Kingma and J. Ba. Adam: A method for stochastic optimization. In International Conference on Learning Representations , 2015. 6
- [21] J. Kirkpatrick, R. Pascanu, N. Rabinowitz, J. Veness, G. Desjardins, A. A. Rusu, K. Milan, J. Quan, T. Ramalho, A. Grabska-Barwinska, D. Hassabis, C. Clopath, D. Kumaran, and R. Hadsell. Overcoming catastrophic forgetting in neural networks. Proceedings of the National Academy of Sciences , 2017. 1, 2, 3
- [22] M. Kochurov, T. Garipov, D. Podoprikhin, D. Molchanov, A. Ashukha, and D. Vetrov. Bayesian incremental learning for deep neural networks. In International Conference on Learning Representations Workshop , 2018. 2
- [23] A. Krizhevsky and G. Hinton. Learning multiple layers of features from tiny images. Technical report, University of Toronto, 2009. 3
- [24] Y. LeCun, C. Cortes, and C. Burges. MNIST handwritten digit database. AT&amp;T Labs [Online]. Available: http://yann. lecun. com/exdb/mnist , 2, 2010. 3
- [25] S.-W. Lee, J.-H. Kim, J. Jun, J.-W. Ha, and B.-T. Zhang. Overcoming catastrophic forgetting by incremental moment matching. In Advances in Neural Information Processing Systems , pages 4652-4662, 2017. 1, 2
- [26] Z. Li and D. Hoiem. Learning without forgetting. IEEE Transactions on Pattern Analysis and Machine Intelligence , 40(12):2935-2947, 2018. 1, 2
- [27] D. Lopez-Paz and M. Ranzato. Gradient episodic memory for continual learning. In Advances in Neural Information Processing Systems , pages 6467-6476, 2017. 1, 2, 3
- [28] M. McCloskey and N. J. Cohen. Catastrophic interference in connectionist networks: The sequential learning problem. Psychology of Learning and Motivation , 24:109-165, 1989. 1
- [29] C. V. Nguyen, Y. Li, T. D. Bui, and R. E. Turner. Variational continual learning. In International Conference on Learning Representations , 2018. 1, 2, 3, 4, 5, 6
- [30] G. I. Parisi, R. Kemker, J. L. Part, C. Kanan, and S. Wermter. Continual lifelong learning with neural networks: A review. arXiv:1802.07569 , 2018. 2
- [31] R. Ratcliff. Connectionist models of recognition memory: Constraints imposed by learning and forgetting functions. Psychological Review , 97(2):285, 1990. 1
- [32] M. B. Ring. CHILD: A first step towards continual learning. Machine Learning , 28(1):77-104, 1997. 1
- [33] H. Ritter, A. Botev, and D. Barber. Online structured Laplace approximations for overcoming catastrophic forgetting. arXiv:1805.07810 , 2018. 1, 2
- [34] A. Robins. Catastrophic forgetting, rehearsal and pseudorehearsal. Connection Science , 7(2):123-146, 1995. 1
- [35] S. Ruder and B. Plank. Learning to select data for transfer learning with Bayesian optimization. In Conference on Empirical Methods in Natural Language Processing , pages 372-382, 2017. 2, 6, 7
- [36] A. A. Rusu, N. C. Rabinowitz, G. Desjardins, H. Soyer, J. Kirkpatrick, K. Kavukcuoglu, R. Pascanu, and R. Hadsell. Progressive neural networks. arXiv:1606.04671 , 2016. 1, 3
- [37] J. C. Schlimmer and D. Fisher. A case study of incremental concept induction. In AAAI Conference on Artificial Intelligence , volume 86, pages 496-501, 1986. 1
- [38] J. Schwarz, D. Altman, A. Dudzik, O. Vinyals, Y. W. Teh, and R. Pascanu. Towards a natural benchmark for continual learning. In Continual Learning Workshop @ NeurIPS , 2018. 3
- [39] J. Schwarz, J. Luketina, W. M. Czarnecki, A. GrabskaBarwinska, Y. W. Teh, R. Pascanu, and R. Hadsell. Progress &amp;compress: A scalable framework for continual learning. In International Conference on Machine Learning , 2018. 1, 3
- [40] J. Serr`, D. Sur´ ıs, M. Miron, and A. Karatzoglou. Overcoming a catastrophic forgetting with hard attention to the task. In International Conference on Machine Learning , 2018. 1
- [41] A. Sharif Razavian, H. Azizpour, J. Sullivan, and S. Carlsson. CNN features off-the-shelf: An astounding baseline for recognition. In IEEE Conference on Computer Vision and Pattern Recognition Workshops , pages 806-813, 2014. 1
- [42] H. Shin, J. K. Lee, J. Kim, and J. Kim. Continual learning with deep generative replay. In Advances in Neural Information Processing Systems , pages 2990-2999, 2017. 1, 2
- [43] A. J. Smola, S. Vishwanathan, and E. Eskin. Laplace propagation. In Advances in Neural Information Processing Systems , pages 441-448, 2004. 2
- [44] R. S. Sutton and S. D. Whitehead. Online learning with random representations. In International Conference on Machine Learning , pages 314-321, 1993. 1
- [45] S. Swaroop, C. V. Nguyen, T. D. Bui, and R. E. Turner. Improving and understanding variational continual learning. In Continual Learning Workshop @ NeurIPS , 2018. 2
- [46] A. T. Tran, C. V. Nguyen, and T. Hassner. Transferability and hardness of supervised classification tasks. In IEEE/CVF International Conference on Computer Vision , 2019. 3
- [47] H. Tseran, M. E. Khan, T. Harada, and T. D. Bui. Natural variational continual learning. In Continual Learning Workshop @NeurIPS , 2018. 2
- [48] H. Xiao, K. Rasul, and R. Vollgraf. Fashion-MNIST: A novel image dataset for benchmarking machine learning algorithms. arXiv:1708.07747 , 2017. 3
- [49] J. Yoon, E. Yang, J. Lee, and S. J. Hwang. Lifelong learning with dynamically expandable networks. In International Conference on Learning Representations , 2018. 1
- [50] A. R. Zamir, A. Sax, W. Shen, L. Guibas, J. Malik, and S. Savarese. Taskonomy: Disentangling task transfer learn-

ing. In IEEE Conference on Computer Vision and Pattern Recognition , 2018. 3

- [51] F. Zenke, B. Poole, and S. Ganguli. Continual learning through synaptic intelligence. In International Conference on Machine Learning , pages 3987-3995, 2017. 1, 2, 3, 4, 5, 6
- [52] C. Zeno, I. Golan, E. Hoffer, and D. Soudry. Task agnostic continual learning using online variational Bayes. In Bayesian Deep Learning Workshop @ NeurIPS , 2018. 2
- [53] X. Zhang, Y. Cui, Y. Song, H. Adam, and S. Belongie. The iMaterialist challenge 2017 dataset. In FGVC Workshop at CVPR , volume 2, 2017. 3
Cuong V. Nguyen Alessandro Achille Michael Lam Tal Hassner Vijay Mahadevan Stefano Soatto
Toward Understanding Catastrophic Forgetting in Continual Learning
We study the relationship between catastrophic forgetting and properties of
task sequences. In particular, given a sequence of tasks, we would like to
understand which properties of this sequence influence the error rates of
continual learning algorithms trained on the sequence. To this end, we propose
a new procedure that makes use of recent developments in task space modeling as
well as correlation analysis to specify and analyze the properties we are
interested in. As an application, we apply our procedure to study two
properties of a task sequence: (1) total complexity and (2) sequential
heterogeneity. We show that error rates are strongly and positively correlated
to a task sequence's total complexity for some state-of-the-art algorithms. We
also show that, surprisingly, the error rates have no or even negative
correlations in some cases to sequential heterogeneity. Our findings suggest
directions for improving continual learning benchmarks and methods.