Design of the duration discrimination task in Experiment 1. Each trial consisted of an isochronous standard sequence of five sounds (four intervals), followed by silence and another pair of sounds. The comparison duration was either shorter or longer than the standard intervals and took on one of ten values (DEV) that were proportional to the inter-onset interval (IOI) between tones making up the standard sequence. The task was to press the S or L key to indicate whether the comparison interval was shorter or longer than the standard IOI. Over the course of 400 unique trials of a single session, IOI ranged from 2 to 998 ms. In linear-order sessions, IOI increased in each trial in the first 200 trials and decreased in the other half of the trials (or vice versa, counter-balanced across participants) in steps of 4 ms. In random-order sessions, change in stimulus rate between a given trial n and immediately preceding trial n-1 (ΔIOI) was maximized, and the distribution of ΔIOI ranged from -778 ms to +770 ms.

Descriptive statistics and test results for comparison of Beta estimates against null distributions in Experiment 1 analyses.

Main findings of Experiment 1. A Left: Each circle represents a single participant’s estimate from random-order session (x axis) and linear-order session (y axis). The histograms across the axes show the distributions of estimates for each session type. The dotted lines represent 1:2 and 2:1 ratio between the axes, and the straight line represents one-to-one correspondence. Right: permutation test results. The distribution of summed residuals (distance of data points to the closest y=x, y=2*x and y=x/2 lines) of shuffled data over 1000 iterations, and the summed residual from original data (dashed line) which fell below .008 of the permutation distribution. B Average accuracy from random-order (left, green) and linear-order (right, blue) sessions. C Effects of |±ΔIOI| on responses in Experiment 1 duration discrimination paradigm. Left: Flexibility estimates. Each circle represents an individuals’ slope (β) obtained from logistic models, fitted separately to conditions where |-ΔIOI| (left, green) or |+ΔIOI| (right blue) predicted accuracy, with greater values (arrow’s direction) indicating better oscillator flexibility. The distribution β from both conditions were smaller than zero, indicating a negative effect of between-trial absolute rate change on accuracy. Right: Participants’ average bias from |-ΔIOI| (left, green), and |+ΔIOI| (right, blue) conditions. Box plots in B and C panels show median (black vertical line), 25th and 75th percentiles (box edges) and extreme datapoints (whiskers).

Descriptive statistics of unpaced tapping measures in first and second experiments, and test results for pairwise comparisons.

Experiment 2 (A) timeline, and illustrations of the duration discrimination (B), paced tapping (C), slider (D) and keypress (E) tasks.

Results of Experiment 2 preferred rate analyses. A Top: Estimates of preferred rate from each task condition. Box plots show median (black vertical line), 25th and 75th percentiles (box edges) and remaining data range (whiskers). The horizontal dashed lines represent the minimum and maximum stimulus rates presented in the experiment. Bottom: Pairwise correlations between preferred rates across tasks. For the slider and key-press tasks, boxes are colored to indicate fast-start (blue) and slow-start (pink) conditions. Correlations and p-values are reported for significant correlations only. B Relationship between the preferred rate estimates from the paced tapping, duration discrimination, and SMT (first measurement) tasks. Left: participants’ estimates from the three tasks. Each circle represents an individual’s preferred rate estimate, connected by lines between the tasks. Both circles and lines are color-sorted by individuals’ SMT, ranging from magenta (fast) to blue (slow). Right: correlation between the difference scores. Each circle represents a single participant’s difference score, namely, how different the estimates from SMT (x axis) and paced tapping (y axis) tasks were than those from the duration discrimination task. Straight black line represents the regression line, dashed lines represent 95% confidence intervals.

Results of Experiment 2 flexibility analyses. A-D Effects of between-trial absolute rate-change on performance in Experiment 2 duration discrimination (A, C) and paced tapping (B,D) tasks. In A-B, each circle represents an individuals’ slope (β) obtained from models, fitted separately to conditions where |-ΔIOI| (left, green) or |+ΔIOI| (right, blue) predicted accuracy in duration discrimination (A) or |TME| in paced tapping (B) task. In C-D, box plots show average bias in duration discrimination (C) and average TME in paced tapping (B) tasks in |-ΔIOI| (left, green) and|+ΔIOI| (right, blue) conditions. In all panels, box plots show median (black vertical line), 25th and 75th percentiles (box edges) and extreme datapoints (whiskers). E-F Correlations between individuals’ age and the flexibility estimates from duration discrimination (E) and paced tapping (F) tasks. Straight black lines represent the regression line, dashed lines represent 95% confidence intervals. Histograms above each plot show the distribution of participants’ age after outlier corrections.

Descriptive statistics and test results for comparison of Beta estimates against null distributions in Experiment 2 analyses.

Illustration of the optimization procedure and parameter choices for smoothing accuracy in Experiment 1. A Bottom: An example participant’s linear-order session dataset. Each color represents an output of the smoothing function that uses a window size, ranging from 10 (yellow) to 50 (dark blue). Top: The number of maximum values on the smoothed accuracy for each window size. B Participants’ average number of curve maxima for random-order (blue) and linear-order (magenta) sessions. Arrows show the optimized window sizes for the session types, where each individual’s dataset had only one curve maximum (dashed line).

Results of the bootstrapping analysis. In A, each circle shows harmonic difference between preferred rate estimates from the original and downsampled datasets for session 1 (blue) and session 2 (magenta) and their average (dotted black line) at the respective step size. In B, each circle shows harmonic difference between preferred rate estimates from the downsampled session 1 and session 2 datasets at the respective step size.