Autistic traits, but not schizotypy, predict increased weighting of sensory information in Bayesian visual integration

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Data availability
References
Article and author information
Metrics

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.001

Introduction

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 low-level prior expectations might lead to forming compensatory strong and idiosyncratic high-level 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.

Figure 1

Download asset Open asset

Alternative hypotheses for ASD impairments within the Bayesian inference framework.

In Bayesian terms, the percept can be described as a posterior distribution, which is a combination of sensory information (likelihood) and prior expectations (prior). Two contrasting hypotheses have been proposed to underlie behavioral differences in ASD: enhanced sensory precision, that is, smaller σ_sens (left) vs. attenuated priors, that is, larger σ_exp (right). Both hypotheses predict a reduced influence (bias) of the prior on the location of the posterior distribution (posterior mean). However, these alternatives differ in their predictions for perceptual variability, which is determined by the posterior width: the enhanced sensory precision hypothesis should lead to reduced variability while the attenuated prior hypothesis should lead to increased variability. By measuring both bias and variability, our experimental paradigm can distinguish between these two hypotheses.

https://doi.org/10.7554/eLife.34115.002

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 slow-speed 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 light-from-above 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 decision-making task (Manning et al., 2017) but a follow-up 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 two-tone 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.

Figure 2

Download asset Open asset

The moving dots task.

(a) Sequence of events on a single trial. First, a fixation point is presented. Next, a field of coherently moving dots is presented along with an estimation bar (extending from the fixation point) which participants are required to move to indicate perceived motion direction. Lastly, in a two-alternative forced choice, participants are asked to report whether they saw the dots during the estimation part (detection task). (b) The probability of different motion directions being presented: directions at ±32° are presented more often than other directions. Motion direction is plotted relative to a central reference angle (at 0°), which was randomly set for each participant.

https://doi.org/10.7554/eLife.34115.003

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

Download asset Open asset

Average group performance on low-contrast trials (**a–d**) and on trials with no stimulus (e).

(a) Mean estimation bias, (b) standard deviation of estimations, (c) estimation reaction time and (d) fraction of trials in which the stimulus was detected. (f) Probability distribution of estimation responses on trials without stimulus. The solid line denotes the estimation responses when participants reported detecting a stimulus (hallucinations). The dash-dot line denotes estimation distributions when participants correctly reported not detecting a stimulus. (e) Distribution of hallucinations for high and low AQ groups (median split). The vertical dashed lines correspond to the two most frequently presented motion directions (±32°). Error bars and shaded areas represent within-subject standard error.

https://doi.org/10.7554/eLife.34115.004

Figure 3—source data 1 This zip archive contains .csv files with all of the data that was used to produce plots in Figure 3. EstimationBias.csv contains estimation biases at each of the nine presented angles. EstimationVariability.csv contains standard deviation of estimations at each of the nine presented angles. NostimDetected.csv and NostimUndetected.csv contain estimation responses when stimulus was detected and not detected, respectively, on no-stimulus trials. Traits.csv contains AQ scores of each individual (column 3) as well as all other traits. SourceData_Readme.txt contains more detailed description of each data file. The plots can be reproduced from MATLAB script master.m which is available in the provided Source code 1. SourceCode_Readme.txt contains more detailed description of the source code.: https://doi.org/10.7554/eLife.34115.005
Download elife-34115-fig3-data1-v2.zip

No-stimulus 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:

p_{r e l} = p (θ_{e s t} = \pm 32 (\pm 8)^{\circ}) \cdot N_{b i n s},

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, dash-dot 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

Download asset Open asset

Correlations between AQ scores and task performance on low contrast trials (**a, b**) and when no stimulus is presented (c).

(a) Mean absolute bias (r = −0.175, p=0.053), (b) mean standard deviation (i.e. variability) of estimations (r = −0.327, p<0.001), and (c) the total number of hallucinations (r = −0.238, p=0.010). The blue lines are robust regression slopes.

https://doi.org/10.7554/eLife.34115.006

Figure 4—source data 1 This zip archive contains .csv files with all of the data that was used to produce plots in Figure 4. EstimationBias.csv contains estimation biases at each of the nine presented angles. EstimationVariability.csv contains standard deviation of estimations at each of the nine presented angles. NostimDetected.csv contains the number of hallucinations at different directions. Traits.csv contains AQ scores of each individual (column 3) as well as all other traits. SourceData_Readme.txt contains more detailed description of each data file. The plots were produced with MATLAB script analyze_data.m which is available in the provided Source code 1. SourceCode_Readme.txt contains more detailed description of the source code.: https://doi.org/10.7554/eLife.34115.007
Download elife-34115-fig4-data1-v2.zip

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.

