Abstract
Recent theories propose that schizophrenia/schizotypy and autistic spectrum disorder are related to impairments in Bayesian inference that is, how the brain integrates sensory information (likelihoods) with prior knowledge. However existing accounts fail to clarify: (i) how proposed theories differ in accounts of ASD vs. schizophrenia and (ii) whether the impairments result from weaker priors or enhanced likelihoods. Here, we directly address these issues by characterizing how 91 healthy participants, scored for autistic and schizotypal traits, implicitly learned and combined priors with sensory information. This was accomplished through a visual statistical learning paradigm designed to quantitatively assess variations in individuals’ likelihoods and priors. The acquisition of the priors was found to be intact along both traits spectra. However, autistic traits were associated with more veridical perception and weaker influence of expectations. Bayesian modeling revealed that this was due, not to weaker prior expectations, but to more precise sensory representations.
https://doi.org/10.7554/eLife.34115.001Introduction
In recent years Bayesian inference has come to be regarded as a general principle of brain function that underlies not only perception and motor execution, but hierarchically extends all the way to higher cognitive phenomena, such as belief formation and social cognition. Impairments of Bayesian inference have been proposed to underlie deficits observed in mental illness, particularly schizophrenia (Fletcher and Frith, 2009; Corlett et al., 2009; Adams et al., 2013; Hemsley and Garety, 1986; Friston, 2005; Stephan et al., 2006) and autistic spectrum disorder (ASD) (Pellicano and Burr, 2012a; Van de Cruys et al., 2014; Lawson et al., 2014; Palmer et al., 2017). The general hypothesis for both disorders is that the weight, also called ‘precision’, ascribed to sensory evidence and prior expectations is imbalanced, resulting in sensory evidence having relatively too much influence on perception.
In schizophrenia, overweighting of sensory information could explain the decreased susceptibility to perceptual illusions (Notredame et al., 2014), as well as the peculiar tendency to jump to conclusions (Speechley et al., 2010). Moreover, the systematically weakened lowlevel prior expectations might lead to forming compensatory strong and idiosyncratic highlevel priors (beliefs), which would explain the emergence and persistence of delusions as well as reoccurring hallucinations (Fletcher and Frith, 2009; Corlett et al., 2009; Adams et al., 2013).
In ASD, the relatively stronger influence of sensory information could explain hypersensitivity to sensory stimuli and extreme attention to details. The weaker influence of prior expectations would also result in more variability in sensory experiences. The desire for sameness and rigid behaviors could then be understood as an attempt to introduce more predictability in one’s environment (Pellicano and Burr, 2012a). Furthermore, this could lead to prior expectations which are too specific and which do not generalize across situations (Van de Cruys et al., 2014). While all theories agree that the relative influence of prior expectations is weaker in ASD, the primary source of this imbalance is debated: does it arise from increased sensory precision (i.e. sharper likelihood) or from reduced precision of prior expectations? (Brock, 2012; Pellicano and Burr, 2012b; van Boxtel and Lu, 2013) (Figure 1). Some authors argue for attenuated priors (Pellicano and Burr, 2012a; Pellicano and Burr, 2012b), while others argue for increased sensory precision (Lawson et al., 2014; Palmer et al., 2017; Brock, 2012; Van de Cruys et al., 2013) but conclusive experimental evidence is lacking.
A number of studies have aimed at testing Bayesian theories, either in a clinical population, or by studying individual differences in the general population (Powell et al., 2016; Skewes et al., 2015; Teufel et al., 2015; Schmack et al., 2013) under the hypothesis of a continuum between autistic/schizotypal traits and ASD/schizophrenia (Nelson et al., 2013; van Os et al., 2009; Constantino and Todd, 2003).
Attenuated slowspeed priors were reported in a motion perception task in individuals with ASD traits (Powell et al., 2016). Autistic children also showed attenuated central tendency prior in temporal interval reproduction (Karaminis et al., 2016). Attenuated priors were also reported in perceptual tasks that incorporate probabilistic reasoning (Skewes et al., 2015; Skewes and Gebauer, 2016). However, the direction of gaze priors (Pell et al., 2016) and the lightfromabove priors (Croydon et al., 2017) were found to be intact. Autistic children also demonstrated intact ability to update their priors in a volatile environment in a decisionmaking task (Manning et al., 2017) but a followup study in ASD adults showed that they overestimate volatility in a changing environment (Lawson et al., 2017).
In schizophrenia/schizotypal traits, Teufel et al. (2015) reported increased influence of prior expectations when disambiguating twotone images, while Schmack et al. (2015, 2017) reported weakened influence of stabilizing predictions when observing a bistable rotating sphere.
Overall, the existing findings are not only mixed, but also employ very different paradigms, which makes their direct comparison difficult. Further, a critical limitation of most studies (except for Karaminis et al., 2016) is the lack of formal computational models that can test whether behavioral differences originate from different priors or from different likelihoods. Moreover, to our knowledge, despite the similarity of the Bayesian theories proposed for ASD and schizophrenia, there is no previous work investigating both autistic and schizotypal traits within the same experimental paradigm so as to test their differences.
We here address these questions empirically in a context of visual motion perception. We used a previously developed statistical learning task (Chalk et al., 2010) in which participants have to estimate the direction of motion of coherently moving clouds of dots (Figure 2). Chalk et al. (2010) found that in this task healthy participants rapidly and implicitly develop prior expectations for the most frequently presented motion directions. This in turn alters their perception of motion on low contrast trials resulting in attractive estimation biases towards the most frequent directions. In addition, prior expectations lead to reduced estimation variability and reaction times, as well as increased detection performance for the most frequently presented directions. When no stimulus is presented, the acquired expectations sometimes lead to false alarms (hallucinations), again, mostly in the most frequent directions. Importantly, such biases were well described using a Bayesian model, where participants acquired a perceptual prior for the visual stimulus that is combined with sensory information and influences their perception. As such, this paradigm is well suited to quantitatively model variations in likelihoods and priors in individuals with ASD or schizotypal traits.
Results
Here, we investigated individual differences in statistical learning in relation to autistic and schizotypal traits in a sample of 91 healthy participants. Eight participants failed to perform the task satisfactorily and were excluded from the analysis (see Materials and methods), leaving 83 participants in the study (41 women and 42 men, age range: 18–69; mean: 25.7).
Task behavior at low contrast
First, we investigated whether participants acquired priors on the group level. We discarded the first 170 trials as that is how long it took for the 2/1 and 4/1 staircases contrast levels to converge (Appendix 1—figure 2) and for prior effects to become significant (Appendix 1—figures 3, 4 and 5). We analyzed task performance at low contrast levels (converged 2/1 and 4/1 staircases contrast levels) where sensory uncertainty is high. Replicating findings of Chalk et al. (2010), we found that on the group level people acquired priors that approximated the statistics of the task. Such priors were indicated by: attractive biases towards ±32° (Figure 3a), less variability in estimations at ±32° (Figure 3b; standard deviation of estimations 11.9 ± 0.30° at ±32° versus 13.84 ± 2.38° over all other motion directions; signed rank test: p<0.001), shorter estimation reaction times at ±32° as compared to all other motion directions (Figure 3c; average reaction time was 201.87 ± 2.47 ms at ±32° versus 207.75 ± 2.60 ms over all other motion directions; signed rank test: p<0.001) and better detection at ±32° as compared to all other motion directions (Figure 3d; detected 75.57 ± 0.65% at±32° versus 66.70 ± 0.83% over all other motion directions; signed rank test: p<0.001).

Figure 3—source data 1
 https://doi.org/10.7554/eLife.34115.005
Nostimulus performance
Another indicator of acquired priors is the distribution of estimation responses on trials when no actual stimulus was presented. We found that participants sometimes still reported seeing dots (experienced hallucinations) but mostly so around ±32° (Figure 3f, solid line). To quantify the statistical significance of hallucinations around ±32°, the space of possible motion directions was divided into 45 bins of 16° and the probability of estimation within 8° of ±32° was multiplied by the total number of bins:
where N_{bins} is the number of bins (45), each of size 16°. This probability ratio would be equal to one if participants were equally likely to estimate within 8° of ±32°, as they were to estimate within other bins. We found that the median of p_{rel} was significantly greater than 1 (median(p_{rel})=1.6, p<0.001, signed rank test). Furthermore, the estimation distribution when no dots where detected (Figure 3f, dashdot line) was found to be significantly flatter (median(p_{rel})=0, p<0.001, signed rank test comparing with the median of p_{rel} for hallucinations), suggesting that the hallucinations were indeed of perceptual nature (rather than related to a response bias).
Task performance and autistic/schizotypy traits
Participants were prescreened to make sure they covered a wide range of autistic and schizotypy scores. The AQ scores in our sample ranged from 6 to 41 with a mean (±SD) of 20.3 (±8.3). The RISC scores ranged from 8 to 55 with a mean of 31.7 (±11.9), and the SPQ scores ranged from 4 to 59 with a mean of 26.4 (±13.8).
We found that on low contrast trials autistic traits lead to less variability in estimations (Figure 4b; mean standard deviation of estimations: r = −0.327, p<0.001), which remained significant after Bonferroni correction (p=0.002). Moreover, there was a negative relationship between autistic traits and estimation bias, which was trending according to robust regression (Figure 4a; mean absolute estimation bias: r = −0.175, p=0.053) and significant according to Kendall’s correlation (τ_{b} = −0.163, p=0.032), however, it did not survive Bonferroni correction (p=0.212). In the Bayesian framework, less bias could arise either due to wider priors or narrower sensory likelihoods, while less variability could be a result of either narrower priors or narrower likelihoods (see Figure 1). Thus, observing less bias and less variability together suggests that the effects are driven by narrower likelihoods. An alternative is that the differences in variability could be due to differences in motor precision, which we further assess via modeling (below).

