Simulation method: All distributions were generated in Matlab. For both no-distractor and distractor conditions, 10,000 SRTs were drawn randomly from a distribution in the Pearson system with mean equal to 183 ms, standard deviation equal to 52.83 ms, skewness equal to 0.78, and kurtosis equal to 3.89. These parameters were chosen to simulate the baseline distribution obtained empirically by Buonocore and McIntosh (
2008). Distributions (10,000 values) for neural delay, SI duration, and saccadic delay were generated by drawing values randomly from Gaussian distributions with means of 70, 70, and 62 ms and standard deviations of 10, 25, and 25 ms, respectively. In the distractor condition, the simulated SRTs were modulated following a simple set of rules. For each simulated saccade, an inhibitory window was defined with a lower limit set to the current Stimulus Onset Asynchrony (SOA) plus a value extracted from the neural delay distribution, and an upper limit set to the current lower limit plus a value extracted from the SI duration distribution. Each SRT falling within the inhibitory window had an associated value ranging from 0 to 1 taken from a discrete uniform distribution. If this value was smaller than the magnitude of inhibition (i.e., the probability to be inhibited; see below), a value extracted from the saccadic delay distribution was added to the current SRT. Otherwise the current SRT remained unchanged. In the simultaneous condition the upper limit was fixed at 155 ms after target onset. The probability of inhibiting a saccade was fixed at 0.28 (28% of the maximum inhibition). From the resulting set of simulated SRTs, percentage frequency histograms were created (bin width 4 ms) and lightly smoothed using a Gaussian kernel with 16-ms window and 2 ms
SD and then interpolated to obtain 1-ms precision. The absolute change for each point in time was computed by subtracting the baseline condition from the distractor condition (cf. Buonocore & McIntosh,
2008). The proportional change for each point in time in the distractor distribution relative to the no-distractor distribution was computed from the formula: (no-distractor − distractor)/no-distractor (cf. Bompas & Sumner,
2011). To simulate the effect of SOA, the SOA was incremented by 25 ms each loop, starting from a value of zero, for a total of seven iterations. Each simulation was run 100 times. At the end of the simulation the output of each condition was averaged across runs.