No-stimulus 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 no-stimulus 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 different-sized 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₀₁ = 7.24; RISC - τ_b = -0.050, BF₀₁ = 3.73; SPQ - τ_b = 0.101, BF₀₁ = 8.72; within 16° of ±32°: AQ - τ_b = -0.068, BF₀₁ = 2.86; RISC - τ_b = -0.129, BF₀₁ = 0.84; SPQ - τ_b = 0.018, BF₀₁ = 5.45; within 24° of ±32°: AQ - τ_b = 0.057, BF₀₁ = 11.67; RISC - τ_b = -0.078, BF₀₁ = 2.40; SPQ - τ_b = 0.006, BF₀₁ = 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).

Figure 5

Download asset Open asset

Bayesian model of estimation response for a single trial.

The actual motion direction (θ_act) is corrupted by sensory uncertainty (σ_sens), and then combined with prior expectations (mean θ_exp and uncertainty σ_exp) to form a posterior distribution. The perceptual estimate (θ_perc) is defined as the mean of the posterior distribution. Finally, motor precision ( $1 / σ_{m}^{2}$ ) and a probability of random response (α) are incorporated to generate the response (θ_est). This results in four free model parameters: σ_sens, σ_exp, θ_exp and α. The motor precision is estimated from high contrast trials and is used as a fixed parameter.

https://doi.org/10.7554/eLife.34115.008

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 p-values within Figure 6a).

Figure 6

Download asset Open asset

Modelling results.

(a) Model comparison for all participants using Bayesian Information Criterion (BIC). y-axis measures the relative difference between BIC of each model (as indicated on the x-axis) and BIC of BAYES model. Values greater than zero on the y-axis indicate that the BAYES model provided a better fit. Each dot represents a participant. Red horizontal lines denote median values; blue horizontal lines denote 25th and 75th percentiles. p-values above the plot indicate whether the median of the difference was significantly different from zero for each model (signed rank test). Panels (a) and (c) present task performance at different motion directions as predicted by BAYES model: (b) estimation bias, (c) standard deviation of estimations. Error bars represent within-subject standard error. (d) Population averaged prior as recovered via BAYES model. The vertical dashed lines correspond to the two most frequently presented motion directions (±32°).

https://doi.org/10.7554/eLife.34115.009

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₀₁ = 6.99).

Figure 7

Download asset Open asset

Correlations between AQ scores and BAYES model parameters.

(a) θ_exp - mean of the prior expectations (r = 0.031, p=0.820), (b) σ_exp - uncertainty of the prior distribution (r = 0.018, p=0.962), (c) σ_sens - uncertainty in the sensory likelihood (r = −0.185, p=0.011) and (d) α - fraction of random estimations (r = −0.135, p=0.238). The blue lines are robust regression slopes.

https://doi.org/10.7554/eLife.34115.010

Figure 7—source data 1 This zip archive contains .csv files with all of the data that was used to produce plots in Figure 7. BayesEstimatedParams.csv contains BAYES model parameter estimates. Traits.csv contains AQ scores of each individual (column 3) as well as all other traits. SourceData_Readme.txt contains more detailed description of each data file. The plots were produced with MATLAB script analyze_params.m which is available in the provided Source code 1. The SourceCode_Readme.txt contains more detailed description of the source code.: https://doi.org/10.7554/eLife.34115.011
Download elife-34115-fig7-data1-v2.zip

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₀₁ = 6.90; SPQ: τ_b = 0.071, BF₀₁ = 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).

Figure 8

Download asset Open asset

Comparison of actual (x-axis) vs. recovered (y-axis) parameters using the BAYES’ model.

(a) θ_exp - mean of the prior expectations (r = 0.90), (b) σ_exp - uncertainty of the prior distribution (r = 0.92), (c) σ_sens - uncertainty in the sensory likelihood (r = 0.95), (d) α - fraction of random estimations (r = 0.98). The dashed diagonal line is a reference line indicating perfect parameter recovery.

https://doi.org/10.7554/eLife.34115.012

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 two-tone versions of images before and after seeing the actual images in full color and had to report whether the presented two-tone image contains a face. The low-level 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 low-level (perceptually-induced - stabilizing) priors but stronger high-level (cognitively-induced) 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 first-order (luminance-defined) 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 top-down 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 first-order stimulus in autism, while for second-order (texture-defined) 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 direction-of-gaze 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 context-sensitive in autism. This is in line with the view that generalization is impaired in autism (Plaisted, 2015). Furthermore, such over-specificity 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 protocol

ASD was assessed using 50-item version Autism Spectrum Quotient (AQ) (Baron-Cohen et al., 2001), which is commonly used for assessing milder variants of autistic-like 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 sub-group of 41 participants also completed Schizotypal Personality Questionnaire (SPQ) (Raine, 1991). Finally, all participants were also asked to complete the Warwick-Edinburgh Mental Well-being Scale (WEMWBS) (Tennant et al., 2007) in order to control for potential depression-induced 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² 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 protocol

The 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²) 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², 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 protocol

