<?xml version="1.0" encoding="UTF-8"?>
<article xmlns:xlink="http://www.w3.org/1999/xlink">
<front><article-meta>
<title-group><article-title>Replay-Based Continual Learning on Sequential Digit Tasks: Buffer Size, Task Similarity, and the Effect of</article-title></title-group>
<contrib-group>
<contrib contrib-type="author"><name><surname>Marchetti</surname><given-names>Anna L.</given-names></name><xref ref-type="aff" rid="aff1"/></contrib>
<contrib contrib-type="author"><name><surname>Novak</surname><given-names>Jan H.</given-names></name><xref ref-type="aff" rid="aff2"/></contrib>
<contrib contrib-type="author"><name><surname>Tanabe</surname><given-names>Hiroshi</given-names></name><xref ref-type="aff" rid="aff3"/></contrib>
<aff id="aff1">Dipartimento di Informatica, Universita di Pisa, Pisa, Italy.</aff>
<aff id="aff2">Faculty of Information Technology, Czech Technical University in Prague, Praha, Czech Republic.</aff>
<aff id="aff3">RIKEN Center for Computational Science, Kobe, Japan.</aff>
</contrib-group>
<abstract><p>We study replay-based mitigation of forgetting on five binary digit tasks using the 1797-sample subset of UCI Optical Recognition of Handwritten Digits distributed by scikit-learn. A shared 64-16-16 ReLU network with task-specific sigmoid heads is trained under sequential learning, Experience Replay (ER), past-task-balanced Interleaved Replay (IR), and joint training with all tasks available from the outset. Reported means over ten seeds give Backward Transfer (BWT) -0.170 and final accuracy 0.861 for naive training, versus joint-reference accuracy 0.996. In the natural-layout ER condition, nominal per-class capacities 25 to 200 improve BWT from -0.035 to -0.010 and accuracy from 0.969 to 0.988; capacity 200 exceeds the available training examples per class and stores the full observed class data. IR does not significantly outperform ER in BWT at any capacity by the reported unadjusted tests (p ≥ 0.058). ER and IR also differ in current-task inclusion and expected past-task rehearsal dose, so this contrast does not isolate temporal interleaving alone. Per-task feature permutations reduce measured mean-image correlation and are associated with less forgetting in several exploratory contrasts. The twenty-six reported tests are unadjusted; applying Holm to their rounded p-values leaves none below 0.05. Results characterize one small task-incremental benchmark, and mechanism, execution, and reproducibility claims require the original code and per-run outputs.</p></abstract>
<kwd-group kwd-group-type="author"><kwd>Continual learning</kwd><kwd>Catastrophic forgetting</kwd><kwd>Experience replay</kwd><kwd>Complementary learning systems</kwd><kwd>Task</kwd></kwd-group>
</article-meta></front>
<body>
<sec id="sec-1">
<p>similarity, Optdigits.</p>
</sec>
<sec id="sec-2"><label>1</label><title>Introduction</title>
<p>Artificial learning systems are increasingly expected to acquire knowledge sequentially: new users appear, sensor calibrations drift, and novel categories must be learned while earlier ones remain in service. Networks trained with gradient descent handle this demand poorly. When the training distribution shifts from one task to the next, the updates that fit the new task partially overwrite the parameters that encoded the previous one, so that competence acquired earlier degrades rapidly. This failure mode, known as catastrophic interference, was demonstrated convincingly in the early connectionist literature: McCloskey and Cohen <xref ref-type="bibr" rid="ref-1">[1]</xref> showed that a network trained by back-propagation of errors <xref ref-type="bibr" rid="ref-2">[2]</xref> on one paired-associate list and then a second lost almost all retention of the first, and they argued that the problem was intrinsic to the way distributed representations are updated rather than an artifact of a particular architecture.</p>
<p>The most influential theoretical account of how biological brains sidestep this problem is the Complementary Learning Systems (CLS) framework of McClelland, McNaughton, and O&apos;Reilly <xref ref-type="bibr" rid="ref-3">[3]</xref>. CLS theory posits two interacting memories: a fast, sparse hippocampal store that binds the specifics of new episodes, and a slow neocortical system that acquires the statistical structure of the environment through gradual, interleaved re-exposure to old and new material. The computational translation of this account is rehearsal: retain some record of past experience and re-present it, interleaved with new data, while learning continues. Robins <xref ref-type="bibr" rid="ref-4">[4]</xref> demonstrated that rehearsal with a small store of true exemplars — and, more remarkably, pseudorehearsal with internally generated pseudo-exemplars — substantially reduces catastrophic forgetting in small networks. Rehearsal-style buffers have since become a standard tool, most visibly in reinforcement learning, where Experience Replay (ER) stabilizes learning from non-stationary streams <xref ref-type="bibr" rid="ref-5">[5]</xref>.</p>
<p>Contemporary continual learning research has broadened this toolkit considerably. Regularization approaches such as elastic weight consolidation identify parameters important to old tasks and penalize their modification <xref ref-type="bibr" rid="ref-6">[6]</xref>; architectural approaches allocate or grow task-specific capacity; and replay-based approaches store, condense, or generate past data and mix it into current training. Recent surveys <xref ref-type="bibr" rid="ref-7">[7]</xref>, <xref ref-type="bibr" rid="ref-8">[8]</xref> and overview chapters <xref ref-type="bibr" rid="ref-9">[9]</xref> catalogue dozens of variants and report that replay methods with even modest memory remain among the most robust baselines across benchmarks. At the same time, the empirical picture is fragmented: reported comparisons typically change several design factors at once — the amount of stored data, the rule by which stored samples are drawn into training, and the relatedness of the tasks themselves — so that the individual contribution of each factor is difficult to isolate.</p>
<p>The study uses the 1797-sample digit subset bundled by scikit-learn <xref ref-type="bibr" rid="ref-10">[10]</xref>, <xref ref-type="bibr" rid="ref-11">[11]</xref>, partitioned into five tasks discriminating digit pairs 2i and 2i+1. ER draws uniformly from a class-balanced buffer including the current task; IR uses the same buffer but samples only past tasks. Thus the compared distributions also change past-task rehearsal dose and current-task duplication. Joint training provides a future-data reference, not a mathematical upper bound. Per-task input permutations manipulate layout statistics, and ten seeded replications support conditional exploratory comparisons.</p>
<p>The design varies nominal buffer capacity, replay sampling distribution and input layout. Shared splits, initialization and buffer snapshots are reported controls, subject to code confirmation. The evidence concerns aggregate dose response and exploratory layout effects; it neither proves a pure scheduling effect nor establishes a biological memory mechanism.</p>
</sec>
<sec id="sec-3"><label>2</label><title>Related Work</title>
<p>Catastrophic interference was characterized long before it became a benchmarking concern. McCloskey and Cohen <xref ref-type="bibr" rid="ref-1">[1]</xref> trained back-propagation networks <xref ref-type="bibr" rid="ref-2">[2]</xref> on paired-associate lists and observed that learning a second list drove retention of the first toward zero, in sharp contrast to human subjects, who showed moderate proactive and retroactive interference but no wholesale erasure. Their analysis located the problem in the overlap of internal representations: because the same hidden units participate in encoding both lists, gradients from the new task move weights that the old task depends on. Interference is therefore expected to scale with representational overlap between sequentially learned tasks — an expectation that directly motivates our third hypothesis.</p>
<p>The CLS framework <xref ref-type="bibr" rid="ref-3">[3]</xref> motivates retaining episodes and re-exposing a shared learner to them. Robins <xref ref-type="bibr" rid="ref-4">[4]</xref> studied literal and pseudo-rehearsal. The present ER-IR comparison varies current-task inclusion and past-task sample dose along with task balancing, so it tests a combined replay distribution rather than pure temporal scheduling.</p>
<p>Regularization methods such as elastic weight consolidation protect important parameters <xref ref-type="bibr" rid="ref-6">[6]</xref>; surveys review replay and other continual-learning families <xref ref-type="bibr" rid="ref-7">[7]</xref>, <xref ref-type="bibr" rid="ref-8">[8]</xref>. A recent chapter places these methods within learning and memory <xref ref-type="bibr" rid="ref-9">[9]</xref>. Reinforcement-learning replay serves additional purposes such as decorrelation <xref ref-type="bibr" rid="ref-5">[5]</xref>. This experiment investigates exemplar memory and task sampling in a supervised task-incremental network, rather than testing all rehearsal mechanisms.</p>
<p>Continual learning involves sequential exposure to changing task data. Active learning concerns selective acquisition of labels and need not assume a stationary pool <xref ref-type="bibr" rid="ref-12">[12]</xref>; classical pooled prediction <xref ref-type="bibr" rid="ref-13">[13]</xref> addresses another evaluation setting. These distinctions motivate specifying the task sequence and memory-access assumptions directly.</p>
<p>This is a controlled empirical comparison rather than a new algorithm. Capacity, current-task inclusion, past-task sampling and layout statistics are explicit factors. Shared buffer contents improve the comparison, but differing replay doses and future-data access must be accounted for before attributing effects solely to scheduling. The present review does not establish that this particular factorization has never been studied.</p>
</sec>
<sec id="sec-4"><label>3</label><title>Research Gap and Questions</title>
<p>The study asks how nominal capacity relates to retention, how two replay distributions differ, and how per-task layout permutations affect interference. These are operational questions in a small network; the design does not equate layout correlation with semantic task similarity or measured hidden-representation overlap. RQ1 (capacity). How does replay buffer capacity, from 25 to 200 stored samples per class, govern Backward Transfer (BWT) and final average accuracy on sequential digit tasks?</p>
<p>H1. forgetting decreases monotonically — BWT rises toward zero — as capacity grows.</p>
<p>RQ2 (scheduling). Does distributing a fixed replay budget uniformly over past tasks (Interleaved Replay (IR)) reduce forgetting relative to drawing the same number of samples uniformly from the entire buffer (ER)?</p>
<p>H2. IR achieves higher BWT than same-budget ER, with the benefit expected at moderate buffer sizes (50– 100 samples per class).</p>
<p>RQ3 (similarity). Does the degree of input-layout overlap between tasks modulate interference?</p>
<p>H3. streams whose tasks share a common input layout (natural digit images) suffer more interference than streams whose tasks are decorrelated by per-task feature permutations, consistent with overlap-based accounts of catastrophic interference <xref ref-type="bibr" rid="ref-1">[1]</xref>, <xref ref-type="bibr" rid="ref-3">[3]</xref>.</p>
<p>The scope is one dataset subset and a compact task-incremental network. Task identity is supplied at test time. Shared controls and paired seeds support conditional comparisons, while the full pipeline and run-level statistics require artifact confirmation.</p>
</sec>
<sec id="sec-5"><label>4</label><title>Materials and Methods</title>
<sec id="sec-6"><label>4.1</label><title>Dataset and Task Construction</title>
<p>The scikit-learn load_digits benchmark contains 1797 images with 64 intensity features in an 8x8 grid and ten classes <xref ref-type="bibr" rid="ref-11">[11]</xref>. It is the test subset of the larger UCI Optical Recognition of Handwritten Digits collection, not the complete 5620-example corpus <xref ref-type="bibr" rid="ref-10">[10]</xref>. The source reports per-class counts 174-183, adjacent-pair tasks of 354-363 samples, 80/20 per-task splits with about 283-290 training and 71-73 test examples, and identical splits across conditions. Pixel intensities 0-16 are scaled to 0-1. Task identity is provided during training and evaluation.</p>
</sec>
<sec id="sec-7"><label>4.2</label><title>Network Architecture and Training Objective</title>
<p>The learner is a task-incremental multilayer perceptron implemented from scratch in NumPy <xref ref-type="bibr" rid="ref-13">[13]</xref>. The 64- dimensional input feeds a shared body of two fully connected hidden layers of 16 rectified-linear units each, followed by T = 5 task-specific readout heads, each a single logistic unit. Training a task activates only that task&apos;s head, so cross-task interference arises exclusively through the shared body — the component that continual learning must protect. The loss is binary cross-entropy computed on the active head; when a training batch contains replayed samples from earlier tasks, each replayed sample contributes through its own task&apos;s head while gradients still flow into the shared body. Optimization is mini-batch stochastic gradient descent with a fixed learning rate of 0.3 and a current-stream batch of 32 samples; no momentum, weight decay, or adaptive step sizes are used, keeping the dynamics transparent. Each task is trained for 1,000 gradient steps, which at the given batch size corresponds to roughly 111 epochs over the ~287 training samples of a task; this deliberately over-trained regime drives the current task to near-zero loss before the stream advances, maximizing the opportunity for interference. Weights are initialized with He scaling from a per-seed generator that is shared by all strategies, so that every strategy within a seed starts from an identical model.</p>
</sec>
<sec id="sec-8"><label>4.3</label><title>Replay Memory</title>
<p>The source reports class-balanced reservoirs built during the first exposure epoch, with nominal capacities 25, 50, 100, 200 per class and identical snapshot sequences for ER and IR. A reservoir stores at most min(c, n_class_train) distinct examples; with roughly 139-146 train examples per class, c=200 stores all available examples rather than 200 distinct examples. Actual retained counts and insertion timing must be reported. The standard replacement probability c/n_seen applies after capacity is reached; code must confirm that repeated epochs do not create duplicated reservoir entries.</p>
</sec>
<sec id="sec-9"><label>4.4</label><title>Continual-Learning Strategies</title>
<p>Naive training uses current-task batches. ER appends 32 samples drawn with replacement from the whole buffer, including the current task. IR draws 32 samples from past tasks only, balancing the chosen past task before selecting its stored example. With i completed tasks and equal full class buffers, ER&apos;s expected past-task dose is 32i/(i+1), while IR&apos;s is 32; actual partially filled buffers change the ER fraction. At the first task IR has no past examples, so its fallback and batch size require specification. Joint training draws from all tasks from the outset. Equal 5000 update counts do not imply equal sample computation or information access, and shared current-task batch sequences apply to the sequential conditions rather than joint pooled batches. Replay supplies extra sample evaluations per update relative to naive learning. The two replay distributions alter the proportion of old and current samples, and may differ at the empty-history boundary. Joint training has access to future-task data and changes the data distribution from the start. These are treatment differences that constrain an interpretation as pure temporal ordering or equal-compute superiority.</p>
</sec>
<sec id="sec-10"><label>4.5</label><title>Task-Similarity Manipulation</title>
<p>Natural tasks share the canonical coordinate frame, with reported mean pairwise correlation 0.82 between task-average images. Independent fixed task permutations lower this summary to -0.002. An invertible coordinate permutation preserves each task&apos;s information, but does not eliminate every cross-task relation or guarantee identical optimization difficulty. The observed layout manipulation is not a direct measurement of gradient conflict or hidden representation overlap.</p>
</sec>
<sec id="sec-11"><label>4.6</label><title>Evaluation Protocol</title>
<p>After task i, A[i, j] is held-out accuracy on task j. For zero-based task indices, BWT averages final accuracy minus just-trained accuracy over the four earlier tasks, while final average accuracy averages the final row. Negative BWT indicates forgetting; zero BWT may also mask compensating task gains and losses. Above-diagonal accuracy in non-joint runs calibrates untrained heads, whereas joint training has already accessed all tasks.</p>
<p>∑T j= − 02[AT−1,j − Aj,j] ∑T j= − 01 AT−1,j</p>
<p>BWT = , AAfinal = (1)</p>
<p>T − 1 T</p>
</sec>
<sec id="sec-12"><label>4.7</label><title>Fairness and Controls</title>
<p>Reported controls use common splits, initialization and reservoir snapshots with separate per-seed generators. Update counts are shared, but sample evaluations and future-task access differ. ER-versus-IR comparisons therefore concern the stated sampling distributions rather than identical doses of past experience. First-task fallback and actual batch normalization require code confirmation.</p>
</sec>
</sec>
<sec id="sec-13"><label>5</label><title>Experimental Setup</title>
<p>The source reports 200 runs: ten seeds 42-51, two layouts and ten conditions per layout, with 186 s runtime. NumPy <xref ref-type="bibr" rid="ref-13">[13]</xref>, scikit-learn <xref ref-type="bibr" rid="ref-11">[11]</xref>, SciPy <xref ref-type="bibr" rid="ref-14">[14]</xref> and Matplotlib <xref ref-type="bibr" rid="ref-15">[15]</xref> are named components. The original script and results file are needed to verify execution, controls, timings, and programmatic insertion of values.</p>
<p>Reported checks include first-task loss changing from 0.551 to 0.00016 and mean untrained accuracy 0.49 over non-joint runs. These are plausible sanity summaries, but pooled mean accuracy does not prove unbiased heads or absence of forward transfer on each task. Run-level loss and accuracy matrices are required to confirm the checks.</p>
<p>The source reports a balanced two-layout by ten-condition grid per seed, yielding 200 runs. Naive, ER/IR at four caps, and joint have ten replications per cell. Shared splits, initialization and snapshots support pairing among sequential strategies; original artifacts must confirm completeness, exclusions and exact controls. Paired two-sided Wilcoxon signed-rank tests summarize differences over ten seeds; rank-biserial effects and percentile bootstrap intervals are reported. The source uses 2000 bootstrap draws with generator 1234. The twenty-six comparisons are exploratory and unadjusted in the original tables. A Holm sensitivity calculation on the displayed rounded p-values gives minimum adjusted p = 0.052 across the complete 26-test family, so none reaches 0.05. Exact nonzero counts, ties and rounding require the raw results before assigning exact adjusted probabilities.</p>
<table-wrap id="tab1"><label>Table 1</label><caption><p>Architecture, training hyperparameters, and experimental design as reported in the source manuscript; execution requires the original code and run-level record.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:4:66,127.5,537.1,295.5"/><graphic content-type="pdf-region" xlink:href="pdf-region:5:66,613.1,537.1,755.6"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
</sec>
<sec id="sec-14"><label>6</label><title>Results</title>
<sec id="sec-15"><label>6.1</label><title>Catastrophic Forgetting under Naive Sequential Training</title>
<p><xref ref-type="fig" rid="fig1">Fig. 1</xref> shows strong retention losses in naive sequential training. The reported natural-layout BWT is -0.170 and final accuracy 0.861, versus joint-reference 0.996; task 0 falls from 1.000 just-trained to 0.754 after the fifth task. The other earlier-task losses are 13.1, 18.8 and 11.7 percentage points. Above-diagonal values are near binary chance in aggregate, which does not exclude individual-task forward transfer. Permuted naive BWT is -0.135 with larger seed variation.</p>
<fig id="fig1"><label>Figure 1</label><caption><p>Mean test-accuracy matrices A[i,j] over 10 seeds — accuracy on task j (columns) after completing task i (rows) — for naive sequential training, ER (100 samples per class), and IR (100 per class), in the natural (top row) and permuted (bottom row) regimes.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:5:66,214.5,537.1,480.9"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>Lower-triangular degradation relative to the diagonal shows forgetting; above-diagonal entries sit at the binary chance level (<xref ref-type="fig" rid="fig1">Fig. 1</xref>).</p>
</sec>
<sec id="sec-16"><label>6.2</label><title>Effect of Replay Buffer Size</title>
<p>Natural-layout ER means improve with nominal capacity: BWT -0.035, -0.024, -0.015, -0.010 at 25, 50, 100, 200, and accuracy 0.969, 0.977, 0.984, 0.988. Capacity 100 stores about 1000 examples, roughly 70 percent of the full training stream, and closes approximately 91 percent of the naive-to-joint accuracy gap. Capacity 200 saturates the available class training data. Monotonicity is an aggregate ER pattern, not a universal law: permuted IR BWT changes from -0.010 at 25 to -0.011 at 50 and its accuracy from 0.987 to 0.986. No formal paired monotone-trend test is provided.</p>
<table-wrap id="tab2"><label>Table 2</label><caption><p>Reported final accuracy and BWT over ten seeds.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:6:66,353.6,537.1,721.1"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<fig id="fig2"><label>Figure 2</label><caption><p>BWT versus replay buffer capacity (mean over 10 seeds; shaded bands are 95% bootstrap confidence intervals) for ER and IR in the natural (left) and permuted (right) regimes.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:6:66,171.6,537.1,362.6"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>Dashed and dotted horizontal lines mark the naive and joint-training BWT levels (<xref ref-type="fig" rid="fig2">Fig. 2</xref>).</p>
</sec>
<sec id="sec-17"><label>6.3</label><title>Learning Trajectories Across the Stream</title>
<p><xref ref-type="fig" rid="fig3">Fig. 3</xref> reports mean accuracy over encountered tasks. Naive means decline from 1.000 after the first task to 0.861 after the fifth, with intermediate 0.947, 0.939, 0.903. Replay curves remain near ceiling but do not all rise: seen-task averages can decrease as retained tasks accumulate. These averages mix current acquisition and past retention; a separate plasticity endpoint is required to claim replay does not impede new-task learning.</p>
<fig id="fig3"><label>Figure 3</label><caption><p>Natural-regime learning curves: Mean accuracy over all tasks encountered so far, measured after completing each training task, for naive training, replay at three capacities, IR (100 per class), and joint training (mean over 10 seeds).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:7:66,536.1,537.1,739.1"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
</sec>
<sec id="sec-18"><label>6.4</label><title>Interleaved Versus Uniform Replay</title>
<p><xref ref-type="table" rid="tab3">Table 3</xref> reports small ER-versus-IR differences. None of the eight BWT comparisons reaches an unadjusted 0.05 threshold; the smallest reported p is 0.058. One of sixteen interleaving comparisons, permuted accuracy at capacity 25, has unadjusted p = 0.036 but does not survive multiplicity control. These are failures to detect benefit under the reported protocol, not equivalence. Since ER includes current-task samples while IR excludes them, the contrast combines old-task dose with sampling allocation rather than isolating temporal scheduling.</p>
<table-wrap id="tab3"><label>Table 3</label><caption><p>Paired Wilcoxon signed-rank tests, IR versus ER, over 10 seeds per comparison.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:7:66,117.9,537.1,341.3"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<p>W = test statistic; p = two-sided p-value; r = matched-pairs rank-biserial correlation</p>
</sec>
<sec id="sec-19"><label>6.5</label><title>Effect of Task Similarity</title>
<p>Per-task permutation is associated with higher BWT in several unadjusted comparisons. Natural ER 25 has BWT -0.035 versus -0.014 under permutation (reported p = 0.002); naive and IR 200 do not reach unadjusted 0.05. Seven comparisons are below 0.05 before correction. Across the whole 26-test family none survives Holm on the rounded table values; treating the 10 layout tests as a separately prespecified family would retain only ER 25, IR 25 and ER 100. Such family choice is a sensitivity convention here, not retrospective preregistration. The observed effect is specific to coordinate layout and does not identify a general semantic-similarity or gradient mechanism.</p>
<fig id="fig4"><label>Figure 4</label><caption><p>Effect of task similarity on BWT (left) and final average accuracy (right) for selected conditions (mean over 10 seeds; error bars are 95% bootstrap confidence intervals).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:8:66,462.1,537.1,641.8"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>Natural = shared input layout; permuted = per-task feature permutations (natural and permuted variants of each condition shown side by side) (<xref ref-type="fig" rid="fig4">Fig. 4</xref>).</p>
<table-wrap id="tab4"><label>Table 4</label><caption><p>Task-similarity effect on BWT: Natural versus permuted regime (paired two-sided Wilcoxon signed-rank tests over 10 seeds).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:8:66,220.4,537.1,369.4"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
</sec>
</sec>
<sec id="sec-20"><label>7</label><title>Statistical Analysis</title>
<p><xref ref-type="table" rid="tab3">Tables 3</xref> and 4 retain the twenty-six reported unadjusted paired tests: sixteen ER-IR comparisons and ten layout comparisons. One ER-IR accuracy contrast has p = 0.036, whereas all BWT contrasts have p ≥ 0.058. Holm correction on displayed values across 26 tests gives no adjusted p below 0.05 (minimum 0.052); separating the exploratory questions into 16- and 10-test families leaves no ER-IR result and only three layout contrasts below 0.05. Exact adjustment should use unrounded raw probabilities and a justified family definition.</p>
<p>Unadjusted layout p-values range from 0.002 to 0.020 for the seven favourable contrasts, with large reported rank-biserial magnitudes. Separate marginal bootstrap intervals are not a paired confidence interval for the contrast, and nonoverlap cannot replace a multiplicity-aware paired test. Conversely, overlapping intervals do not prove equivalence or absence of a practically relevant difference.</p>
<p>Ten nonzero, untied pairs permit a minimum exact two-sided signed-rank probability 2/2^10=0.001953125. Zeros or ties alter the effective sample and implementation; the raw comparisons are needed to reconcile p-values and effect sizes. No formal equivalence margin, TOST, or preregistered testing family is supplied. Small observed mean differences do not establish a bound on the true paired effect.</p>
</sec>
<sec id="sec-21"><label>8</label><title>Discussion</title>
<p>Memory capacity is associated with most of the measured accuracy recovery in natural-layout ER, with diminishing increments at larger nominal caps. Capacity 100 retains roughly 70 percent of the training examples; capacity 200 retains the complete available class data. These are sizeable memory fractions on a small benchmark, so they should not be described as an arbitrary 200-sample reservoir comparison or a deployment provisioning rule.</p>
<p>The replay distributions differ in both current-task duplication and past-task dose. In a filled class-balanced buffer after i prior tasks, ER includes the current task with probability 1/(i+1); that fraction declines as tasks accumulate. IR samples only prior tasks and assigns them equal task probabilities. At early stages and while reservoirs fill, these differences can be larger. The absence of a detectable ER-IR benefit does not test pure temporal ordering at an identical past-task dose or prove a specific complementary-learning-systems mechanism. A matched past-dose and current-duplication experiment would isolate those questions.</p>
<p>Input permutations reduce task-average image correlation and are associated with improved retention in some exploratory conditions. Gradient conflict is a possible explanation, but no gradient dot products, hidden-state overlaps, or parameter-importance trajectories are measured. Invertible task encodings preserve information but can change optimization and representation geometry; both may contribute. The evidence therefore supports a benchmark-specific association with layout rather than a verified interference mechanism.</p>
<p>The layout contrast is not significant for naive training, where seed variation is largest. A formal interaction test is needed before concluding that layout effects require partial memory. Separate significances across capacities do not by themselves establish an interaction.</p>
<p>This task-incremental study compares exemplar replay distributions only. It does not benchmark against regularization methods <xref ref-type="bibr" rid="ref-6">[6]</xref>, generative replay, larger networks or class-incremental inference. Shared reported controls improve comparison, but implementation verification and the replay-dose distinction remain necessary for strong internal-validity claims.</p>
</sec>
<sec id="sec-22"><label>9</label><title>Implications</title>
<p>The experiment separates nominal memory capacity, a combined replay-distribution treatment, and coordinate layout. It finds a marked capacity association and no multiplicity-supported ER-IR advantage. Joint future-data access has a different information regime, so its performance cannot identify biological interleaving as the causal mechanism. Layout is a useful controlled benchmark variable, while generalization and mechanistic attribution require new measurements.</p>
<p>On this small digit subset, a nominal cap 100 retains about 70 percent of the training stream and yields ER accuracy 0.984 versus joint-reference 0.996; retaining all class training examples at cap 200 yields 0.988. These fractions and task-identity assumptions limit direct use for edge-memory provisioning. Per-task remapping may affect optimization and interface requirements, so a deployment recommendation cannot follow from this single layout experiment.</p>
</sec>
<sec id="sec-23"><label>10</label><title>Limitations</title>
<p>The learner is small and uses task-specific heads with known task identity. The dataset is the 1797-example scikit-learn subset rather than the full UCI corpus. Coordinate permutations modify layout statistics without changing task labels, but do not remove every shared representation. Nominal cap 200 exceeds available class training examples. The ER-IR contrast changes past-task sample dose and current-task inclusion, and the first-task fallback needs confirmation.</p>
<p>The seed-level tests are exploratory, contain ties and zeros, and lack a preregistered family, paired-effect interval, power analysis or equivalence margin. Multiplicity sensitivity limits strong significance claims. The source reports 1000 updates per task, around 111 current-sample exposures in epoch units, with additional replay work; normalization and actual counts require code. Hyperparameters are shared, and no regularization or generative baselines are evaluated.</p>
</sec>
<sec id="sec-24"><label>11</label><title>Reproducibility and Data Availability</title>
<p>The source specifies code/experiment.py, results/results.json and code/build_docx.py as the generation workflow, with seeds 42-51 and bootstrap seed 1234. It reports Python 3.12.14, NumPy 2.1.3, SciPy 1.14.1, scikit-learn 1.5.2 and Matplotlib 3.9.2, and 186 s runtime. A stable repository or supplied code/result package is required to verify these commands, environments and generation claims; it is not identified here.</p>
<p>The public dataset is the bundled scikit-learn digit subset <xref ref-type="bibr" rid="ref-10">[10]</xref>, <xref ref-type="bibr" rid="ref-11">[11]</xref>. Exact splits, permutations, reservoirs, fallback batches and the claimed 200 accuracy matrices require the seeded original script and machine-readable output. No new human or animal participants are described, but applicable dataset-use conditions and declarations should be confirmed by the authors.</p>
</sec>
<sec id="sec-25"><label>12</label><title>Conclusion</title>
<p>The reported task-incremental experiment compares exemplar replay on a small real digit subset. Natural-layout ER has improving mean retention with larger nominal capacity, from BWT -0.170 without replay to - 0.010 when the available training examples are retained, while final accuracy rises from 0.861 to 0.988 against a joint future-data reference of 0.996. These means are descriptive of the stated protocol rather than a monotonicity theorem.</p>
<p>No reported BWT contrast detects an ER-IR advantage at unadjusted 0.05, and none survives whole-family multiplicity control. Because past-task dose and current-task inclusion differ, this is a comparison of sampling distributions rather than a pure interleaving test. Layout permutations have favourable exploratory contrasts, but no whole-family Holm significance or measured gradient mechanism is established. The results provide a bounded empirical baseline whose reproduction and broader relevance require code, raw matrices and larger, matched-treatment studies.</p>
</sec>
<sec id="sec-26"><title>Data Availability</title>
<p>The study uses the 1797-example load_digits subset of UCI Optical Recognition of Handwritten Digits. Splits and permutations are experimental transformations; their exact reconstruction requires the original seed-controlled code.</p>
<p>The study describes no newly recruited human or animal participants. Dataset permissions and any applicable review requirements must be confirmed by the authors.</p>
</sec>
<sec id="sec-27"><title>Funding</title>
<p>This work received no external funding.</p>
</sec>
<sec id="sec-28"><title>Competing Interests</title>
<p>none declared.</p>
</sec>
<sec id="sec-29"><title>Author Contributions</title>
<p>Actual contributions and author approval must be confirmed by the named researchers. A program or AI tool cannot establish the author roles or scientific accountability.</p>
</sec>
<sec id="sec-30"><title>AI Disclosure</title>
<p>The source reports AI-assisted drafting. The authors must confirm the tools and tasks actually used, the origin of experiment code and data, and responsibility for independent scientific verification.</p>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="ref-1"><label>1</label><mixed-citation>McCloskey, M., &amp; Cohen, N. J. (1989). Catastrophic interference in connectionist networks: The sequential learning problem. In Psychology of learning and motivation (Vol. 24, pp. 109–165). Elsevier. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/S0079-7421(08)60536-8">https://doi.org/10.1016/S0079-7421(08)60536-8</ext-link></mixed-citation></ref>
<ref id="ref-2"><label>2</label><mixed-citation>Rumelhart, D. E., Hinton, G. E., &amp; Williams, R. J. (1986). Learning representations by back-propagating errors. Nature, 323(6088), 533–536. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/323533a0%0A%0A">https://doi.org/10.1038/323533a0%0A%0A</ext-link></mixed-citation></ref>
<ref id="ref-3"><label>3</label><mixed-citation>McClelland, J. L., McNaughton, B. L., &amp; O’Reilly, R. C. (1995). Why there are complementary learning systems in the hippocampus and neocortex: Insights from the successes and failures of connectionist models of learning and memory. Psychological review, 102(3), 419. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1037/0033-295X.102.3.419">https://doi.org/10.1037/0033-295X.102.3.419</ext-link></mixed-citation></ref>
<ref id="ref-4"><label>4</label><mixed-citation>Robins, A. (1995). Catastrophic forgetting, rehearsal and pseudorehearsal. Connection science, 7(2), 123–146. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1080/09540099550039318">https://doi.org/10.1080/09540099550039318</ext-link></mixed-citation></ref>
<ref id="ref-5"><label>5</label><mixed-citation>Sutton, R. S., Barto, A. G., &amp; Barto, A. (1998). Reinforcement learning: An introduction (Vol. 1). MIT press, Cambridge. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1017/S0263574799271172">https://doi.org/10.1017/S0263574799271172</ext-link></mixed-citation></ref>
<ref id="ref-6"><label>6</label><mixed-citation>Kirkpatrick, J., Pascanu, R., Rabinowitz, N., Veness, J., Desjardins, G., &amp; Rusu, A. A. (2017). Overcoming catastrophic forgetting in neural networks. Proceedings of the national academy of sciences, 114(13), 3521–3526. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1073/pnas.1611835114">https://doi.org/10.1073/pnas.1611835114</ext-link></mixed-citation></ref>
<ref id="ref-7"><label>7</label><mixed-citation>De Lange, M., Aljundi, R., Masana, M., Parisot, S., Jia, X., Leonardis, A., … &amp; Tuytelaars, T. (2021). A continual learning survey: Defying forgetting in classification tasks. IEEE transactions on pattern analysis and machine intelligence, 44(7), 3366–3385. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/TPAMI.2021.3057446">https://doi.org/10.1109/TPAMI.2021.3057446</ext-link></mixed-citation></ref>
<ref id="ref-8"><label>8</label><mixed-citation>Wang, L., Zhang, X., Su, H., &amp; Zhu, J. (2024). A comprehensive survey of continual learning: Theory, method and application. IEEE transactions on pattern analysis and machine intelligence, 46(8), 5362–5383. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/TPAMI.2024.3367329">https://doi.org/10.1109/TPAMI.2024.3367329</ext-link></mixed-citation></ref>
<ref id="ref-9"><label>9</label><mixed-citation>de Ven, G. M., Soures, N., &amp; Kudithipudi, D. (2024). Continual learning and catastrophic forgetting. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/B978-0-443-15754-7.00073-0">https://doi.org/10.1016/B978-0-443-15754-7.00073-0</ext-link></mixed-citation></ref>
<ref id="ref-10"><label>10</label><mixed-citation>Alpaydin, E., &amp; Kaynak, C. (1998). Optical recognition of handwritten digits. UCI Machine Learning Repository. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.24432/C50P49">https://doi.org/10.24432/C50P49</ext-link></mixed-citation></ref>
<ref id="ref-11"><label>11</label><mixed-citation>Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., &amp; Grisel, O. (2011). Scikit-learn: Machine learning in Python. The journal of machine learning research, 12, 2825–2830. <ext-link ext-link-type="uri" xlink:href="https://jmlr.org/papers/v12/pedregosa11a.html">https://jmlr.org/papers/v12/pedregosa11a.html</ext-link></mixed-citation></ref>
<ref id="ref-12"><label>12</label><mixed-citation>Settles, B. (2009). Active learning literature survey. <ext-link ext-link-type="uri" xlink:href="https://burrsettles.com/pub/settles.activelearning.pdf">https://burrsettles.com/pub/settles.activelearning.pdf</ext-link></mixed-citation></ref>
<ref id="ref-13"><label>13</label><mixed-citation>Harris, C. R., Millman, K. J., Van Der Walt, S. J., Gommers, R., Virtanen, P., &amp; Cournapeau, D. (2020). Array programming with NumPy. Nature, 585(7825), 357–362. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/s41586-020-2649-2">https://doi.org/10.1038/s41586-020-2649-2</ext-link></mixed-citation></ref>
<ref id="ref-14"><label>14</label><mixed-citation>Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., &amp; Cournapeau, D. (2020). SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nature methods, 17(3), 261–272. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/s41592-019-0686-2">https://doi.org/10.1038/s41592-019-0686-2</ext-link></mixed-citation></ref>
<ref id="ref-15"><label>15</label><mixed-citation>Hunter, J. D. (2007). Matplotlib: A 2D graphics environment. Computing in science &amp; engineering, 9(3), 90–95. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/MCSE.2007.55">https://doi.org/10.1109/MCSE.2007.55</ext-link></mixed-citation></ref>
</ref-list></back>
</article>