Contrast manipulation used in the experiment

Top shows two example stimuli illustrating minimum (left) and maximum (right) contrast values. Bottom panel shows the prediction for the Piéron, the Fechner and the linear laws for all contrast levels (C) used in the study for a fixed set of parameters. The y-axis refers the time predicted by each law given a contrast value (x-axis) and the chosen set of parameters. α, β and ν are respectively the estimated participant specific intercept, slope and exponent for the three laws. The Fechner diffusion model additionally includes non-decision and decision threshold parameters (see methods).

Square root of the mean prediction error (milliseconds) from the leave-one-out procedure for the models applied to the intervals between each event (Ev.) including stimulus onset (S.) and response (R.).

Bold indicates the best model per interval. The R2 in parentheses refers to the fit of the predicted vs. observed mean durations. Prop. corr. refers to the R2 for the proportion of correct responses observed vs proportion predicted by the fit of the Fechner diffusion model to each intervals.

Behavioral results partially support Fechner’s law

Left: Mean RT (dot) for each contrast level and averaged predictions of the individual fits (line) with Accuracy (Top) and Speed (Bottom) instructions. Right: Mean proportion of correct responses averaged over trials and participant for each contrast level (triangles) along the average predicted proportions for the Fechner diffusion models (line) in Accuracy (Top) and Speed (Bottom).

HMP estimation reveals five sequential cognitive events

Top: Trial and participant averaged cumulative event occurrence probabilities (grey lines) for the five detected hidden multivariate events. Overlaid is the average representation of the time location (vertical lines) and electrode contribution of the five identified events based on their trial-by-trial maximum probability. Bottom: t-values from one-sample t-tests on the by-participant differences for each symmetric electrode pair for each event. White dots indicate electrode pairs that deviate significantly from 0 after a Bonferroni-correction with a base alpha level of 0.05 (corrected p-value = 0.00071), black dots denote non-significant electrode pairs.

Summary of the linear mixed models coefficients with the mean and the 95% credible intervals from the posterior distributions for the linear mixed models on intervals for each event (Ev.) including stimulus onset (S.) and response (R.).

Inter-event interval as a function of experimental factors

A) The panels represent the average duration between events for each contrast level, averaged across participants and trials (stimulus and response respectively as first and last events) for accuracy (top) and speed instructions (bottom). The lines represent the fits for the linear (black), Piéron (green) and Fechner (yellow) models. Winning model lines are represented as thicker than the other for each panel. B) Same as in A but with similar intervals summed together to form encoding (Ev.1-Ev2 + Ev.2-Ev.3) and decision (Ev.3-Ev4 + Ev.4-Ev.5) times. C) Observed proportion of correct responses (triangles) compared to the prediction (lines) from the Fechner diffusion model fit to the decision interval in B.

Electrophysiological analysis of the time leading to each event

A) Stimulus centered ERPs obtained by averaging three electrodes in occipital, parietal and frontal areas across trials and participants for 10 binned contrast levels in both speed-accuracy instructions. Vertical lines indicate the average peak time of each HMP event (later = darker). To statistically test for an effect of contrast we performed a non-parametric cluster-level paired t-test for high (> 50%) vs. low (< 50%) contrasts on the participant averaged waveforms. Significant clusters are marked as black dots at the corresponding times in each panel. B) Illustration of the event match computation. A single-trial match time-series (right most panel) is obtained by taking the dot product between the single-trial electrode time-series and the weights of an events’ average topography. C) HMP detected events time-courses obtained by using the average topographies (top panel) as a spatial filter for the stimulus onset to 100ms after the single-trial event peak time. For each speed-accuracy instruction the data is presented at the single-trial level (surface plot), with trials sorted by the single-trial most likely peak time of the event, and averaged across trials for 10 contrast bins after aligning to the most likely peak time of the event (waveform plot). On the surface plot a z-scoring was performed after applying a Gaussian window over 50 trials and with a standard deviation of 20. Black lines represent the peak(s) of previous or following events, relative to the event used as spatial filter, obtained from a rolling mean over 50 trials.

Simulating a ramp or a half-sine and expecting a half-sine

A: Representation of the averaged times and topographies for the three events detected by HMP on the downsampled real data. B: Trial and participant averaged topography (top panels) and time-courses (bottom panels) of the electrodes matched to the third HMP event for the ramp simulation (left panels) vs. the half-sine simulation (right panels). As in the main text, the trials were first centered on the peak of the third event before computing the average ERP for the different contrast values binned into 10 bins.