The stimuli were presented at four different levels of contrast: 0 contrast (no-stimulus 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ía-Pé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² 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 protocol

Responses 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:

(1 - α) \cdot V (θ | μ, κ) + α,

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 no-stimulus 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 no-stimulus trials. We interpreted these false alarms as a simple form of perceptual hallucination.

Modelling

Bayesian models

Request a detailed protocol

Bayesian 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

p_{s e n s} (θ_{s e n s} | θ_{a c t}) = V (θ_{t}, κ_{s e n s}) .

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:

p_{\exp} (θ) = 0.5 [V (- θ_{\exp}, κ_{\exp}) + V (θ_{\exp}, κ_{\exp})]

Combining these via Bayes’ rule gives a posterior probability that the stimulus is moving in a direction θ:

p_{p o s t} (θ | θ_{s e n s}) \propto p_{\exp} (θ) \cdot p_{s e n s} (θ_{s e n s} | θ)

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:

p (θ_{e s t} | θ_{p e r c}) = (1 - α) \cdot V (θ_{p e r c}, κ_{m}) + α,

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 protocol

We 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 protocol

We 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 log-likelihood.

Model comparison

Request a detailed protocol

To compare the model fits we used Bayesian Information Criterion (BIC), which approximates the log of model evidence (Burnham and Anderson, 2004):

- 2 \cdot l o g (P (D | M)) \approx B I C = - 2 \cdot l o g (P (D | M, \hat{Θ})) + k \cdot l o g (n),

where M is model, D is observed data and P (D|M, $\hat{Θ}$ ) is the likelihood of generating the experimental data given the most likely set of parameters, $\hat{Θ}$ ; 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 protocol

To 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 protocol

Due 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₀₁). This was done using JASP (Team, 2017) (Version 0.8.6). Due to the presence of outliers, this analysis was carried out using the non-parametric Kendall’s Tau-b correlation coefficient.

Source code and data

Request a detailed protocol

The 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.

Appendix 1—figure 1

Download asset Open asset

Task performance at the highest contrast level and exclusion Criteria.

Left panel: fraction of detected high contrast trials - quantified as the fraction of trials in which participants both validated their choice with a click within 3000 ms in the estimation part and reported seeing dots (clicked ‘DOTS’) in the detection part. Right panel: root mean square error of estimations on high contrast trials. The dashed lines represent minimum performance criteria (more than 80% detection and less than 30^° RMS error of estimations). Excluded participants are denoted by cross markers.

https://doi.org/10.7554/eLife.34115.016

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² for the other sub-group of 41 (right panel) it was set to 5.18 cd/m². 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² and 4/1 being 0.22 ± 0.04 cd/m² above the 1.16 cd/m² background luminance; and on average 2/1 being 0.42 ± 0.05 cd/m² and 4/1 being 0.46 ± 0.05 cd/m² above the 5.18 cd/m² background luminance. Thus, the two intermediate contrasts were combined for all further data analysis.

Appendix 1—figure 2

Download asset Open asset

Population averaged stimulus contrast relative to the background contrast for the 2/1 (red) and 4/1 (black) staircased contrast levels.

Standard deviation is denoted by shaded areas with corresponding colors. The vertical dashed line marks 170 trials. Left panel: 44 participants (remaining after exclusion) that performed the task with the background luminance set to 1.16 cd/m². Right panel: 39 participants (remaining after exclusion) that performed the task with the background luminance set to 5.18 cd/m².

https://doi.org/10.7554/eLife.34115.017

Combining the different background luminance levels

To compare the two sub-groups that performed the task at different background luminance levels, we performed Wilcoxon two-tailed 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 p-values within).

Appendix 1—figure 3

Download asset Open asset

Cumulative moving average of ratio of estimation reaction times at ±32^° vs average reaction times at all other directions.

Red bars indicate median values and blue bars indicate 25th and 75th percentiles. p-values indicate whether RTs at ±32^° are significantly shorter than average RTs over all other directions (one-tailed Wilcoxon signed rank test).

https://doi.org/10.7554/eLife.34115.018

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 p-values within).

Appendix 1—figure 4

Download asset Open asset

Cumulative moving average of ratio of fraction of detected stimuli at ±32^° vs average fraction detected at all other directions.

Red bars indicate median values and blue bars indicate 25th and 75th percentiles. p-values indicate whether fraction detected at ±32^° are significantly larger than average fraction detected over all other directions (one-tailed Wilcoxon signed rank test).

https://doi.org/10.7554/eLife.34115.019

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:

p_{r e l} = p (θ_{e s t} = \pm 32 (\pm 8)^{\circ}) \cdot N_{b i n s},

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 p-values within).

Appendix 1—figure 5

Download asset Open asset

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.

Appendix 1—figure 6

Download asset Open asset

Average group performance on low-contrast trials (**a–d**) and on trials with no stimulus (e) by groups split by median RISC score.

(a) Mean estimation bias, (b) standard deviation of estimations, (c) estimation reaction time and (d) fraction of trials in which the stimulus was detected. (e) Distribution of hallucinations. The vertical dashed lines correspond to the two most frequently presented motion directions (±32°). Error bars and shaded areas represent within-subject standard error.

