(I) Adaptation
Figure 3 shows that the speed profiles of the drawer opening movement were perturbed by the spontaneous increase in drawer resistance, but adapted in subsequent trials.
Initial Error between Baseline and First Adaptation Trial
The 3 (age) x 2 (initial error: baseline vs. first adaptation trial) repeated measures ANOVAs for each dependent variable (peak speed, movement time, number of movement units) revealed significant initial errors of all variables across all age groups.
In detail, peak speed showed a main effect of initial error, F(1, 54) = 65.99, p < .001, η2p = 0.55, with a significant decrease from baseline (M = 354.57 ± 10.19 mm/s) to the first adaptation trial (M = 266.12 ± 10.28 mm/s). No main effect of age, F(2, 54) = 1.17, p = .317, η2p = 0.04, or interaction between age and initial error, F(2, 54) = 1.46, p = .242, η2p = 0.05, was found (Fig. 4a).
Movement time revealed a significant main effect of initial error, F(1, 52) = 30.14, p < .001, η2p = 0.37, with an increase from baseline (M = 990.89 ± 35.38 ms) to the first adaptation trial (M = 1338.91 ± 77.10 ms). A main effect of age was detected, F(2, 52) = 6.08, p = .004, η2p = 0.19, with shorter movement time in adults (M = 945.28 ± 37.82 ms) compared to 3-year-olds (M = 1242.21 ± 85.66 ms), pbonf = .035, d = ˗0.71, and compared to 1.5-year-olds (M = 1319.41 ± 129.74 ms), pbonf = .005, d = ˗0.89. There was no significant interaction between age and initial error, F(2, 52) = 1.06, p = .355, η2p = 0.04 (Fig. 4b).
Movement units displayed a significant main effect of initial error, F(1, 54) = 18.02, p < .001, η2p = 0.25 with an increase from baseline (M = 1.33 ± 0.04) to the first adaptation trial (M = 1.65 ± 0.08). A main effect of age was identified, F(2, 54) = 18.49, p < .001, η2p = 0.41. Post-hoc tests revealed fewer movement units in adults (M = 1.13 ± 0.10) compared to 3-year-olds (M = 1.65 ± 0.08), pbonf < .001, d = ˗1.23, and compared to 1.5-year-olds (M = 1.69 ± 0.11), pbonf < .001, d = ˗1.33. No significant interaction between age and initial error was determined, F(2, 55) = 0.23, p = .799, η2p = 0.01 (Fig. 4c). A histogram for the number of movement units for all age groups is attached in the Supplementary Figure S2 online.
Residual Errors in Adaptation Block
The 3 (age) x 2 (residual error: baseline vs. mean last three adaptation trial) repeated measures ANOVAs revealed non-significant residual errors of all variables across all age groups.
Hence, peak speed displayed no main effect of residual error, F(1, 55) = 3.70, p = .059, η2p = 0.06), indicating no significant differences between end of adaptation (M = 341.34 ± 11.36 mm/s) and baseline (M = 360.20 ± 11.49 mm/s). A main effect of age was identified, F(2, 55) = 7.10, p = .002, η2p = 0.21), with higher peak speed in adults (M = 360.47 ± 11.36 mm/s) compared to 3-year-olds (M = 303.10 ± 15.15 mm/s), pbonf = .044, d = 0.71, and in 3-year-olds compared to 1.5-year-olds (M = 388.22 ± 25.75 mm/s), pbonf = .002, d = ˗1.06. The interaction between age and residual error was not significant, F(2, 55) = 1.43, p = .249, η2p = 0.05 (Fig. 4a).
Movement time obtained no main effect of residual error, F(1, 55) = 1.21, p = .276, η2p = 0.02, suggesting that movement time did not differ significantly between baseline (M = 989.64 ± 34.02 ms) and end of adaptation (M = 1031.58 ± 44.42 ms). A main effect of age was detected, F(2, 55) = 10.91, p < .001, η2p = 0.28, with shorter movement time in adults (M = 831.72 ± 25.09 ms) compared to 3-year-olds (M = 1152.61 ± 60.24 ms), pbonf < .001, d = ˗1.19, and compared to 1.5-year-olds (M = 1052.34 ± 84.54 ms), pbonf = .008, d = ˗0.82). No significant interaction between age and residual error was observed, F(2, 55) = 2.32, p = .108, η2p = 0.08 (Fig. 4b).
Movement units revealed no main effect of residual error, F(1, 55) = 0.09, p = .763, η2p < 0.01. Thus, the number of movement units did not differ significantly between baseline (M = 1.33 ± 0.04) and end of adaptation (M = 1.32 ± 0.05). We found a main effect of age, F(2, 55) = 34.55, p < .001, η2p = 0.56, with adults (M = 1.00 ± 0.00), having fewer movement units compared to 3-year-olds (M = 1.53 ± 0.07), pbonf < .001, d = ˗1.99, and compared to 1.5-year-olds (M = 1.46 ± 0.08), pbonf < .001, d = ˗1.72. We found a significant interaction between age and residual error, F(2, 55) = 4.83, p = .012, η2p = 0.15. Compared to the baseline, adults had equal movement units at the end of adaptation (Mdiff = 0.00 ± 0.00), 3-year-olds had more movement units at the end of adaptation (Mdiff = 0.14 ± 0.01), and 1.5-year-olds had fewer movement units at the end of adaptation (Mdiff = ˗0.17 ± 0.04) (Fig. 4c).
Linear Mixed Model Trial-wise Adaptation Analysis
Linear mixed models were applied to each dependent variable, with trials (13–24) and age as fixed effects, and subject as a random factor, to analyze the trial-to-trial adaptation process between age groups (cf., statistical analyses). Model fits and estimates are presented in the Supplementary Tables S4-6 online.
A significant main effect of trials on peak speed was found, F(1, 618.29) = 29.72, p < .001, ηp2 = 0.05. This supports the adaptation process (i.e. an increase) of peak speed with subsequent trials. Peak speed was not significantly different between age groups across all trials, F(2, 504.80) = 0.20, p = .822, ηp2 < 0.01. There was no significant interaction effect between age and trials, F(2, 618.29) = 1.77, p = .171, ηp2 < 0.01 (Fig. 5a).
We found a significant main effect of trials on movement time, F(1, 602.46) = 33.72, p < .001, ηp2 = 0.05. Thus, movement time adapted (i.e. increased) with subsequent trials. Further, a main effect of age was identified, F(2, 570.34) = 7.66, p < .001, ηp2 = 0.03. Customized contrasts showed that adults had shorter movement time compared to 3-year-olds (p < .001), and compared to 1.5-year-olds (p = .008). We found a significant interaction between age and trials, F(2, 602.46) = 3.30, p = .038, ηp2 = 0.03. Post-hoc analysis with customized contrasts and trend function revealed that adults exhibited smaller reduction in movement time across trials with an estimated slope of ˗10.1 ms per trial. In contrast, 3-year-olds showed a significant reduction in movement time of ˗21.2 ms per trial, while 1.5-year-olds had the greatest reduction in movement time with ˗33.6 ms per trial. Statistically adults reduced their movement time significantly less than 1.5-year-olds (p = .028). Descriptively, adults reduced their movement time to baseline level already in the second trial, while children needed more trials to adapt their movement time and had a higher trial-to-trial variability (see Fig. 5b). To check the trial-by-trial variability statistically we conducted a 1 x 3(age) ANOVA on the participants individual standard deviation of the movement time divided by the individual mean of the movement time10. A significant main effect of age on trial-to-trial variability was found, F(2, 55) = 15.48, p < .001, ηp2 = 0.36. Post-hoc tests revealed a significantly lower trial-to-trial variability in adults (M = 0.13 ± 0.02 ms) compared to 3-year-olds (M = 0.31 ± 0.03 ms), pbonf < .001, d = ˗1.23, and compared to 1.5-year-olds (M = 0.39 ± 0.05 ms), pbonf < .001, d = ˗1.73.
For movement units, we found a main effect of trials, F(1, 610.83) = 20.85, p < .001, ηp2 = 0.03. Movement units were reduced with subsequent trials. A main effect of age, F(2, 663.56) = 3.79, p = .023, ηp2 = 0.01 was found, with adults having fewer movement units compared to 3-year-olds (p < .001) and compared to 1.5-year-olds (p < .001). No interaction between age and trials was detected, F(2, 610.82) = 0.21, p = .813, ηp2 < 0.01 (Fig. 5c).
(II) Deadaptation
Figure 6 presents the averaged speed profiles of the drawer movement within the 12 trails of the deadaptation block compared to the averaged baseline trials and the averaged last three adaptation trials for all three age groups.
Aftereffects Between First Deadaptation Trial and Baseline
To investigate whether aftereffects of the adaptation block existed when reducing the drawer resistance back to baseline level, we conducted for each dependent variable a 3 (age) x 2 (aftereffect: baseline vs. first deadaptation trial) repeated measures ANOVA.
Peak speed showed a significant main effect of aftereffects, F(1, 53) = 10.09, p = .002, η2p = 0.16. Thus, peak speed in the first trial of deadaptation (M = 424.97 ± 19.84 mm/s) overshot peak speed of the baseline (M = 358.57 ± 11.84 mm/s) significantly. A main effect of age was observed, F(2, 53) = 5.12, p = .009, η2p = 0.16, with 3-year-olds (M = 346.96 ± 18.56 mm/s) having lower peak speed than 1.5-year-olds (M = 439.40 ± 38.96 mm/s), pbonf = .007, d = ˗0.79. No significant interaction between age and aftereffect was revealed, F(2, 53) = 0.63, p = .537, η2p = 0.02 (Fig. 4a).
Movement time exhibited a main effect of aftereffects, F(1, 52) = 5.32, p = .025, η2p = 0.09. Movement time at the first deadaptation trial (M = 844.73 ± 44.60 ms) undershot the baseline average (M = 971.80 ± 33.43 ms). A main effect of age was shown, F(2, 52) = 7.56, p = .001, η2p = 0.23, indicating that adults had shorter movement time (M = 774.30 ± 30.26 ms) compared to 3-year-olds (M = 957.58 ± 44.45 ms), pbonf = .014, d = ˗0.66, and compared to 1.5-year-olds (M = 1000.35 ± 97.98 ms), pbonf = .002, d = ˗0.81. No significant interaction between age and aftereffect was detected, F(2, 52) = 0.31, p = .734, η2p = 0.01 (Fig. 4b).
Movement units did not reveal a main effect of aftereffect, F(1, 54) = 1.14, p = .290, η2p = 0.02. However, a main effect of age was present, F(2, 54) = 11.49, p < .001, η2p = 0.30, with fewer movement units in adults (M = 1.10 ± 0.09) compared to 3-year-olds (M = 1.47 ± 0.11), pbonf = .001, d = ˗0.88, and compared to 1.5-year-olds (M = 1.54 ± 0.10), pbonf < .001, d = ˗1.06. There was no significant interaction between age and aftereffect, F(2, 54) = 0.65, p = .527, η2p = 0.02 (Fig. 4c).
Residual Error in Deadaptation Block
To investigate whether participants deadapted their movements to the baseline level at the end of deadaptation, we used again 3 (age) x 2 (residual error: mean of last three deadaptation trials vs. baseline) repeated measures ANOVAs for each dependent variable.
Peak speed showed a main effect of residual error, F(1, 55) = 10.45, p = .002, η2p = 0.16, with still higher peak speed at the end of deadaptation (M = 392.71 ± 11.88 mm/s) compared to baseline (M = 360.20 ± 11.49 mm/s). A significant main effect of age was shown, F(2, 55) = 4.67, p = .013, η2p = 0.15, with lower peak speed in 3-year-olds (M = 335.54 ± 14.91 mm/s) than in 1.5-year-olds (M = 409.56 ± 27.29), pbonf = .012, d = ˗0.87. No significant interaction was revealed between age and residual error, F(2, 55) = 0.25, p = .780, η2p < 0.01 (Fig. 4a).
For movement time, we found a significant main effect of residual error, F(1, 55) = 17.61, p < .001, η2p = 0.24, with still higher movement time at the end of deadaptation (M = 986.64 ± 34.02 ms) compared to baseline (M = 855.32 ± 28.83 ms). A main effect of age, F(2, 55) = 9.69, p < .001, η2p = 0.26, showed that adults (M = 777.18 ± 25.82 ms) had shorter movement time than 3-year-olds (M = 1011.86 ± 48.54 ms), pbonf < .001, d = ˗1.08, and shorter movement time than 1.5-year-olds (M = 981.47 ± 67.08 ms), pbonf = .003, d = ˗0.94. The interaction between age and residual error was not significant, F(2, 55) = 1.33, p = .272, η2p = 0.05 (Fig. 4b).
Movement units displayed a significant main effect of residual error, F(1, 55) = 24.20, p < .001, η2p = 0.31, with fewer movement units at the end of deadaptation (M = 1.17 ± 0.03) compared to baseline (M = 1.33 ± 0.04). A main effect of age was determined, F(2, 55) = 40.34, p < .001, η2p = 0.60, with fewer movement units in adults (M = 1.00 ± 0.00) compared to 3-year-olds (M = 1.33 ± 0.05), pbonf < .001, d = ˗1.62, and compared to 1.5-year-olds (M = 1.43 ± 0.06), pbonf < .001, d = ˗2.16. The interaction between age and residual error was significant, F(2, 55) = 6.34, p = .003, η2p = 0.19. Adults and 1.5-year-olds had fewer movement units at the end of adaptation compared to baseline, while 3-year-olds had slightly more movement units (Mdiff_adults = 0.00 ± 0.00, Mdiff_3−yo = 0.27 ± 1.33, Mdiff_1.5−yo = ˗0.22 ± 1.434) (Fig. 4c).
Linear Mixed Model Trial-wise Deadaptation Analysis
Linear mixed models were applied to each dependent variable, with trials (25–36) and age as fixed effects, and subject as a random factor.
For peak speed, we found a significant main effect of trials, F(1, 613.52) = 12.27, p < .001, ηp2 = 0.02, indicating that peak speed was significantly deadapted. No significant main effect of age, F(2, 653.55) = 3.01, p = .050, ηp2 < 0.01 was detected. The interaction between age and trials, F(2, 613.52) = 1.22, p = .295, ηp2 < 0.01 was not significant. This suggests that peak speed was not reduced differently with subsequent trials between age groups (Fig. 5a).
For movement time, we found no main effect of trials, F(1, 599.77) = 0.23, p = .634, ηp2 < 0.01. This suggests that movement time was not deadapted. We also found no main effect of age, F(2, 644.76) = 0.90, p = .408, ηp2 < 0.01. Furthermore, the interaction between trials and age was not significant, F(2, 599.74) = 0.75, p = .475, ηp2 < 0.01 (Fig. 5b).
However, a main effect of trials on the number of movement units, F(1, 616.97) = 9.49, p = .002, ηp2 = 0.02 was found. Thus, the number of movement units decreased with subsequent trials. No main effect of age, F(1, 651.05) = 2.06, p = .128, ηp2 < 0.01 was revealed as well as and no interaction between age and trials, F(2, 616.96) = 0.49, p = .610, ηp2 < 0.01 (Fig. 5c).