Trials in which the participant did not answer within time limit (3,000 ms) were excluded from the analysis. For each participant, we calculated the proportion of trials in which the participant reported that the reference stimulus was larger (or wider) than the standard. We then fitted two models to the data. The first two-parameter model was the sigmoid function
\(\frac{1}{{1 + {e^{ - \frac{{x - A}}{B}}}}}\). This model assumes that performance in the task solely relies on perceptual judgments. The second model included another “lapse rate” parameter. Following the logic of
Wichmann and Hill (2001), we applied a mixture model in which the lapse rate was added to the model as a free parameter. Addended trials followed a psychophysical function, whereas the responses in lapsed trials were determined by a “coin flip”—namely, randomly choosing between the left/right responses with equal probability. The mixture parameter
p—namely, the probability of lapse responses—was estimated as part of the model. Model comparison was based on the BIC (Bayesian information criterion) statistic, which penalizes for the number of parameters (BIC = –2 * log-likelihood + log(
N)*
K, where
N is the number of observations and
K is the number of parameters). Based on the sigmoid function of each model, we extracted the values of PSE, constant error (CE), JND, and GOF, being the squared correlation between observed and fitted values. The CE represents the magnitude of the illusion and was computed by subtracting the value of the PSE (50% “larger” responses) from the value of the standard stimulus. The JND represents a perceptual resolution to size differences in the context of the illusion and was calculated by dividing the range between 25% and 75% of the function by 2 (
Figure 1). For clarity, we transformed the CE and the JND raw scores to percentage scores for each participant, representing the magnitude of the illusion, and the magnitude of the JND in percentages compared to the standard stimulus. Reaction times (RTs) were also measured in each trial, and the mean RT was calculated for each participant and in each illusion.
A test–retest reliability was assessed by the correlation between the two sessions of each illusion. In addition, we calculated the correlations between the average CEs and JNDs in each illusion. To examine if the exposure time to the illusions affected the susceptibility to the illusion, we calculated the correlation between RTs (which also corresponds to exposure time in
Experiment 1) and between CEs for each task. Comparisons between the reliabilities of the CEs and JNDs of the illusions for nonoverlapping groups were performed using the “cocor” package in R (
Diedenhofen & Musch, 2015).