https://doi.org/10.7554/eLife.34115.021

Appendix 1—figure 7

Download asset Open asset

Average group performance on low-contrast trials (**a–d**) and on trials with no stimulus (e) by groups split by median SPQ score.

(a) Mean estimation bias, (b) standard deviation of estimations, (c) estimation reaction time and (d) fraction of trials in which the stimulus was detected. (e) Distribution of hallucinations. The vertical dashed lines correspond to the two most frequently presented motion directions (±32°). Error bars and shaded areas represent within-subject standard error.

https://doi.org/10.7554/eLife.34115.022

Appendix 1—figure 8

Download asset Open asset

Correlations between personality traits, RISC (top row) and SPQ (bottom row) and task performance.

There were no significant correlations with any of the measures: mean absolute bias (left column), mean estimation variability (middle column) and total number of hallucinations (right column). Robust correlation coefficients and p-values are indicated above each plot. The blue lines denote robust regression.

https://doi.org/10.7554/eLife.34115.023

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₀₁ = 6.90; SPQ: τ_b = 0.071, BF₀₁ = 3.97).

Appendix 1—figure 9

Download asset Open asset

Correlations with the BAYES model parameter values and schizotypy traits (as measured by both RISC and SPQ).

First column: *θ_exp* - mean of the prior expectations, second column: *σ_exp* - uncertainty of the prior distribution, third column: *σ_sens* - uncertainty in the sensory likelihood and fourth column: α - fraction of random estimations. Robust correlation coefficients and p-values are indicated above each plot. The blue lines denote robust regression.

https://doi.org/10.7554/eLife.34115.024

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 1—figure 10

Download asset Open asset

A representative sample of prior expectations for each individual as reconstructed via ‘BAYES’ model.

The dashed lines correspond to the two most frequently presented motion directions (±32^°).

https://doi.org/10.7554/eLife.34115.025

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:

p_{s e n s} (θ_{s e n s} | θ_{a c t}) = V (θ_{a c t}, κ_{s e n s})

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:

p (θ_{e s t} | θ_{p e r c}) = (1 - α) * V (θ_{p e r c}, κ_{m}) + α .

Bringing all this together, the distribution of estimation responses for a single participant is given by:

p (θ_{e s t} | θ_{a c t}) = (1 - α) [(1 - a (θ)) p_{s e n s} (θ_{s e n s} = θ_{e s t} | θ_{a c t}) + a (θ) p_{e x p} (θ)] * V (0, κ_{m}) + α .

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:

\begin{array}{ll} p (θ_{e s t} | θ_{a c t}) = & (1 - α) [(1 - a (θ) - b (θ)) p_{s e n s} (θ_{s e n s} = θ_{e s t} | θ_{a c t}) \\ + a (θ) p_{e x p N} (θ) + b (θ) p_{e x p P} (θ)] * V (0, κ_{m}) + α . \end{array}

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).

Non-symmetric 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:

p_{e x p} (θ) = (1 - π) \cdot V (θ_{p}, κ_{e x p}) + π \cdot V (θ_{n}, κ_{e x p})

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.

Appendix 2—figure 1

Download asset Open asset

Comparison of actual and recovered parameters via ‘BAYES_π’ model.

*θ_p* and *θ_n* - positive and negative modes of the bimodal distribution of prior expectations, *σ_exp* - uncertainty of the prior distribution, *σ_sens* uncertainty in the sensory likelihood, α - fraction of random estimations, π - mixing parameter responsible for the degree of bimodality. Actual parameters are scattered along x-axis and recovered parameters are scattered along y-axis. The dashed diagonal line is a reference line indicating perfect parameter recovery. Pearson’s correlation coefficients are indicated above each plot.

https://doi.org/10.7554/eLife.34115.027

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:

