In order to quantify the choice repetition bias and its modulation by previous evidence, response time, and confidence and to qualitatively compare our findings to those of previous studies, we first employed a psychometric function fitting approach. We estimated choice repetition independently for each condition, following the analytical approach of earlier studies (
Akaishi et al., 2014;
Urai et al., 2017).
We first expressed the probability of a higher coherence response,
P(
rt = 1), as a function of the signed coherence difference between the reference and test stimulus (
\(\widetilde {{s_t}}\)) and fit a psychometric function (
Figure 2a) (
Wichmann & Hill, 2001) of the form
\begin{equation*}P\left( {{r_t} = 1|\widetilde {{s_t}}} \right) = \lambda + \left( {1 - 2\lambda } \right)g\left( {\delta + \alpha \widetilde {{s_t}}} \right)\end{equation*}
where
λ is the probability of stimulus-independent errors (lapses),
g is the logistic function,
α is perceptual sensitivity, and
δ is a bias term. The free parameters
λ,
α, and
δ were estimated by minimizing the negative log-likelihood of the data (using the MATLAB fminsearchbnd function).
For the quantification of serial choice bias, we first split the data into two bins corresponding to the previous choice such that one bin contained all trials for which the participant previously reported higher coherence, and the other bin contained all trials for which they reported lower coherence. For each level of previous absolute evidence strength (st–1) within these bins, we further split the data by previous response time (rt, based on a median split) or previous confidence (low ratings, 1 or 2; high ratings, 3 or 4). For each of those subsets of trials, we fit the psychometric function as described above. In order to compute the choice repetition bias, the resulting bias terms, δ, were transformed from log-odds into probabilities by the inverse logit function P = eδ/(1 + eδ). This probability reflects \(P( {{r_t} = 1|\widetilde {{s_t}} = 0} )\), which is the probability of choosing higher coherence in a hypothetical ambiguous trial (no evidence) in the current trial. To compute choice repetition probabilities, for each bin we averaged the probability to repeat the previous choice across the two previous choice options: \(p( {repeat} ) = ( {p({r_t} = 1|\;{r_{t - 1}} = 1,\;\widetilde {{s_t}} = 0} )\; + \;p({r_t} = 0|\;{r_{t - 1}} = 0,\;\widetilde {{s_t}} = 0))\;/\;2\). Finally, to test for differences in choice repetition probability, or p(repeat), across bins, we performed repeated-measures analyses of variance (ANOVAs) using SPSS Statistics 23 for Windows (IBM, Armonk, NY).