Figure 4—source data 1
 https://doi.org/10.7554/eLife.34115.007
Schizotypy traits (RISC and SPQ scores) did not show any effect on task performance at low contrast as indicated by the absence of correlations with mean absolute estimation bias (RISC: r = 0.140, p=0.197; SPQ (N = 39): r = −0.160, p=0.204) and with mean estimation variability (RISC: r = 0.197, p=0.092; SPQ (N = 39): r = −0.229, p=0.171); see Appendix 1—figures 6, 7 and 8.
Nostimulus trials and autistic/schizotypal traits
We also investigated how the traits affected performance on trials when no actual stimulus was presented. First, we looked at the total number of estimations. We found that autistic traits were associated with less hallucinations (Figure 4c; r = −0.238, p=0.010), while schizotypal traits were found to have no effect on the number of hallucinations (RISC: r = 0.126, p=0.163; SPQ (N = 39): r = −0.010, p=0.959). Secondly, we looked for relationships between the traits and how the estimations on nostimulus trials were distributed. Specifically, we were interested in whether the traits predicted how densely hallucinations were distributed around ±32°, as this could be considered to reflect the differences in the width of the underlying acquired prior distribution. For weaker priors we would expect a more spread out distribution of hallucinations. To test this hypothesis, we looked at the fraction of total hallucinations in the region around ±32° for three differentsized windows: Within 8°, within 16° and within 24° of ±32°. Bayesian Kendall correlation analysis on these measures provided positive evidence that none of the traits had any effect on how hallucinations were distributed, suggesting no differences in the acquired prior distributions (fraction of hallucinations within 8° of ±32°: AQ  τ_{b} = 0.003, BF_{01} = 7.24; RISC  τ_{b} = 0.050, BF_{01} = 3.73; SPQ  τ_{b} = 0.101, BF_{01} = 8.72; within 16° of ±32°: AQ  τ_{b} = 0.068, BF_{01} = 2.86; RISC  τ_{b} = 0.129, BF_{01} = 0.84; SPQ  τ_{b} = 0.018, BF_{01} = 5.45; within 24° of ±32°: AQ  τ_{b} = 0.057, BF_{01} = 11.67; RISC  τ_{b} = 0.078, BF_{01} = 2.40; SPQ  τ_{b} = 0.006, BF_{01} = 5.02).
Modeling results
Group level results
To quantitatively evaluate the relationships between underlying perceptual mechanisms and task performance we fitted a range of generative models. One class of models was Bayesian  it was based on the assumption that participants combine prior expectations with uncertain sensory information on a single trial basis (Figure 5).
To account for the possibility that the bimodal probability distribution of the stimuli, in addition to inducing prior expectations, has also affected the sensory likelihood, we constructed three variations of the Bayesian model: BAYES, where the sensory precision was constrained to be the same across all presented motion directions, 'BAYES_varmin, where the sensory precision was allowed to be different for the most frequently presented motion directions, but was the same across all other directions, and BAYES_var where sensory precision was allowed to be different across all motion directions. Another class of models was based on the assumption that task performance can be explained by response strategies that do not involve Bayesian inference. That is, on any given trial participants responded based on the prior expectations or sensory information alone. We considered four variations of response strategy models: ADD1, ADD2, ADD1_m, and ADD2_m (see Methods for details).
To compare the models, we computed BIC values for each individual for each model; we used individual BIC values as a summary statistic and compared the models using signed rank test in order to preserve individual variability, which corresponds to a random effects Bayesian model selection procedure. We found that the BAYES model had significantly smaller BIC values than the remaining models (see the pvalues within Figure 6a).
To determine how the best fitting model compared to the actual data, we analyzed the estimation biases and variation in estimation responses as predicted by BAYES (Figure 6b,c). As in the experimental data analysis, we computed estimation distributions predicted by the model by assuming occasional random estimations (see Equation 2). Finally, using the BAYES model, we reconstructed the priors acquired by participants. While on the individual level there was a considerable variation in the shape of acquired priors (see Appendix 1—figure 10), on the group level, it approximated the statistics of the task (Figure 6d).
Model parameters and autistic/schizotypal traits
Correlational analysis of BAYES model parameters showed that there was no correlation between AQ and the precision of the prior σ_{exp} (Figure 7b; r = 0.018, p=0.962). That autistic traits had no effect on the precision of the prior was confirmed by Bayesian Kendall correlation, which provided positive evidence (τ_{b} = 0.001, BF_{01} = 6.99).

Figure 7—source data 1
 https://doi.org/10.7554/eLife.34115.011