p_{e x p} (θ, s) = {\begin{matrix} (1 - b) . \frac{1}{2 π}, & i f s = 0 \\ b . \frac{1}{2} [V (- θ_{e x p}, κ_{e x p}) + V (θ_{e x p}, κ_{e x p})], & i f s = 1 \end{matrix}

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 trial-by-trial basis, such that:

p (θ | θ_{a c t}) = V (θ_{a c t}, κ_{s e n s}),

where κ_sens is sensory precision. Given the estimate of the mean θ, the sensory input θ_sens is represented with the associated uncertainty via:

p_{s e n s} (θ_{s e n s} | θ) = V (θ, κ_{s e n s}) .

Putting all this together, the sensory likelihood was expressed as:

p_{s e n s} (θ_{s e n s}, s_{s e n s} | θ, s) = p_{s e n s} (θ_{s e n s} | θ, s) p (s_{s e n s} | s)

where p_sens(θ_sens|θ_act,s) was parameterized as:

p_{s e c s} (θ_{s e n s} | θ_{a c t}, s) = {\begin{matrix} \frac{1}{2 π}, & i f s = 0 \\ V (θ, κ_{s e n s}), & i f s = 1 \end{matrix}

where we assumed that sensory likelihood is uniform when no stimulus is presented. Finally, p_sens(s_sens|s) was parameterized as:

p_{s e n s} (s_{s e n s} = {0, 1} | s) = {\begin{matrix} {1 - c, c}, & i f s = 0 \\ {1 - d, d}, & i f s = 1 \end{matrix}

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:

p_{p o s t} (θ, s | θ_{s e n s}, s_{s e n s}) α p_{s e n s} (θ_{s e n s} | θ, s) \cdot p_{s e n s} (s_{s e n s} | s) \cdot p_{e x p} (θ, s),

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):

s_{p e r c} = a r g m a x_{s} [p_{p o s t} (s | θ_{s e n s}, s_{s e n s})]

Another way is to perform probability matching and choose in accordance to the size of the probabilities:

s_{p e r c} = {\begin{matrix} 0, & i f p_{p o s t} (s = 0 | θ_{s e n s}, s_{s e n s}) > η \\ 1, & i f p_{p o s t} (s = 0 | θ_{s e n s}, s_{s e n s}) < η \end{matrix}

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:

θ_{p e r c} = \int θ . p_{p o s t} (θ | θ_{s e n s}, s_{s e n s}) d θ = \frac{1}{Z} \int θ . \sum_{s} p_{e x p} (θ) . p_{s e n s} (θ_{s e n s} | θ, s) . p_{s e n s} (s_{s e n s} | s) d θ,

As previously, we accounted for motor precision and the lapse responses via:

p (θ_{e s t} | θ_{p e r c}) = (1 - α) \cdot V (θ_{p e r c}, κ_{m o t o r}) + α \cdot p_{e x p} (θ) * V (0, κ_{m o t o r}) .

In total, ‘BAYES_full’ model had seven free parameters. To fit the model, in addition to intermediate contrast trials, we also used no-stimulus 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.

Appendix 2—figure 2

Download asset Open asset

Task performance as predicted by the BAYES_full model.

Left panel: mean estimation bias at different motion directions. Middle panel: standard deviation of estimations at different motion directions. Right panel: fraction of detected stimuli at different motion directions. The dashed lines correspond to the two most frequently presented motion directions (±32^°). Error bars represent within-subject standard error.

https://doi.org/10.7554/eLife.34115.028

Appendix 2—figure 3

Download asset Open asset

Comparison of actual and recovered parameters via ‘BAYES_full’ model.

*θ_exp* - the mean of prior expectations of motion direction, *σ_exp* - uncertainty of the prior expectations of motion direction, *σ_sens* - uncertainty in the sensory likelihood, α - fraction of random estimations, b - prior expectation for dots being presented, c likelihood of detecting the dots when they are not presented, d - likelihood of detecting the dots when they are presented. Actual parameters are scattered along x-axis and recovered parameters are scattered along y-axis. The dashed diagonal line is a reference line indicating perfect parameter recovery.

https://doi.org/10.7554/eLife.34115.029

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

1. Adams RA
2. Stephan KE
3. Brown HR
4. Frith CD
5. Friston KJ
(2013) The computational anatomy of psychosis
Frontiers in Psychiatry 4:47.

https://doi.org/10.3389/fpsyt.2013.00047
- PubMed
- Google Scholar
(2001) Cognitive deficits in depression
The British Journal of Psychiatry 178:200–206.

https://doi.org/10.1192/bjp.178.3.200
- Google Scholar
(2001) The autism-spectrum quotient (AQ): evidence from Asperger syndrome/high-functioning autism, males and females, scientists and mathematicians
Journal of Autism and Developmental Disorders 31:5–17.

https://doi.org/10.1023/A:1005653411471
- PubMed
- Google Scholar
(2003) Motion perception in autism: a "complex" issue
Journal of Cognitive Neuroscience 15:218–225.

https://doi.org/10.1162/089892903321208150
- PubMed
- Google Scholar
(2005) Enhanced and diminished visuo-spatial information processing in autism depends on stimulus complexity
Brain 128:2430–2441.

https://doi.org/10.1093/brain/awh561
- PubMed
- Google Scholar
1. Brainard DH
(1997) The Psychophysics Toolbox
Spatial Vision 10:433–436.

https://doi.org/10.1163/156856897X00357
- PubMed
- Google Scholar
1. Brock J
(2012) Alternative Bayesian accounts of autistic perception: comment on Pellicano and Burr
Trends in Cognitive Sciences 16:573–574.

https://doi.org/10.1016/j.tics.2012.10.005
- PubMed
- Google Scholar
1. Burnham KP
2. Anderson DR
(2004) Multimodel inference
Sociological Methods & Research 33:261–304.

https://doi.org/10.1177/0049124104268644
- Google Scholar
(2010) Rapidly learned stimulus expectations alter perception of motion
Journal of Vision 10:2.

https://doi.org/10.1167/10.8.2
- PubMed
- Google Scholar
1. Constantino JN
2. Todd RD
(2003) Autistic traits in the general population: a twin study
Archives of General Psychiatry 60:524–530.

https://doi.org/10.1001/archpsyc.60.5.524
- PubMed
- Google Scholar
(2009) From drugs to deprivation: a Bayesian framework for understanding models of psychosis
Psychopharmacology 206:515–530.

https://doi.org/10.1007/s00213-009-1561-0
- PubMed
- Google Scholar
1. Croydon A
2. Karaminis T
3. Neil L
4. Burr D
5. Pellicano E
(2017) The light-from-above prior is intact in autistic children
Journal of Experimental Child Psychology 161:113–125.

https://doi.org/10.1016/j.jecp.2017.04.005
- PubMed
- Google Scholar
(2014) Oblique orientation discrimination thresholds are superior in those with a high level of autistic traits
Journal of Autism and Developmental Disorders 44:2844.

https://doi.org/10.1007/s10803-014-2147-1
- PubMed
- Google Scholar
1. Fletcher PC
2. Frith CD
(2009) Perceiving is believing: a Bayesian approach to explaining the positive symptoms of schizophrenia
Nature Reviews Neuroscience 10:48–58.

https://doi.org/10.1038/nrn2536
- PubMed
- Google Scholar
1. Friston K
(2005) A theory of cortical responses
Philosophical Transactions of the Royal Society B: Biological Sciences 360:815–836.

https://doi.org/10.1098/rstb.2005.1622
- Google Scholar
1. García-Pérez MA
(1998) Forced-choice staircases with fixed step sizes: asymptotic and small-sample properties
Vision Research 38:1861–1881.

https://doi.org/10.1016/S0042-6989(97)00340-4
- PubMed
- Google Scholar
1. Happé F
2. Frith U
(2006) The weak coherence account: detail-focused cognitive style in autism spectrum disorders
Journal of Autism and Developmental Disorders 36:5–25.

https://doi.org/10.1007/s10803-005-0039-0
- PubMed
- Google Scholar
1. Harris H
2. Israeli D
3. Minshew N
4. Bonneh Y
5. Heeger DJ
6. Behrmann M
7. Sagi D
(2015) Perceptual learning in autism: over-specificity and possible remedies
Nature Neuroscience 18:1574.

https://doi.org/10.1038/nn.4129
- PubMed
- Google Scholar
1. Hemsley DR
2. Garety PA
(1986) The formation of maintenance of delusions: a Bayesian analysis
The British Journal of Psychiatry 149:51–56.

https://doi.org/10.1192/bjp.149.1.51
- PubMed
- Google Scholar
1. Holland PW
2. Welsch RE
(1977) Robust regression using iteratively reweighted least-squares
Communications in Statistics - Theory and Methods 6:813–827.

https://doi.org/10.1080/03610927708827533
- Google Scholar
1. Huber PJ
(1964) Robust Estimation of a Location Parameter
The Annals of Mathematical Statistics 35:73–101.

https://doi.org/10.1214/aoms/1177703732
- Google Scholar
(2016) Central tendency effects in time interval reproduction in autism
Scientific Reports 6:28570.

https://doi.org/10.1038/srep28570
- PubMed
- Google Scholar
1. Kass RE
2. Raftery AE
(1995) Bayes Factors
Journal of the American Statistical Association 90:773–795.

https://doi.org/10.1080/01621459.1995.10476572
- Google Scholar
(2017) Adults with autism overestimate the volatility of the sensory environment
Nature Neuroscience 20:1293–1299.

https://doi.org/10.1038/nn.4615
- PubMed
- Google Scholar
(2014) An aberrant precision account of autism
Frontiers in Human Neuroscience 8:302.

https://doi.org/10.3389/fnhum.2014.00302
- PubMed
- Google Scholar
1. Ma WJ
2. Beck JM
3. Latham PE
4. Pouget A
(2006) Bayesian inference with probabilistic population codes
Nature Neuroscience 9:1432–1438.

https://doi.org/10.1038/nn1790
- PubMed
- Google Scholar
1. Manning C
2. Kilner J
3. Neil L
4. Karaminis T
5. Pellicano E
(2017) Children on the autism spectrum update their behaviour in response to a volatile environment
Developmental Science 20:e12435.

https://doi.org/10.1111/desc.12435
- Google Scholar
(2013) Evidence of a dimensional relationship between schizotypy and schizophrenia: a systematic review
Neuroscience & Biobehavioral Reviews 37:317–327.

https://doi.org/10.1016/j.neubiorev.2013.01.004
- PubMed
- Google Scholar
1. Notredame CE
2. Pins D
3. Deneve S
4. Jardri R
(2014) What visual illusions teach us about schizophrenia
Frontiers in Integrative Neuroscience 8:.

https://doi.org/10.3389/fnint.2014.00063
- PubMed
- Google Scholar
(2017) Bayesian approaches to autism: Towards volatility, action, and behavior
Psychological Bulletin 143:521–542.

https://doi.org/10.1037/bul0000097
- PubMed
- Google Scholar
(2017) The Importance of Falsification in Computational Cognitive Modeling
Trends in Cognitive Sciences 21:425–433.

https://doi.org/10.1016/j.tics.2017.03.011
- PubMed
- Google Scholar
(2016) Intact priors for gaze direction in adults with high-functioning autism spectrum conditions
Molecular Autism 7:25.

https://doi.org/10.1186/s13229-016-0085-9
- PubMed
- Google Scholar
1. Pellicano E
2. Burr D
(2012a) When the world becomes 'too real': a Bayesian explanation of autistic perception
Trends in Cognitive Sciences 16:504–510.

https://doi.org/10.1016/j.tics.2012.08.009
- PubMed
- Google Scholar
1. Pellicano E
2. Burr D
(2012b) Response to Brock: noise and autism
Trends in Cognitive Sciences 16:574–575.

https://doi.org/10.1016/j.tics.2012.10.004
- Google Scholar
1. Plaisted KC
(2015)
Reduced Generalization in Autism: An Alternative to Weak Central Coherence

Reduced Generalization in Autism: An Alternative to Weak Central Coherence.
- Google Scholar
(2016) Bayesian models of individual differences
Psychological Science 27:1562–1572.

https://doi.org/10.1177/0956797616665351
- PubMed
- Google Scholar
(2017) Pavlovian conditioning-induced hallucinations result from overweighting of perceptual priors
Science 357:596–600.

https://doi.org/10.1126/science.aan3458
- PubMed
- Google Scholar
1. Raine A
(1991) The SPQ: a scale for the assessment of schizotypal personality based on DSM-III-R criteria
Schizophrenia Bulletin 17:555–564.

https://doi.org/10.1093/schbul/17.4.555
- PubMed
- Google Scholar
1. Rust J
(1988) The rust inventory of schizotypal cognitions (RISC)
Schizophrenia Bulletin 14:317–322.

https://doi.org/10.1093/schbul/14.2.317
- PubMed
- Google Scholar
1. Sato Y
2. Kording KP
(2014) How much to trust the senses: likelihood learning
Journal of Vision 14:13.

https://doi.org/10.1167/14.13.13
- PubMed
- Google Scholar
(2013) Delusions and the role of beliefs in perceptual inference
Journal of Neuroscience 33:13701–13712.

https://doi.org/10.1523/JNEUROSCI.1778-13.2013
- PubMed
- Google Scholar
(2017) Enhanced predictive signalling in schizophrenia
Human Brain Mapping 38:1767–1779.

https://doi.org/10.1002/hbm.23480
- PubMed
- Google Scholar
(2015) Perceptual instability in schizophrenia: probing predictive coding accounts of delusions with ambiguous stimuli
Schizophrenia Research: Cognition 2:72–77.

https://doi.org/10.1016/j.scog.2015.03.005
- PubMed
- Google Scholar
1. Seriès P
2. Seitz AR
(2013) Learning what to expect (in visual perception)
Frontiers in Human Neuroscience 7:668.

https://doi.org/10.3389/fnhum.2013.00668
- PubMed
- Google Scholar
1. Skewes JC
2. Gebauer L
(2016) Brief report: suboptimal auditory localization in autism spectrum disorder: support for the bayesian account of sensory symptoms
Journal of Autism and Developmental Disorders 46:2539–2547.

https://doi.org/10.1007/s10803-016-2774-9
- PubMed
- Google Scholar
(2015) Perceptual inference and autistic traits
Autism 19:301–307.

https://doi.org/10.1177/1362361313519872
- PubMed
- Google Scholar
(2010) The contribution of hypersalience to the "jumping to conclusions" bias associated with delusions in schizophrenia
Journal of Psychiatry and Neuroscience 35:7–17.

https://doi.org/10.1503/jpn.090025
- PubMed
- Google Scholar
(2006) Synaptic plasticity and dysconnection in schizophrenia
Biological Psychiatry 59:929–939.

https://doi.org/10.1016/j.biopsych.2005.10.005
- PubMed
- Google Scholar
1. Stocker AA
2. Simoncelli EP
(2006)
Advances in Neural Information Processing Systems

1289–1296, Sensory adaptation within a Bayesian framework for perception, Advances in Neural Information Processing Systems.
- Google Scholar
Software
1. Team J
(2017) JASP, version 0.8.6
JASP.

https://jasp-stats.org/faq/how-do-i-cite-jasp/
1. Tennant R
2. Hiller L
3. Fishwick R
4. Platt S
5. Joseph S
6. Weich S
7. Parkinson J
8. Secker J
9. Stewart-Brown S
(2007) The Warwick-Edinburgh Mental Well-being Scale (WEMWBS): development and UK validation
Health and Quality of Life Outcomes 5:.

https://doi.org/10.1186/1477-7525-5-63
- PubMed
- Google Scholar
1. Teufel C
2. Subramaniam N
3. Dobler V
4. Perez J
5. Finnemann J
6. Mehta PR
7. Goodyer IM
8. Fletcher PC
(2015) Shift toward prior knowledge confers a perceptual advantage in early psychosis and psychosis-prone healthy individuals
PNAS 112:13401–13406.

https://doi.org/10.1073/pnas.1503916112
- PubMed
- Google Scholar
1. van Boxtel JJ
2. Lu H
(2013) A predictive coding perspective on autism spectrum disorders
Frontiers in Psychology 4:19.

https://doi.org/10.3389/fpsyg.2013.00019
- PubMed
- Google Scholar
(2013) Weak priors versus overfitting of predictions in autism: Reply to Pellicano and Burr (TICS, 2012)
i-Perception 4:95–97.

https://doi.org/10.1068/i0580ic
- PubMed
- Google Scholar
(2014) Precise minds in uncertain worlds: predictive coding in autism
Psychological Review 121:649–675.

https://doi.org/10.1037/a0037665
- PubMed
- Google Scholar
(2009) A systematic review and meta-analysis of the psychosis continuum: evidence for a psychosis proneness-persistence-impairment model of psychotic disorder
Psychological Medicine 39:179–195.

https://doi.org/10.1017/S0033291708003814
- PubMed
- Google Scholar

Article and author information

Author details

Povilas Karvelis

IANC, School of Informatics, University of Edinburgh, Edinburgh, United Kingdom

Contribution
Formal analysis, Investigation, Visualization, Writing—original draft

Competing interests
No competing interests declared
Aaron R Seitz

Department of Psychology, UC Riverside, Riverside, United States

Contribution
Software, Funding acquisition, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Stephen M Lawrie
1. Division of Psychiatry, University of Edinburgh, Edinburgh, United Kingdom
2. Patrick Wild Centre, University of Edinburgh, Edinburgh, United Kingdom
Contribution
Methodology, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2444-5675
Peggy Seriès

IANC, School of Informatics, University of Edinburgh, Edinburgh, United Kingdom

Contribution
Conceptualization, Supervision, Project administration, Writing—review and editing

For correspondence
pseries@inf.ed.ac.uk

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8580-7975

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.

Copyright

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

4,949

views
656

downloads
71

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Article PDF

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Povilas Karvelis
Aaron R Seitz
Stephen M Lawrie
Peggy Seriès

(2018)

Autistic traits, but not schizotypy, predict increased weighting of sensory information in Bayesian visual integration

eLife 7:e34115.

https://doi.org/10.7554/eLife.34115

Categories and tags

Research organism

Human

Share this article

Cite this article

Alternative hypotheses for ASD impairments within the Bayesian inference framework.

The moving dots task.

Average group performance on low-contrast trials (a–d) and on trials with no stimulus (e).

Figure 3—source data 1

Correlations between AQ scores and task performance on low contrast trials (a, b) and when no stimulus is presented (c).

Figure 4—source data 1

Bayesian model of estimation response for a single trial.

Modelling results.

Correlations between AQ scores and BAYES model parameters.

Figure 7—source data 1

Comparison of actual (x-axis) vs. recovered (y-axis) parameters using the BAYES’ model.

Task performance at the highest contrast level and exclusion Criteria.

Population averaged stimulus contrast relative to the background contrast for the 2/1 (red) and 4/1 (black) staircased contrast levels.

Cumulative moving average of ratio of estimation reaction times at ±32° vs average reaction times at all other directions.

Cumulative moving average of ratio of fraction of detected stimuli at ±32° vs average fraction detected at all other directions.

Cumulative moving average of ratio of fraction of detected stimuli at ±32° vs average fraction detected at all other directions.

Average group performance on low-contrast trials (a–d) and on trials with no stimulus (e) by groups split by median RISC score.

Average group performance on low-contrast trials (a–d) and on trials with no stimulus (e) by groups split by median SPQ score.

Correlations between personality traits, RISC (top row) and SPQ (bottom row) and task performance.

Correlations with the BAYES model parameter values and schizotypy traits (as measured by both RISC and SPQ).

A representative sample of prior expectations for each individual as reconstructed via ‘BAYES’ model.

Comparison of actual and recovered parameters via ‘BAYES_π’ model.

Task performance as predicted by the BAYES_full model.

Comparison of actual and recovered parameters via ‘BAYES_full’ model.

Author details

Povilas Karvelis

Contribution

Competing interests

Aaron R Seitz

Contribution

Competing interests

Stephen M Lawrie

Contribution

Competing interests

Peggy Seriès

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism

Further reading

Cumulative moving average of ratio of estimation reaction times at ±32^° vs average reaction times at all other directions.

Cumulative moving average of ratio of fraction of detected stimuli at ±32^° vs average fraction detected at all other directions.

Cumulative moving average of ratio of fraction of detected stimuli at ±32^° vs average fraction detected at all other directions.