Skip to article frontmatterSkip to article content

Supplementary Figures

Explained variance in state energy by first two principal components. See supplemental_material.ipynb for details.

Figure 1:Explained variance in state energy by first two principal components. See supplemental_material.ipynb for details.

Cross-validation classification accuracy of the fcHNN, when predicting the attractor state from state
activation.  See supplemental_material.ipynb for details.

Figure 2:Cross-validation classification accuracy of the fcHNN, when predicting the attractor state from state activation. See supplemental_material.ipynb for details.

Parameter sweep of fcHNN parameters threshold and beta. the number of attractor states is color-coded. See supplemental_material.ipynb for details.

Figure 3:Parameter sweep of fcHNN parameters threshold and beta. the number of attractor states is color-coded. See supplemental_material.ipynb for details.

HNNs initialized with the empirical connectome has better convergence properties than permutation-based null models. Histograms show the number of iterations until convergence with the real (unpermuted) connectome, for 100 random initializations, with various beta values. HNN models based on the permuted connectivity matrix (with retaining symmetry) did not reach convergence in 10000 iterations in more than 98% of the same random initializations. See convergence-analysis.ipynb for details.

Figure 4:HNNs initialized with the empirical connectome has better convergence properties than permutation-based null models. Histograms show the number of iterations until convergence with the real (unpermuted) connectome, for 100 random initializations, with various beta values. HNN models based on the permuted connectivity matrix (with retaining symmetry) did not reach convergence in 10000 iterations in more than 98% of the same random initializations. See convergence-analysis.ipynb for details.

Statistical inference of the fcHNN state occupancy prediction with different null models.
A Results with a spatial autocorrelation-preserving null model for the empirical activity patterns. See null_models.ipynb for more details.
B Results where simulated samples are randomly sampled from a multivariate normal distribution, with the functional connectome as the covariance matrix, and compared to the fcHNN performance. See supplemental_material.ipynb for details.

Figure 5:Statistical inference of the fcHNN state occupancy prediction with different null models. A Results with a spatial autocorrelation-preserving null model for the empirical activity patterns. See null_models.ipynb for more details. B Results where simulated samples are randomly sampled from a multivariate normal distribution, with the functional connectome as the covariance matrix, and compared to the fcHNN performance. See supplemental_material.ipynb for details.

FcHNN can reconstruct the pain “ghost attractor”.
Signal-to-noise values range from 0.003 to 0.009. Asterisk denotes the location of the simulated “ghost attractor”. P-values are based on permutation testing, by randomly changing the conditions in a per-participant basis. See main_analyses.ipynb for more details.

Figure 6:FcHNN can reconstruct the pain “ghost attractor”. Signal-to-noise values range from 0.003 to 0.009. Asterisk denotes the location of the simulated “ghost attractor”. P-values are based on permutation testing, by randomly changing the conditions in a per-participant basis. See main_analyses.ipynb for more details.

FcHNN can reconstruct the changes in brain dynamics caused by the voluntary donwregulation of pain (as contrasted to upregulation)
Signal-to-noise values range from 0.001 to 0.005. P-values are based on permutation testing, by randomly changing the conditions in a per-participant basis. See main_analyses.ipynb for more details.

Figure 7:FcHNN can reconstruct the changes in brain dynamics caused by the voluntary donwregulation of pain (as contrasted to upregulation) Signal-to-noise values range from 0.001 to 0.005. P-values are based on permutation testing, by randomly changing the conditions in a per-participant basis. See main_analyses.ipynb for more details.

Robustness of the fcHNN weights to noise.
We set the temperature of the fcHNN, so that two attractor states emerge and iteratively add noise to the connectome.
To account for the change in dynamics, we adjust the temperature (beta) of the noisy fcHNN so that exactly two states emerge. We then highlight the decrease in nodal strength of the noisy connectome (the fcHNN weights) as a reference metric
vs the correlation of the attractor states that emerge from the noisy connectome. See supplemental_material.ipynb for details.

Figure 8:Robustness of the fcHNN weights to noise. We set the temperature of the fcHNN, so that two attractor states emerge and iteratively add noise to the connectome. To account for the change in dynamics, we adjust the temperature (beta) of the noisy fcHNN so that exactly two states emerge. We then highlight the decrease in nodal strength of the noisy connectome (the fcHNN weights) as a reference metric vs the correlation of the attractor states that emerge from the noisy connectome. See supplemental_material.ipynb for details.

All significant differences of the mean state activation analysis on the ABIDE dataset; label denotes the region
in the BASC122 atlas. See supplemental_material.ipynb for details.

Figure 9:All significant differences of the mean state activation analysis on the ABIDE dataset; label denotes the region in the BASC122 atlas. See supplemental_material.ipynb for details.

Supplementary Tables

Table 1:Neurosynth meta-analyses. The table includes details about the term used for the automated meta-analyses, as well as the number of studies included in the meta-analysis, the total number of reported activations and the maximal Z-statistic from the meta-analysis.

search termnum. studiesnum. activationsmax. Z
pain5162329514.8
motor256510949122.5
auditory12524655725.3
visual311011572615.4
face8963184226.8
autobiographical143725115.7
theory of mind181776115.1
sentences3561346116.5

Supplementary Methods

Study 4 instructions for upregulation. “During this scan, we are going to ask you to try to imagine as hard as you can that the thermal stimulations are more painful than they are. Try to focus on how unpleasant the pain is, for instance, how strongly you would like to remove your arm from it. Pay attention to the burning, stinging and shooting sensations. You can use your mind to turn up the dial of the pain, much like turning up the volume dial on a stereo. As you feel the pain rise in intensity, imagine it rising faster and faster and going higher and higher. Picture your skin being held up against a glowing hot metal or fire. Think of how disturbing it is to be burned, and visualize your skin sizzling, melting and bubbling as a result of the intense heat.”

Study 4 instructions for downregulation. “During this scan, we are going to ask you to try to imagine as hard as you can that the thermal stimulations are less painful than they are. Focus on the part of the sensation that is pleasantly warm, like a blanket on a cold day. You can use your mind to turn down the dial of your pain sensation, much like turning down the volume dial on a stereo. As you feel the stimulation rise, let it numb your arm, so any pain you feel simply fades away. Imagine your skin is very cool, from being outside, and think of how good the stimulation feels as it warms you up.”