Importantly, autistic traits were found to be strongly associated with less uncertainty in the sensory likelihood, σ_{sens} (Figure 7c; r = −0.185, p=0.011), which also remained significant after Bonferroni correction (p=0.044). Finally, there was no correlation with the amount of random estimations (Figure 7d; r = −0.135, p=0.238). Motor precision, which was estimated from high contrast trials, separately from all other parameters (see Methods), was also correlated with autistic traits (r = 0.245, p=0.012). On the other hand, consistent with the absence of differences in the behavioral findings, schizotypal traits were not associated with any difference in the BAYES model parameter values (Appendix 1—figure 9), and in particular, were found to have no effect on prior precision (RISC: τ_{b} = 0.012, BF_{01} = 6.90; SPQ: τ_{b} = 0.071, BF_{01} = 3.97).
Parameter recovery for BAYES
Finally, to further investigate that in our experimental paradigm the influence of stronger likelihoods can be distinguished from that of weaker priors (Brock, 2012; Pellicano and Burr, 2012b) we performed parameter recovery for the winning BAYES model. Parameter recovery involves generating synthetic data with different sets of parameters ('actual parameters') and then fitting the same model to estimate the parameters ('recovered parameters') that are most likely to have produced the data. If actual and recovered parameters are in a good agreement, it means that the effects of different parameters can be reliably distinguished. At the same time, parameter recovery is also affected by the parameter estimation methods and even more so by the amount of data used for model fitting. Therefore, parameter recovery provides an overall check for the reliability of modelling results and is recommended as an essential step in computational modelling approaches (Palminteri et al., 2017).
We found that overall BAYES model (and MLE parameter estimation using simplex optimization function) recovered parameters very well, which was reflected in Pearson’s correlation between actual and recovered estimates being r > 0.9 for all model parameters (Figure 8).
Discussion
In this study, we investigated whether autistic and schizotypal traits are associated with differences in the implicit Bayesian inference performed by the brain. Specifically, we wanted to know whether autistic and schizotypal traits are accompanied by (1) differences in how the priors are updated and/or in their precision and/or by (2) differences in the precision with which the sensory information (the likelihood) is represented. We used a visual motion estimation task (Chalk et al., 2010) that induces implicit prior expectations via more frequent exposure of two motion directions (±32°). We found that on the group level (N = 83) participants acquired prior expectations towards ±32° motion directions. This was indicated by shorter estimation reaction times and better detection at ±32°, as well as attractive biases towards ±32° and reduced estimation variability at ±32°. Moreover, when no stimulus was presented, participants sometimes still reported seeing the stimulus, mostly around ±32°. Performance was best explained by a simple Bayesian model, which provided a good fit to the data and captured the characteristic features of perceptual bias and variability. This model provided estimates of Bayesian priors and sensory likelihoods for each participant, which were then analyzed in relation to participants’ schizotypal and autistic traits.
Schizotypal traits were found to have no measurable effect on perceptual biases in our task and, therefore, were not associated with any differences in the precision ascribed to priors and likelihoods. This finding challenges recent accounts of positive symptoms of schizophrenia that predict impaired updating of priors and an imbalance in precision ascribed to sensory information and prior expectations (Fletcher and Frith, 2009; Corlett et al., 2009; Adams et al., 2013). An immediate explanation might be that the influence of schizotypal traits in the healthy population is not strong enough to lead to behavioral differences, even if the dimensionality assumption holds. This would need to be addressed by further research investigating clinical populations. Another possibility is that the aberrant perception subconstruct of schizotypal traits, for which we did not acquire explicit measures, is more relevant for the hypothesized effects then the entire construct as a whole. For example, a recent study by Powers et al. (2017) found that overweighing of perceptual priors was specifically linked to hallucinatory propensity and not to the diagnostic status of psychosis itself. Furthermore, Teufel et al. (2015) also found that stronger influence of prior knowledge was primarily associated with hallucinatory propensity and not with delusional propensity. Another possible difference between Teufel et al. (2015) study and ours might be the level at which the priors operate. In Teufel et al. (2015) participants were presented with ambiguous twotone versions of images before and after seeing the actual images in full color and had to report whether the presented twotone image contains a face. The lowlevel prior for basic perceptual features (as induced in our task) might function at a hierarchically lower level than prior knowledge related to complex collection of features and semantic content (faces). The level at which prior expectations are induced has indeed been shown to matter. A series of studies by Schmack et al. (2013, 2015, 2017) using 3D rotating cylinders report weaker lowlevel (perceptuallyinduced  stabilizing) priors but stronger highlevel (cognitivelyinduced) priors in both schizophrenia and schizotypal traits. It is difficult to compare and reconcile these findings with ours. One possibility is that the priors induced in our task lie in between their perceptual and cognitive levels. The taxonomy of priors in relation to their place in the computational hierarchy or to their complexity or specificity is still far from being established (Seriès and Seitz, 2013) and thus the potential relevance of such distinctions is still not known.
Autistic traits were associated with significant behavioral differences: weaker biases and lower variability of direction estimation on low contrast trials. Modeling revealed that this was because of increased sensory precision as well as higher motor precision, while there was no attenuation of acquired priors. Parameter recovery analysis confirmed that our methodology provides reliable parameter estimates and, in particular, allows disentangling variations in priors and likelihoods.
Autistic traits were also found to be associated with less false detections (hallucinations) on trials when no stimulus was presented, consistent with the idea that prior expectations had less influence in individuals with higher AQ. In an attempt to measure those individual differences, we fitted a more sophisticated Bayesian model that could account not only for the estimation performance but also for the detection data (see Appendix 2). This model provided a good fit to both estimation and detection data, and preserved the correlation between ASD traits and the precision of the motion direction likelihood (r = −0.202, p=0.029). However, parameter recovery was not as good as for the BAYES model presented above (see Appendix 2—figure 3) and for this reason we focused on the simpler model in this paper.
Overall, our findings are in agreement with most of the recent Bayesian theories of ASD, namely, that autistic traits are associated with a relatively weaker influence of prior expectations. However, we find that this is due to enhanced sensory precision (Lawson et al., 2014; Palmer et al., 2017; Brock, 2012; Van de Cruys et al., 2013), rather than attenuated priors per se (Pellicano and Burr, 2012a). Other empirical studies inspired by the Bayesian accounts have reported either attenuated or intact priors, but most are subject to methodological limitations, either because they did not use computational modeling (Skewes et al., 2015; Skewes and Gebauer, 2016; Croydon et al., 2017) or because their model could not extract likelihoods and quantify their variations (Powell et al., 2016; Lawson et al., 2017).
The idea that sensory processing could be enhanced in autism has long been proposed outside the Bayesian framework. Autistic traits have been associated with enhanced orientation discrimination (Dickinson et al., 2014), but only for firstorder (luminancedefined) stimulus (Bertone et al., 2005). This enhancement has been proposed to be a result of either enhanced lateral (Bertone et al., 2005), or a failure to attenuate sensory signals via topdown gain control (Lawson et al., 2014), both of which could be directly related to narrower likelihoods in the Bayesian framework (Ma et al., 2006). However, in motion perception, previous research did not find improved discrimination for firstorder stimulus in autism, while for secondorder (texturedefined) stimulus, the autistic group was found to underperform (Bertone et al., 2003). Our findings challenge these results and call for more research in this area.
In ASD as in schizotypy, prior integration might function differently at different levels of sensory processing. For example, Pell et al. (2016) reported intact directionofgaze priors for healthy individuals with high autistic traits and for highly functional individuals with a clinical diagnosis. The authors did not directly investigate differences in sensory precision, but the lack of behavioral differences suggests that there was none. Arguably, their paradigm involves more complex stimuli than used in our task, which are also strongly associated with semantic content (faces). It would not be surprising if increased sensory precision does not extend to such stimuli. In fact, autistic individuals are known to exhibit differential performance based on the complexity of the stimulus (Bertone et al., 2005), which also lies at the foundation of some theoretical accounts, such as the ‘Weak Central Coherence’ (Happé and Frith, 2006).
In our paradigm people acquire prior expectations very quickly, within 200 trials (see Appendix 1), which did not allow us to study individual differences in the rate at which the priors are acquired. Bayesian accounts predict differences in the dynamical updating of the priors, namely, that both autistic and schizotypal traits should be associated with increased learning rate  which is the ratio of likelihood and posterior precisions (Palmer et al., 2017). Our findings of increased sensory precision in autistic traits also suggest that their learning rate should be faster. However, this prediction might need to be more nuanced for volatile environments when there are multiple (hierarchical) levels of uncertainty that need to be updated simultaneously. A recent study by Lawson et al. (2017) found that when transitioning from stable to volatile environments, autistic adults showed larger change in the learning rate about volatility and smaller change in the learning rate about the environmental probabilities, while the average learning rates were found to not be different from those of controls.
Another aspect that our paradigm could not test is the specificity of the acquired priors (Seriès and Seitz, 2013). Some Bayesian accounts (Van de Cruys et al., 2014) predict that priors may be overly contextsensitive in autism. This is in line with the view that generalization is impaired in autism (Plaisted, 2015). Furthermore, such overspecificity is thought to be stronger with more repetitive stimuli (Harris et al., 2015). Future research could address this using statistical learning paradigms that incorporate increasingly distinct contexts or stimuli.
Conclusion
We investigated statistical learning and Bayesian inference in a visual motion perception task along autistic and schizotypal traits. To our knowledge, this study is the first to investigate differences in Bayesian inference along both trait spectra in a single task. Furthermore, this study is the first visual study to computationally disentangle and quantitatively assess the variations in individuals’ likelihoods and priors. Surprisingly, schizotypal traits were found to have no effect on task performance and thus were not associated with any differences in the underlying statistical learning and Bayesian inference. For autistic traits, however, significant behavioral differences in prior integration were found, which were due to an increase in the precision of internal sensory representations in participants with higher AQ. Whether the current results extend to clinical populations will have to be examined in the future.
Materials and methods
Participants
91 (47 females, 44 males, age range: 18–69) naïve participants with no motor disabilities and with normal (or corrected to normal) vision were recruited from the general population. We advertised for participants using posters and the internet across University of Edinburgh locations and other sites across Edinburgh. All participants gave informed written consent and received monetary compensation for participation. The study was approved by the University of Edinburgh School of Informatics Ethics Panel.
Questionnaires
Request a detailed protocolASD was assessed using 50item version Autism Spectrum Quotient (AQ) (BaronCohen et al., 2001), which is commonly used for assessing milder variants of autisticlike traits within the general population. Schizotypal traits were assessed using The Rust Inventory of Schizotypal Cognitions (RISC) (Rust, 1988). RISC is specifically developed to measure schizotypal traits in the general population. In addition, a subgroup of 41 participants also completed Schizotypal Personality Questionnaire (SPQ) (Raine, 1991). Finally, all participants were also asked to complete the WarwickEdinburgh Mental Wellbeing Scale (WEMWBS) (Tennant et al., 2007) in order to control for potential depressioninduced differences in performance (Austin et al., 2001).
Apparatus
The visual stimuli were generated using Matlab Psychophysics Toolbox (Brainard, 1997). Participants viewed the display in a dark room at a distance of 80–100 cm. The stimuli consisted of a cloud of dots with a density of 2 dots/deg^{2} moving coherently (100%) at a speed of 9°/sec. Dots appeared within a circular annulus with minimum diameter of 2.2° and maximum diameter of 7°. The stimuli were displayed on a Dell P790 monitor running at 1024 × 768 at 100 Hz. The display luminance was calibrated using a Cambridge Research Systems Colorimeter (ColorCal MKII).
The task
Request a detailed protocolThe task was developed previously in our laboratory (Chalk et al., 2010). Participants have to: (i) estimate the direction of coherently moving simple stimuli (dots) that are presented at low contrast levels (estimation task) and then (ii) indicate whether they have actually perceived the stimulus or not (detection task). Since Chalk et al. (2010) had shown that the effects of acquired priors become significant within the first 200 trials, instead of two experimental sessions of 850 trials each as in the original study, we used a single session of 567 trials (lasting around 40 min).
Each trial started by first displaying a fixation point (0.5°, 12.2 cd/m^{2}) for 400 ms, after which a field of moving dots appeared along with an orientation bar (length 1.1°, width 0.03°, luminance 4 cd/m^{2}, extending from the fixation point). Initial angle of the bar was randomized for each trial. Participants had to estimate the direction of motion by aligning the bar (using a computer mouse) to the direction the dots were moving in, and by clicking the mouse button to validate their estimate. The display cleared when either the participant had clicked the mouse or when 3000 ms had elapsed. On trials where no stimulus was presented, the bar still appeared for the estimation task to be completed.
After a 200 ms delay, the participants had to indicate whether they had actually detected the presence of dots in the estimation period (detection task). The display was divided into two parts by a vertical white line across the center of the screen, the left hand side area reading ‘NO DOTS’ and the right hand side area reading ‘DOTS’ (Figure 2a). The cursor appeared in the center of the screen, and participants had to move it to the left or right and click to indicate their response. Immediate feedback for correct or incorrect detection responses was given by a cursor flashing green or red, respectively. The screen was cleared for 400 ms before the start of a new trial. Every 20 trials, participants were presented with feedback on their estimation performance in terms of average estimation error in degrees (e.g., ‘In the last 20 trials, your average estimation error was 23°'). Every 170 trials (i.e. on three occasions) participants were given a chance to ‘have a short break to rest their eyes’, in order to prevent fatigue. Participants clicked when they were ready to continue.
Design
Request a detailed protocolThe stimuli were presented at four different levels of contrast: 0 contrast (nostimulus trials), two low levels contrasts and high contrast, randomly mixed across trials. There were 167 trials with no stimulus. The two low levels of contrast were determined using 4/1 and 2/1 staircases on detection performance (GarcíaPérez, 1998). There were 243 trials following the 4/1 staircase and 90 trials following the 2/1 staircase. The remaining 67 trials were at high contrast, which was set to 3.51 cd/m^{2} above the background luminance.
For the two low contrast levels, there was a predetermined number of possible directions: 0°, ±16°, ±32°, ±48°, and ±64° with respect to a reference direction. The reference direction was randomized for each participant. For the 2/1 staircased contrasts, each predetermined motion direction was presented equally frequently. Unbeknownst to participants, stimuli at high and 4/1 staircase contrasts were presented more frequently at −32° and +32° motion directions, resulting in a bimodal probability distribution (Figure 1b). For the 4/1 staircase contrast level, the dots were moving at ±32° in 173 (∼70%) trials and in all the other predetermined motion directions in the remaining 70 (∼30%) trials equally frequently. At the highest contrast level, 34 (∼50%) trials had the dots moving at ±32° and the remaining 33 (∼50%) trials were at random directions (i.e. not just the predetermined directions).
Data analysis
Request a detailed protocolResponses on high contrast trials were used as a performance benchmark to ensure that participants were performing the task adequately. The predefined inclusion criteria were: (1) at least 80% detection and (2) less than 30° root mean squared error of estimations. 8 out of 91 participants failed to satisfy at least one of the criteria and were excluded from further analysis (Appendix 1—figure 1).
Data analysis on the estimation of motion directions was performed on 4/1 and 2/1 staircased contrast levels only and only on trials where participants both validated their choice with a click within 3000 ms in the estimation part and clicked ‘DOTS’ in the detection part. The first 170 trials of each session were excluded from the analysis, as this was the upper limit for the convergence of the staircases to stable contrast levels (Appendix 1—figure 2).
After removing these trials, the luminance levels achieved by the 2/1 and 4/1 staircases were found to be considerably overlapping (Appendix 1—figure 2). Therefore, the data for both of these contrast levels was combined for all further analysis.
To account for random estimations (either accidental or intentional) that participants made on some trials, we fitted each participant’s estimation responses to the probability distribution:
Where α is the proportion of trials in which participant makes random estimates, and V(θµ,κ) is the probability density function for the estimated angle θ for von Mises (circular normal) distribution with the mean µ and precision κ. The parameters µ and κ of the von Mises distribution were determined by maximizing the likelihood of the distribution in Equation (2) for each presented angle.
To analyze the distribution of estimations in nostimulus trials, we constructed histograms of 16° size bins. These histograms were converted into probability distributions by normalizing over all motion directions. We analyzed the estimation distribution when participants reported seeing dots (clicked ‘DOTS’) within nostimulus trials. We interpreted these false alarms as a simple form of perceptual hallucination.
Modelling
Bayesian models
Request a detailed protocolBayesian models assume that participants combined a learned prior of the stimulus directions with their sensory evidence in a probabilistic manner. We first assume that participants make noisy sensory observations of the actual stimulus motion direction (θ_{act}), with a probability
where θ_{t} itself varies from trial to trial around θact according to p(θ_{t}θ_{act})=V(θ_{act}, κ_{sens}).
While participants cannot access the ‘true’ prior, p(θ), directly, we hypothesized that they learned an approximation of this distribution, denoted p_{exp}(θ). This distribution was parameterized as the sum of two von Mises distributions, centered on motion directions θ_{exp} and θ_{exp}, and each with precision κ_{exp}:
Combining these via Bayes’ rule gives a posterior probability that the stimulus is moving in a direction θ:
The perceived direction, θ_{perc}, was taken to be the mean of the posterior distribution (almost identical results would be obtained by using the maximum instead). Finally, we accounted for motor precision and a possibility of random estimates on some trials via:
where α is the proportion of trials in which participants make random estimates and κ_{m} is the motor precision.
Increased exposure to some motion directions might not only give rise to prior expectations, but also induce learning in the sensory likelihood function itself (Stocker and Simoncelli, 2006; Sato and Kording, 2014). Therefore, we fitted two more model variants: 'BAYES_var’ where κ_{sens} varied with the stimulus direction (i.e. it took five different values for each of the angles: 0°, ±16°, ±32°, ±48°, ±64°) and 'BAYES_varmin’ where κ_{sens} was allowed to be different for ±32° but was the same for all other directions.
Response strategy models
Request a detailed protocolWe wanted to test whether task behavior might be better explained by simple behavioral strategies. This class of models assumed that on trials when participants were unsure about the presented motion direction, they made an estimation based solely on prior expectations, while on the remaining fraction of trials they made unbiased estimates based solely on sensory inputs. The first model, 'ADD1’, assumed that estimations derived from prior expectations were simply sampled from a learnt expected distribution, p_{exp}(θ) (see Chalk et al., 2010 and Appendix 2). The second model, 'ADD2’, was just as 'ADD1’ except when participants were unsure about the stimulus motion direction, instead of sampling from the complete learned probability distribution ranging from −180° to +180°, they effectively truncated this distribution on a trial by trial basis and sampled from only one part of it, negative (−180° to 0°) or positive (0° to +180°), depending on which side of the distribution the actual stimulus occurred (see Chalk et al., 2010) and Appendix 2). We also considered slight variations of the ‘ADD1’ and ‘ADD2’ models, denoted ‘ADD1_m’ and ‘ADD2_m’ respectively. These were identical to ‘ADD1’ and ‘ADD2’ except from setting 1/κ_{exp} to zero; that is, on trials when perceptual estimates were derived only from expectations, they were equal to the mode of the learnt distribution (i.e. no uncertainty).
Parameter estimation
Request a detailed protocolWe used performance in high contrast trials to estimate motor precision, κ_{m}, for each individual. We assumed that, for those trials, sensory uncertainty was close to zero. Motor precision was then determined by fitting estimation responses to the distribution in Equation (2) by replacing µ with the actual motion direction, θ_{act}. The estimated motor precision was used in all subsequent model fitting as a fixed parameter. The rest of the free parameters were estimated by fitting the response data at the two low (staircased) contrast levels. For each model with a set of free parameters M, we computed the probability distribution p(θ_{est}θ_{act}; M) of making an estimate θ_{est} given the actual stimulus direction θ_{act}. For the response strategy models, by definition, the p(θ_{est}θ_{act}; M) corresponds to average behavior in the task.
The parameters were estimated by maximizing the fit of the log likelihood function for the experimental data for each participant individually. The maximum likelihood was found using a simplex algorithm, using fminsearchbnd Matlab function. To avoid convergence at a local maximum we constructed a grid of initial κ_{exp} and κ_{sens} parameter values covering the range found in previous studies. We selected the resulting set of parameters that corresponded to the largest loglikelihood.
Model comparison
Request a detailed protocolTo compare the model fits we used Bayesian Information Criterion (BIC), which approximates the log of model evidence (Burnham and Anderson, 2004):
where M is model, D is observed data and P (DM, $\hat{\mathrm{\Theta}}$) is the likelihood of generating the experimental data given the most likely set of parameters, $\hat{\mathrm{\Theta}}$; k is the number of model parameters and n is the number of data points (or equivalently, the number of trials). BIC evaluates the model by how it fits the data by also penalizing for model complexity (number of parameters); lower BIC score indicates a better model.
Parameter recovery
Request a detailed protocolTo determine whether the BAYES model can distinguish the effects of strong likelihoods from those of weak priors (Brock, 2012; Pellicano and Burr, 2012b) and to evaluate the robustness of our methods, we performed parameter recovery. First, we generated 80 sets of parameters (i.e. 80 synthetic individuals) by randomly sampling each parameter from a Gaussian distribution centered on the mean value of each parameter found in our sample (40° for θ_{exp}, 15° for σ_{exp}, 10° for σ_{sens}, 0.06 for α and 10° for σ_{motor}). Second, for each set of parameters, we simulated data for 200 trials with the Bayesian model by randomly sampling from the estimation probability distribution. We used 200 simulated trials only, to match the empirical data (200 corresponds to the amount of experimental trials used for fitting, after excluding high contrast and zero contrast trials; Simulating more trials would result in a better parameter recovery but the results would no longer be informative about the reliability of parameters estimated from empirical data). Finally, we fitted the BAYES model to the simulated data. To evaluate the goodness of recovered parameters, we computed Pearson’s correlation between the actual parameters and the recovered parameters.
Statistical tests
Request a detailed protocolDue to the presence of outliers in many of the measures, we used robust regression techniques for measuring the presence and strength of the effects in our data. This was done using robustfit function in Matlab, which downweighs the influence of outliers in proportion to their distance from the regression line, which is computed via iteratively reweighted least squares (IRLS) (Holland and Welsch, 1977). For the loss function we used Huber function (Huber, 1964) with a tuning constant of 1.345, which corresponds to 95% estimator efficiency as compared to ordinary least squares.
Furthermore, we applied Bonferroni correction for multiple testing based on the number of independent hypotheses that we tested; that is, whether two personality traits, ASD and schizotypy, were associated with the two variables of interest, acquired priors and sensory likelihoods,  this resulted in four different hypotheses. Note that while the number of null hypothesis significance tests that we performed exceeds this number, the tests within each set concerning the same hypothesis were not independent (each test was based on derivative and/or correlated values to those in the other tests within the same set), and thus would not have met the independence assumption on which Bonferroni correction is based.
Finally, due to the limitations of frequentist statistics for accepting the null hypothesis, we performed Bayesian correlation analysis and computed Bayesian Factors (Kass and Raftery, 1995) for the null hypothesis (BF_{01}). This was done using JASP (Team, 2017) (Version 0.8.6). Due to the presence of outliers, this analysis was carried out using the nonparametric Kendall’s Taub correlation coefficient.
Source code and data
Request a detailed protocolThe source data of the main figures is provided. These include, Figure 3—source data 1, Figure 4—source data 1 and Figure 7—source data 1. Source code 1 contains all the source code necessary to reproduce the figures. More detailed information about the source code is in SourceCode_Readme.txt, while SourceData_Readme.txt contains more details about the source data files.
Appendix 1
Exclusion criteria
In order to ensure that participants performed adequately in the psychophysical task, we used predetermined performance criteria for inclusion into the study. Firstly, participants were required to detect the motion stimuli on more than 80% of trials with the high contrast motion stimuli and also make active estimates of the motion directions by clicking the mouse. Secondly, their average estimation performance on the high contrast stimuli had to be within 30° of the correct angle. 8 out of 91 participants failed to satisfy at least one of the criteria: two participants did not satisfy the first criteria, four did not satisfy the second criteria and two did not satisfy both of the criteria (Appendix 1—figure 1). These participants were excluded from further analysis.
Staircased stimulus contrast levels
Appendix 1—figure 2 describes the average convergence of the contrast staircases. Two groups comprising our sample performed the task at different background contrast levels. For a subgroup of 50 participants (left panel), the background luminance was set to 1.16 cd/m^{2} for the other subgroup of 41 (right panel) it was set to 5.18 cd/m^{2}. For both groups, contrast staircases converged after 170 trials for both intermediate contrast levels, denoted with the vertical dashed line. In both groups, 2/1 and 4/1 staircased contrasts were considerably overlapping: on average 2/1 being 0.20 ± 0.04 cd/m^{2} and 4/1 being 0.22 ± 0.04 cd/m^{2} above the 1.16 cd/m^{2} background luminance; and on average 2/1 being 0.42 ± 0.05 cd/m^{2} and 4/1 being 0.46 ± 0.05 cd/m^{2} above the 5.18 cd/m^{2} background luminance. Thus, the two intermediate contrasts were combined for all further data analysis.
Combining the different background luminance levels
To compare the two subgroups that performed the task at different background luminance levels, we performed Wilcoxon twotailed rank sum test for all of the behavioral measures and none of them indicated any differences: mean absolute estimation bias (z = 0.652; ranksum = 1920; p=0.514), mean variance of estimations (z = −0.406; ranksum = 1803; p=0.685), total number of hallucinations (z = 0.128; ranksum = 1862; p=0.898) number of hallucinations within 8° of ±32° (z = 0.870; ranksum = 1943; p=0.384), mean estimation reaction time (z = 0.479; ranksum = 1901; p=0.632). The two groups were therefore combined.
Temporal emergence of the impact of expectations
We investigated how many trials it took for the acquired prior effects to impact behavior. First, we looked at estimation reaction times (RT) and compared mean RT of each individual at ±32^{°} with mean RT at all other directions; we compared cumulative moving averages at every 30 trials (Appendix 1—figure 3). We found that it took less than 90 trials for RT at ±32^{°} to become significantly shorter than average RT at all other directions (Appendix 1—figure 3 and pvalues within).
Similarly, we looked at average detection performance and compared the fraction of trials in which stimulus was detected at ±32^{°} with the mean fraction detected over all other presented directions; again, we compared cumulative moving averages at every 30 trials (Appendix 1—figure 4). We found that it took less than 90 trials for detection at ±32^{°} to become significantly better than average detection over all other presented directions (Appendix 1—figure 4 and pvalues within).
Lastly, for trials where no stimulus was presented, we looked at how long it took participants to start hallucinating predominantly around ±32° as opposed to all other possible directions. This was quantified as a probability ratio p_{rel}:
where N_{bins} is the number of bins (45), each of size 16°. This probability ratio would be equal to one if participants were equally likely to estimate within 8° of ±32° as they were to estimate within other bins. Again, we computed cumulative moving mean at every 30 trials (Appendix 1—figure 5). For participants who did not report seeing dots at any direction within a given number of trials (i.e. zero total hallucinations) this probability ratio was undefined, therefore, those individuals were omitted from significance test at that point. We found that it took less than 210 trials for p_{rel} to become significantly larger than 1 (Appendix 1—figure 5 and pvalues within).
Schizotypy traits and task performance
Appendix 1—figure 6 and Appendix 1—figure 7 show task performance by groups which were formed by splitting the sample on the median RISC and SPQ scores respectively. Appendix 1—figure 8 shows the correlations between RISC and SPQ scores and the corresponding performance measures. There were no significant correlations with any of the measures.
Schizotypy traits and model parameters
Appendix 1—figure 9 shows the robust correlation analysis results between the BAYES model parameter estimates and schizotypy scores. There was no significant correlation with any of the parameters. Further Bayesian correlation analysis provided positive evidence that schizotypy traits had no effect on prior precision (RISC: τ_{b} = 0.012, BF_{01} = 6.90; SPQ: τ_{b} = 0.071, BF_{01} = 3.97).
Individual priors recovered via BAYES model
Appendix 1—figure 10 shows a representative sample of the priors we extracted for a number of individuals, using the ‘BAYES’ model.
Appendix 2
Response bias models
We wanted to account for the possibility that the task behavior might be better explained by simple behavioral strategies. This class of models assumed that on trials when participants were unsure about the presented motion direction they made an estimation based solely on prior expectations, while on the remaining fraction of trials they made unbiased estimates based solely on sensory input.
ADD1
The first model (‘ADD1’) assumed that when participants were unsure about which motion direction they had perceived, they made an estimate that was close to one of the two most frequently presented motion directions. In this model, on each trial, participants make a sensory observation of the stimulus motion direction, θ_{sens}. We parameterize the probability of observing the stimulus to be moving in a direction θ_{sens} by a von Mises (circular normal) distribution centered on the actual stimulus direction and with width determined by 1/κ_{sens}:
On most trials, we assume that participants make a perceptual estimate of the stimulus motion direction (θ_{perc}) that is based entirely on their sensory observation so that θ_{perc} = θ_{sens}. However, on a certain proportion of trials, when participants are uncertain about whether a stimulus was present or not, they resort to their expectations by making a perceptual estimate that is sampled from a learned distribution, p_{exp}(θ). For simplicity, we parameterize this distribution as the sum of two circular normal distributions, each with width determined by 1/κ_{exp}, and centered on motion directions −θ_{ex}_{p} and θ_{exp}, respectively. Finally, we accommodate for the fact that there will be a certain amount of noise associated with moving the estimation bar to indicate which direction the stimulus is moving in as well as allowing for a fraction of trials α, where participants make estimates that are completely random. Thus, the estimation response θ_{est} is related to the perceptual estimate θ_{perc} via the equation:
Bringing all this together, the distribution of estimation responses for a single participant is given by:
where the asterisk denotes a convolution and a(θ) determines the proportion of trials that participants sampled from the expected distribution, p_{exp}(θ). The resulting ‘ADD1’ model has nine free parameters θ_{exp}, κ_{exp}, a(θ) (which can take a different value for each of the five angles: 0°,±16°,±32°,±48°,±64°), κ_{sens} and α.
ADD2
The second model, ‘ADD2’, was just as ‘ADD1’ except that it had slightly more complex strategy for trials when participants were unsure about the stimulus motion direction: instead of sampling from the complete learned probability distribution ranging from −180^{°} to +180^{°} (Equation (11)), they effectively truncated this distribution on a trial by trial basis and sampled from only one part of it, negative (−180° to 0^{°}) or positive (0° to +180^{°}), depending on which side of the distribution the actual stimulus occurred. Incorporating this into the distribution of estimation responses gives:
where asterisk (∗) denotes convolution; a(θ) and b(θ) determine the proportion of trials in which participants sample from either anticlockwise or clockwise distributions p_{expN}(θ) and p_{expP}(θ), respectively.
In addition, we also considered slight variations of the ‘ADD1’ and ‘ADD2’ models, denoted ‘ADD1_m’ and ‘ADD2_m’ respectively. These were identical to ‘ADD1’ and ‘ADD2’ except from setting 1/κ_{exp} to zero; that is, on trials when perceptual estimates were derived only from expectations, they were equal to the mode of the learnt distribution (i.e. no uncertainty).
Nonsymmetric prior models
The stimulus distribution is multimodal and symmetric. Learning such a distribution might be inherently difficult. We reasoned that some individual differences might lie in asymmetries of the acquired priors. Therefore, we explored an alternative parameterization of the acquired priors which allowed them to be asymmetrical. We allowed the two modes in the prior to have different position with respect to 0^{°} and to have different amount of probability associated with each mode. This resulted in:
where π (∈ [0 1]) is a mixing parameter. Using this parameterization we fitted ‘BAYES’ model as described in the main text (thus, we denoted this alternative model as ‘BAYES_π’). The alternative parameterization did not result in a better BIC as compared to ‘BAYES’ model (p=0.378, signed rank test). In addition, we performed parameter recovery to determine how robust ‘BAYES_π’ is and found that recovering the mixing parameter π was not very reliable (r = 0.4), although other parameters retained most of their previous reliability (Appendix 2—figure 1). We thus focused on the simpler model in the current study.
Full models (estimation + detection)
We have built a Bayesian model that incorporates both estimation and detection performance (‘BAYES_full’) in order to fully account for the task behavior. This time, the acquired priors consisted of both the expectations about the direction of stimuli motion (θ) and the expectations about whether stimulus is presented (s = 1) or not (s = 0). It was parameterized as:
where parameter b accounts for a participant’s average expectation that the stimulus will be presented. Thus, we assumed that expectations about motion direction were uniform for when no stimulus was expected. While the expectations about motion direction when the stimulus was expected followed the bimodal probability distribution just as in the previous models.
On each trial, given the presented motion direction (θ_{act}) and the presence of the stimulus (s), participants made sensory measurements p_{sens}(θ_{sens},s_{sens}θ_{act},s). For simplicity, we assumed that the sensory probability of whether the stimulus was present (p_{sens}(s_{sens}θ_{act},s)) was independent of the sensory input about the motion direction (p_{sens}(θ_{sens}θ_{act},s)). We further assumed that s_{sens} was independent of the presented motion direction θ_{act}, as informed by ‘BAYES_var’ model (that allowed the sensory likelihood to vary based on the presented motion direction), which did not produce a better fit. As before, the mean of the motion direction was allowed to fluctuate on trialbytrial basis, such that:
where κ_{sens} is sensory precision. Given the estimate of the mean θ, the sensory input θ_{sens} is represented with the associated uncertainty via:
Putting all this together, the sensory likelihood was expressed as:
where p_{sens}(θ_{sens}θ_{act},s) was parameterized as:
where we assumed that sensory likelihood is uniform when no stimulus is presented. Finally, p_{sens}(s_{sens}s) was parameterized as:
where parameter c is the average probability of detecting dots when they are not presented, and parameter d is the average probability of detecting dots when they are presented. Putting together prior and likelihood, the resulting posterior probability distribution becomes:
With a given posterior participants could have performed detection task at least in two ways. One way is to maximize the posterior (i.e. to always choose the value of s that has higher probability):
Another way is to perform probability matching and choose in accordance to the size of the probabilities:
where η ∈ [0 1] and is drawn for each trial from a uniform distribution. We considered both of these possibilities and implemented a variant of the model for each. Finally, just as in ‘BAYES’ model, the motion direction percept was formed by taking the mean of the posterior:
As previously, we accounted for motor precision and the lapse responses via:
In total, ‘BAYES_full’ model had seven free parameters. To fit the model, in addition to intermediate contrast trials, we also used nostimulus trial data. The rest of the fitting procedure was the same as in the main text: we built a distribution of 1000 posterior estimations for each presented angle and one more distribution of 1000 posterior estimations for no stimulus trials.
We found that ‘BAYES_full’ provided a good fit and captured the main features of both estimation and detection performance (Appendix 2—figure 2). As before, to test how reliable parameters estimated for ‘BAYES_full’ model are, we performed parameter recovery. Just as for ‘BAYES’ parameter recovery described in the main text, we generated 80 sets of parameters and simulated 200 trials of data with ‘BAYES_full’ model for each of them. Then we fitted ‘BAYES_full’ to the simulated data. The results revealed that parameters d and c had very poor recovery (Appendix 2—figure 3). We thus focused on the simpler model in the current study.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source code has been included, which can be used to reproduce the results figures.
References

The computational anatomy of psychosisFrontiers in Psychiatry 4:47.https://doi.org/10.3389/fpsyt.2013.00047

Cognitive deficits in depressionThe British Journal of Psychiatry 178:200–206.https://doi.org/10.1192/bjp.178.3.200

The autismspectrum quotient (AQ): evidence from Asperger syndrome/highfunctioning autism, males and females, scientists and mathematiciansJournal of Autism and Developmental Disorders 31:5–17.https://doi.org/10.1023/A:1005653411471

Motion perception in autism: a "complex" issueJournal of Cognitive Neuroscience 15:218–225.https://doi.org/10.1162/089892903321208150

Alternative Bayesian accounts of autistic perception: comment on Pellicano and BurrTrends in Cognitive Sciences 16:573–574.https://doi.org/10.1016/j.tics.2012.10.005

Multimodel inferenceSociological Methods & Research 33:261–304.https://doi.org/10.1177/0049124104268644

Autistic traits in the general population: a twin studyArchives of General Psychiatry 60:524–530.https://doi.org/10.1001/archpsyc.60.5.524

The lightfromabove prior is intact in autistic childrenJournal of Experimental Child Psychology 161:113–125.https://doi.org/10.1016/j.jecp.2017.04.005

Oblique orientation discrimination thresholds are superior in those with a high level of autistic traitsJournal of Autism and Developmental Disorders 44:2844.https://doi.org/10.1007/s1080301421471

Perceiving is believing: a Bayesian approach to explaining the positive symptoms of schizophreniaNature Reviews Neuroscience 10:48–58.https://doi.org/10.1038/nrn2536

A theory of cortical responsesPhilosophical Transactions of the Royal Society B: Biological Sciences 360:815–836.https://doi.org/10.1098/rstb.2005.1622

The weak coherence account: detailfocused cognitive style in autism spectrum disordersJournal of Autism and Developmental Disorders 36:5–25.https://doi.org/10.1007/s1080300500390

Perceptual learning in autism: overspecificity and possible remediesNature Neuroscience 18:1574.https://doi.org/10.1038/nn.4129

The formation of maintenance of delusions: a Bayesian analysisThe British Journal of Psychiatry 149:51–56.https://doi.org/10.1192/bjp.149.1.51

Robust regression using iteratively reweighted leastsquaresCommunications in Statistics  Theory and Methods 6:813–827.https://doi.org/10.1080/03610927708827533

Robust Estimation of a Location ParameterThe Annals of Mathematical Statistics 35:73–101.https://doi.org/10.1214/aoms/1177703732

Central tendency effects in time interval reproduction in autismScientific Reports 6:28570.https://doi.org/10.1038/srep28570

Bayes FactorsJournal of the American Statistical Association 90:773–795.https://doi.org/10.1080/01621459.1995.10476572

Adults with autism overestimate the volatility of the sensory environmentNature Neuroscience 20:1293–1299.https://doi.org/10.1038/nn.4615

An aberrant precision account of autismFrontiers in Human Neuroscience 8:302.https://doi.org/10.3389/fnhum.2014.00302

Bayesian inference with probabilistic population codesNature Neuroscience 9:1432–1438.https://doi.org/10.1038/nn1790

Children on the autism spectrum update their behaviour in response to a volatile environmentDevelopmental Science 20:e12435.https://doi.org/10.1111/desc.12435

Evidence of a dimensional relationship between schizotypy and schizophrenia: a systematic reviewNeuroscience & Biobehavioral Reviews 37:317–327.https://doi.org/10.1016/j.neubiorev.2013.01.004

What visual illusions teach us about schizophreniaFrontiers in Integrative Neuroscience 8:.https://doi.org/10.3389/fnint.2014.00063

Bayesian approaches to autism: Towards volatility, action, and behaviorPsychological Bulletin 143:521–542.https://doi.org/10.1037/bul0000097

The Importance of Falsification in Computational Cognitive ModelingTrends in Cognitive Sciences 21:425–433.https://doi.org/10.1016/j.tics.2017.03.011

When the world becomes 'too real': a Bayesian explanation of autistic perceptionTrends in Cognitive Sciences 16:504–510.https://doi.org/10.1016/j.tics.2012.08.009

Response to Brock: noise and autismTrends in Cognitive Sciences 16:574–575.https://doi.org/10.1016/j.tics.2012.10.004

Reduced Generalization in Autism: An Alternative to Weak Central CoherenceReduced Generalization in Autism: An Alternative to Weak Central Coherence.

Bayesian models of individual differencesPsychological Science 27:1562–1572.https://doi.org/10.1177/0956797616665351

The SPQ: a scale for the assessment of schizotypal personality based on DSMIIIR criteriaSchizophrenia Bulletin 17:555–564.https://doi.org/10.1093/schbul/17.4.555

The rust inventory of schizotypal cognitions (RISC)Schizophrenia Bulletin 14:317–322.https://doi.org/10.1093/schbul/14.2.317

How much to trust the senses: likelihood learningJournal of Vision 14:13.https://doi.org/10.1167/14.13.13

Delusions and the role of beliefs in perceptual inferenceJournal of Neuroscience 33:13701–13712.https://doi.org/10.1523/JNEUROSCI.177813.2013

Enhanced predictive signalling in schizophreniaHuman Brain Mapping 38:1767–1779.https://doi.org/10.1002/hbm.23480

Perceptual instability in schizophrenia: probing predictive coding accounts of delusions with ambiguous stimuliSchizophrenia Research: Cognition 2:72–77.https://doi.org/10.1016/j.scog.2015.03.005

Learning what to expect (in visual perception)Frontiers in Human Neuroscience 7:668.https://doi.org/10.3389/fnhum.2013.00668

Brief report: suboptimal auditory localization in autism spectrum disorder: support for the bayesian account of sensory symptomsJournal of Autism and Developmental Disorders 46:2539–2547.https://doi.org/10.1007/s1080301627749

The contribution of hypersalience to the "jumping to conclusions" bias associated with delusions in schizophreniaJournal of Psychiatry and Neuroscience 35:7–17.https://doi.org/10.1503/jpn.090025

Synaptic plasticity and dysconnection in schizophreniaBiological Psychiatry 59:929–939.https://doi.org/10.1016/j.biopsych.2005.10.005

Advances in Neural Information Processing Systems1289–1296, Sensory adaptation within a Bayesian framework for perception, Advances in Neural Information Processing Systems.

The WarwickEdinburgh Mental Wellbeing Scale (WEMWBS): development and UK validationHealth and Quality of Life Outcomes 5:.https://doi.org/10.1186/14777525563

A predictive coding perspective on autism spectrum disordersFrontiers in Psychology 4:19.https://doi.org/10.3389/fpsyg.2013.00019

Precise minds in uncertain worlds: predictive coding in autismPsychological Review 121:649–675.https://doi.org/10.1037/a0037665
Decision letter

Klaas Enno StephanReviewing Editor; University of Zurich and ETH Zurich, Switzerland
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Autistic traits, but not schizotypy, predict overweighting of sensory information in Bayesian visual integration" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom, Klaas Enno Stephan (Reviewer #1), is a member of our Board of Reviewing Editors, and the evaluation has been overseen by a Reviewing Editor and Michael Frank as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity; Albert Powers (Reviewer #2); Rick Adams (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
In this carefully conducted study, the authors use a statistical learning paradigm to quantify the contributions of likelihoods and priors to visual perception in healthy participants with different schizotypal and autistic traits. The authors demonstrate an increased precision of sensory information compared to prior beliefs in participants with autistic traits, but do not find alterations in the precision of sensory information or priors in those with schizotypal traits. While all reviewers found this a compelling and important study, there are several issues that would need to be considered in a revision of the paper.
Essential revisions:
1) The study does not examine patients but exploits variability in questionnairebased constructs of personality traits (autistic and schizotypal traits) in the healthy population. These are often thought of as lying on a joint spectrum with the clinical conditions (ASD, schizophrenia). While not unreasonable, all reviewers expressed concern about the implicit message of the paper that its findings will necessarily extend to clinical populations. So far, there is no firm empirical evidence that the spectrum covered by AQ and RISC/SPQ, respectively, directly connects to the spectrum of patients that fulfill DSM/ICD diagnoses. The reader should not be given the impression that the results presented here will necessarily and without doubt extrapolate to the clinical conditions. Whether or not the current results extend to clinical populations will still have to be examined in the future; this is something worth emphasising in the Abstract, the Introduction and the Discussion.
2) The analysis deploys a large number of null hypothesis tests but no correction for multiple testing. It would be helpful if, for each set of tests concerning the same question (e.g. Figure 4, 7), you reported corrected pvalues, perhaps alongside "native" pvalues.
3) The regression plots in Figures 4, 7 show some prominent outliers. It may be worth considering "robust regression" techniques?
4) The negative results for schizotypal traits are described as an absence of effect or relation. This is a little problematic in a frequentist setting (given that one cannot accept the null hypothesis). It would be conceptually cleaner and more interpretable if you computed the Bayes factor for the null hypothesis vs. the alternative hypothesis. A similar issue arises in the subsection “Nostimulus trials and autistic/schizotypal traits” (also discussed in the subsection “Model parameters and autistic/schizotypal traits”) – "suggesting no differences in the acquired prior distributions". As an alternative to a Bayesian analysis, you could state what effect sizes you were powered to reliably detect. This also pertains to the Discussion, and the degree to which this study might contradict other studies.
5) Figure 6D: do you have an explanation why the maxima of the estimated prior probabilities are different from +/ 32 degrees? Is this a reason for concern?
6) Because of the importance of specific symptom dimensions in the literature, it is difficult to reconcile the results at hand with the extant literature without: a) full reporting of SPQ and RISC scores in that group; and b) specific data on the dimensions of schizotypy most present in this sample. A division of scores based upon a factor analysis or relevance to hallucinations would be particularly helpful. Similarly, a breakdown of SPQ based upon symptom dimensions (especially those most relevant to perception) could be enlightening.
7) Similarly, while it is true that the authors focus mostly on ASD in this paper, it would be nice to see data (broken down similar to Figure 3) for the schizotypal group. This would allow the reader to compare performance across groups, which is critical for interpretation.
8) The authors' idea of testing learning rate in autism is a good one, and the authors are correct in stating that their findings support the hypothesis that there may be an increased learning rate in ASD. However, this is not consistent with one work the authors cite (Lawson et al., 2017). Lawson and colleagues find increased volatility estimates in ASD, driving decreased learning. The authors should attempt to reconcile this finding with theirs and perhaps temper their stated expectations for future investigations into learning rate in ASD.
9) An important contrast with the Teufel et al. study is that they used the CAPS and PDI to assess schizotypal hallucinatory and delusional propensity respectively. They found that the CAPS explained all the relation to increased prior 'precision' – i.e. there was no correlation with PDI once CAPS score was partialed out. This implies that greater prior 'precision' is not really related to schizotypy per se but to the subconstruct of aberrant perception. In the 74 item SPQ, only 9 questions refer to perceptual experiences (one could debate whether ideas of reference are also relevant as they may begin with vivid percepts). I'm not familiar with the RISC but I imagine not many of its questions are perceptual. So it seems doubtful that you can really conclude much about your findings in relation to the Teufel et al. study, given such different measures were used. Please mention this in the Discussion. It would be interesting to know whether you detect effects if you only use RISC/SPQ scores relating to perceptual phenomena.
https://doi.org/10.7554/eLife.34115.031Author response
Essential revisions:
1) The study does not examine patients but exploits variability in questionnairebased constructs of personality traits (autistic and schizotypal traits) in the healthy population. These are often thought of as lying on a joint spectrum with the clinical conditions (ASD, schizophrenia). While not unreasonable, all reviewers expressed concern about the implicit message of the paper that its findings will necessarily extend to clinical populations. So far, there is no firm empirical evidence that the spectrum covered by AQ and RISC/SPQ, respectively, directly connects to the spectrum of patients that fulfill DSM/ICD diagnoses. The reader should not be given the impression that the results presented here will necessarily and without doubt extrapolate to the clinical conditions. Whether or not the current results extend to clinical populations will still have to be examined in the future; this is something worth emphasising in the Abstract, the Introduction and the Discussion.
We agree with the reviewers that the question of whether or not the current results extend to clinical populations is an open one. We have modified the Introduction, the Discussion and the Conclusion to make this point clearer.
2) The analysis deploys a large number of null hypothesis tests but no correction for multiple testing. It would be helpful if, for each set of tests concerning the same question (e.g. Figure 4, 7), you reported corrected pvalues, perhaps alongside "native" pvalues.
We now include Bonferroni correction for multiple comparisons and report it alongside the native values where appropriate. The native values now correspond to robust regression rather than to rank correlations (see the response to the essential revision 3). The Bonferroni correction itself was based on the number of independent hypotheses that we tested; that is, we were interested whether ASD and schizotypy had any effect on acquired priors or sensory likelihoods, which resulted in 4 different hypotheses. Even though the number of the null hypothesis significance tests that we performed was higher, the tests all concerning a particular hypothesis were highly dependent (e.g., behavioral biases/variability and priors/likelihoods as recovered via modelling are directly causally related; RISC and SPQ both aim to measure schizotypy construct). Including each individual test as independent would have resulted in overcorrection of pvalues, as Bonferroni rests on the assumption that the tests are independent. We include these details in the Statistical tests section in the Materials and methods. The main result of our study, the negative relationship between AQ and sensory uncertainty, remained significant after applying the correction. Furthermore, decreased behavioral variability and decreased number of hallucinations in ASD also remained significant. The only relationship that did not survive the correction was the one between AQ and estimation bias.
3) The regression plots in Figures 4, 7 show some prominent outliers. It may be worth considering "robust regression" techniques?
We have rerun the analysis using robust regression instead of rank correlations. This was done using robustfit function in Matlab, which downweighs the influence of outliers via iteratively reweighted least squares (IRLS). For the loss function we used Huber function with a tuning constant of 1.345, which corresponds to 95% estimator efficiency as compared to ordinary least squares. We include these details in the Statistical tests section in the Materials and methods. The results are for the most part consistent with the previous Spearman’s rank correlation analysis. The only difference that is worth noting is the relationship between estimation bias and AQ, which went from significant to trending (i.e., from rho = 0.228, p = 0.039 to r = 0.175, p = 0.053).
4) The negative results for schizotypal traits are described as an absence of effect or relation. This is a little problematic in a frequentist setting (given that one cannot accept the null hypothesis). It would be conceptually cleaner and more interpretable if you computed the Bayes factor for the null hypothesis vs. the alternative hypothesis. A similar issue arises in the subsection “Nostimulus trials and autistic/schizotypal traits” (also discussed in the subsection “Model parameters and autistic/schizotypal traits”) – "suggesting no differences in the acquired prior distributions". As an alternative to a Bayesian analysis, you could state what effect sizes you were powered to reliably detect. This also pertains to the Discussion, and the degree to which this study might contradict other studies.
We thank the reviewers for bringing this up. We now report Bayes Factors for the null hypothesis (BF01) that schizotypy and autistic traits have no effect on the acquired priors. The Bayesian analysis was done using JASP (Version 0.8.6). BF01 were computed for the nonparametric Kendall’s Taub correlation coefficient, which maintains a degree of robustness (as compared to Pearson’s correlation – the only other alternative available in JASP). We provide all the details to the reader in the Statistical tests section in the Materials and methods.
5) Figure 6D: do you have an explanation why the maxima of the estimated prior probabilities are different from +/ 32 degrees? Is this a reason for concern?
This phenomenon is indeed intriguing and, aside from the interindividual variability possibly not completely averaging out, we do not have a good explanation for it. One possibility is that this may be due to how the circular space is sampled around the most frequent directions: there is only a narrow 32 degrees range (from 0 to 32 in between the two frequent directions) where estimates would shift the average “inwards”, in the opposite direction than we currently observe. Whereas estimates following in the remaining – much larger – 32 to 180 range would be leading to “outwards” biases, like what we see now. If the interindividual variability is sufficiently large, such asymmetry could result in less individuals’ estimations within the 0 to 32 range then within the remaining 32 to 180 range, shifting the group mean accordingly. This phenomenon has been consistently present in previous studies using this task, including the original study by Chalk et al. (2010). It is also worth pointing out that the modes being shifted away from +/ 32 are consistent with the estimation biases in the raw behavioral data, which are also nonzero at +/ 32, therefore it is not a parameter estimation issue.
6) Because of the importance of specific symptom dimensions in the literature, it is difficult to reconcile the results at hand with the extant literature without: a) full reporting of SPQ and RISC scores in that group; and b) specific data on the dimensions of schizotypy most present in this sample. A division of scores based upon a factor analysis or relevance to hallucinations would be particularly helpful. Similarly, a breakdown of SPQ based upon symptom dimensions (especially those most relevant to perception) could be enlightening.
The distribution (mean, standard deviation and range) of SPQ and RISC scores in our sample is reported in ‘Task performance and autistic/schizotypy traits’ section in Results. We agree that a breakdown of SPQ based upon symptom dimensions could be interesting. Unfortunately, we are not able to perform this posthoc analysis at this time as, due to an unfortunate clerical error, the subscores of each dimension of SPQ have not been kept on record (only the total score). RISC however does not include any dimensions that would be considered primarily perceptual (it consists of dimensions of suspiciousness, magical ideation, ritualistic thinking, thought isolation and delusional ideation). Whether inconsistencies between our findings and those reported by previous studies might be underlied by differences in the aberrant perception dimension in the studied groups will be studied in future work.
7) Similarly, while it is true that the authors focus mostly on ASD in this paper, it would be nice to see data (broken down similar to Figure 3) for the schizotypal group. This would allow the reader to compare performance across groups, which is critical for interpretation.
We have now included such figures for both RISC and SPQ scores in the Appendix 1 (Figures 6 and 7).
8) The authors' idea of testing learning rate in autism is a good one, and the authors are correct in stating that their findings support the hypothesis that there may be an increased learning rate in ASD. However, this is not consistent with one work the authors cite (Lawson et al., 2017). Lawson and colleagues find increased volatility estimates in ASD, driving decreased learning. The authors should attempt to reconcile this finding with theirs and perhaps temper their stated expectations for future investigations into learning rate in ASD.
We thank the reviewers for pointing this out. We have included a short discussion about the apparent inconsistency between our findings and those of Lawson et al. (2017). With regards to this, the most relevant difference between our paradigm and theirs, is that Lawson et al. (2017) introduce several (hierarchical) levels of uncertainty and their main results concern an effect that we don’t investigate (our environment being static): the change in the learning rates between stable and volatile conditions, rather than the learning rate themselves. Interestingly, for the average learning rates (which we speculate to be increased in ASD), they find no effect.
9) An important contrast with the Teufel et al. study is that they used the CAPS and PDI to assess schizotypal hallucinatory and delusional propensity respectively. They found that the CAPS explained all the relation to increased prior 'precision' – i.e. there was no correlation with PDI once CAPS score was partialed out. This implies that greater prior 'precision' is not really related to schizotypy per se but to the subconstruct of aberrant perception. In the 74 item SPQ, only 9 questions refer to perceptual experiences (one could debate whether ideas of reference are also relevant as they may begin with vivid percepts). I'm not familiar with the RISC but I imagine not many of its questions are perceptual. So it seems doubtful that you can really conclude much about your findings in relation to the Teufel et al. study, given such different measures were used. Please mention this in the Discussion. It would be interesting to know whether you detect effects if you only use RISC/SPQ scores relating to perceptual phenomena.
We thank the reviewers for sharing this clarification. Indeed, very few items in RISC could be construed to be related to aberrant perception, as the questionnaire does not explicitly include such dimension (it consists of dimensions of suspiciousness, magical ideation, ritualistic thinking, thought isolation and delusional ideation). Thus, the subconstruct of aberrant perception could indeed potentially explain the differences in findings. We have revised the Discussion accordingly.
https://doi.org/10.7554/eLife.34115.032Article and author information
Author details
Funding
Brain and Behavior Research Foundation (NARSAD Young investigator grant 19271)
 Peggy Seriès
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Gizem Aras for assisting in data collection, and Katie Richards for assisting with participants’ recruitment. This work was done in part while PS was visiting the Simons Institute for the Theory of Computing.
Ethics
Human subjects: All participants gave informed written consent and received monetary compensation for participation. The study was approved by the University of Edinburgh School of Informatics Ethics Panel.
Reviewing Editor
 Klaas Enno Stephan, University of Zurich and ETH Zurich, Switzerland
Publication history
 Received: December 8, 2017
 Accepted: May 6, 2018
 Accepted Manuscript published: May 14, 2018 (version 1)
 Version of Record published: May 23, 2018 (version 2)
Copyright
© 2018, Karvelis et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 3,773
 Page views

 481
 Downloads

 18
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.