The relationship between structure and function in the human brain is well established, but not yet well characterised. Large-scale biophysical models allow us to investigate this relationship, by ...leveraging structural information (e.g. derived from diffusion tractography) in order to couple dynamical models of local neuronal activity into networks of interacting regions distributed across the cortex. In practice however, these models are difficult to parametrise, and their simulation is often delicate and computationally expensive. This undermines the experimental aspect of scientific modelling, and stands in the way of comparing different parametrisations, network architectures, or models in general, with confidence. Here, we advocate the use of Bayesian optimisation for assessing the capabilities of biophysical network models, given a set of desired properties (e.g. band-specific functional connectivity); and in turn the use of this assessment as a principled basis for incremental modelling and model comparison. We adapt an optimisation method designed to cope with costly, high-dimensional, non-convex problems, and demonstrate its use and effectiveness. Using five parameters controlling key aspects of our model, we find that this method is able to converge to regions of high functional similarity with real MEG data, with very few samples given the number of parameters, without getting stuck in local extrema, and while building and exploiting a map of uncertainty defined smoothly across the parameter space. We compare the results obtained using different methods of structural connectivity estimation from diffusion tractography, and find that one method leads to better simulations.
•Biophysical models provide insight into the brain's structure-function interaction.•Assessing their ability to generate dynamics of interest is a difficult problem.•We propose and demonstrate the use of Bayesian optimisation for this purpose.
Schippers, Renken and Keysers (NeuroImage, 2011) present a simulation of multi-subject lag-based causality estimation. We fully agree that single-subject evaluations (e.g., Smith et al., 2011) need ...to be revisited in the context of multi-subject studies, and Schippers' paper is a good example, including detailed multi-level simulation and cross-subject statistical modelling. The authors conclude that “the average chance to find a significant Granger causality effect when no actual influence is present in the data stays well below the p-level imposed on the second level statistics” and that “when the analyses reveal a significant directed influence, this direction was accurate in the vast majority of the cases”. Unfortunately, we believe that the general meaning that may be taken from these statements is not supported by the paper's results, as there may in reality be a systematic (group-average) difference in haemodynamic delay between two brain areas. While many statements in the paper (e.g., the final two sentences) do refer to this problem, we fear that the overriding message that many readers may take from the paper could cause misunderstanding.
► Group-level FMRI simulations can be useful to test methods such as Granger causality. ► Simulation results need careful evaluation and interpretation. ► There is ample evidence of haemodynamic variability across regions and voxels. ► Lag-based FMRI causality analysis may be biassed by such variation. ► This confound should be considered when reporting lag-based results.
•We propose a registration model with inferred spatially adaptive regularisation.•The effects of this regularisation prior are shown on tensor based morphometry.•The inferred prior leads to better ...localisation of deformations.•We illustrate how this leads to a more realistic description of registration uncertainty.•We demonstrate how Bayesian model comparison can be used in registration.
Display omitted
This paper introduces a novel method for inferring spatially varying regularisation in non-linear registration. This is achieved through full Bayesian inference on a probabilistic registration model, where the prior on the transformation parameters is parameterised as a weighted mixture of spatially localised components. Such an approach has the advantage of allowing the registration to be more flexibly driven by the data than a traditional globally defined regularisation penalty, such as bending energy. The proposed method adaptively determines the influence of the prior in a local region. The strength of the prior may be reduced in areas where the data better support deformations, or can enforce a stronger constraint in less informative areas. Consequently, the use of such a spatially adaptive prior may reduce unwanted impacts of regularisation on the inferred transformation. This is especially important for applications where the deformation field itself is of interest, such as tensor based morphometry. The proposed approach is demonstrated using synthetic images, and with application to tensor based morphometry analysis of subjects with Alzheimer’s disease and healthy controls. The results indicate that using the proposed spatially adaptive prior leads to sparser deformations, which provide better localisation of regional volume change. Additionally, the proposed regularisation model leads to more data driven and localised maps of registration uncertainty. This paper also demonstrates for the first time the use of Bayesian model comparison for selecting different types of regularisation.
Our ability to hold information in mind is strictly limited. We sought to understand the relationship between oscillatory brain activity and the allocation of resources within visual short-term ...memory (VSTM). Participants attempted to remember target arrows embedded among distracters and used a continuous method of responding to report their memory for a cued target item. Trial-to-trial variability in the absolute circular accuracy with which participants could report the target was predicted by event-related alpha synchronization during initial processing of the memoranda and by alpha desynchronization during the retrieval of those items from VSTM. Using a model-based approach, we were also able to explore further which parameters of VSTM-guided behavior were most influenced by alpha band changes. Alpha synchronization during item processing enhanced the precision with which an item could be retained without affecting the likelihood of an item being represented per se (as indexed by the guessing rate). Importantly, our data outline a neural mechanism that mirrors the precision with which items are retained; the greater the alpha power enhancement during encoding, the greater the precision with which that item can be retained.
The techniques available for the interrogation and analysis of neuroimaging data have a large influence in determining the flexibility, sensitivity, and scope of neuroimaging experiments. The ...development of such methodologies has allowed investigators to address scientific questions that could not previously be answered and, as such, has become an important research area in its own right.
In this paper, we present a review of the research carried out by the Analysis Group at the Oxford Centre for Functional MRI of the Brain (FMRIB). This research has focussed on the development of new methodologies for the analysis of both structural and functional magnetic resonance imaging data. The majority of the research laid out in this paper has been implemented as freely available software tools within FMRIB's Software Library (FSL).
Functional magnetic resonance imaging studies often involve the acquisition of data from multiple sessions and/or multiple subjects. A hierarchical approach can be taken to modelling such data with a ...general linear model (GLM) at each level of the hierarchy introducing different random effects variance components. Inferring on these models is nontrivial with frequentist solutions being unavailable. A solution is to use a Bayesian framework. One important ingredient in this is the choice of prior on the variance components and top-level regression parameters. Due to the typically small numbers of sessions or subjects in neuroimaging, the choice of prior is critical. To alleviate this problem, we introduce to neuroimage modelling the approach of reference priors, which drives the choice of prior such that it is noninformative in an information-theoretic sense. We propose two inference techniques at the top level for multilevel hierarchies (a fast approach and a slower more accurate approach). We also demonstrate that we can infer on the top level of multilevel hierarchies by inferring on the levels of the hierarchy separately and passing summary statistics of a noncentral multivariate
t distribution between them.
Brainwave viscosity in propofol anaesthesia Fabus, M.S.; Woolrich, M.W.; Warnaby, C.E.
British journal of anaesthesia : BJA,
February 2022, 2022-Feb, 2022-02-00, 20220201, Letnik:
128, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Human EEG during propofol anaesthesia shows large-scale changes including traveling slow waves
. Slow-wave saturation is a potentially individualised marker of loss of perception
. However, much ...remains unclear about the dynamics of slow waves. Iterated empirical mode decomposition (itEMD
) is a novel data-driven method for segregating data into physiologically relevant oscillatory modes. We used itEMD to identify spectral modes and their sources / sinks in propofol EEG. Viscosity is a physical quantity expressing the magnitude of resistance to flow. Considering traveling electric potentials in the brain as a flow, we extended the notion of viscosity to traveling brainwaves. Using this, we explored how brainwave viscosity changes in volunteer propofol anaesthesia.
The therapeutic effects of centrally acting pharmaceuticals can manifest gradually and unreliably in patients, making the drug discovery process slow and expensive. Biological markers providing early ...evidence for clinical efficacy could help prioritize development of the more promising drug candidates. A potential source of such markers is functional magnetic resonance imaging (fMRI), a noninvasive imaging technique that can complement molecular imaging. fMRI has been used to characterize how drugs cause changes in brain activity. However, variation in study protocols and analysis techniques has made it difficult to identify consistent associations between subtle modulations of brain activity and clinical efficacy. We present and validate a general protocol for functional imaging-based assessment of drug activity in the central nervous system. The protocol uses machine learning methods and data from multiple published studies to identify reliable associations between drug-related activity modulations and drug efficacy, which can then be used to assess new data. A proof-of-concept version of this approach was developed and is shown here for analgesics (pain medication), and validated with eight separate studies of analgesic compounds. Our results show that the systematic integration of multistudy data permits the generalized inferences required for drug discovery. Multistudy integrative strategies of this type could help optimize the drug discovery and validation pipeline.
Mixture models are commonly used in the statistical segmentation of images. For example, they can be used for the segmentation of structural medical images into different matter types, or of ...statistical parametric maps into activating and nonactivating brain regions in functional imaging. Spatial mixture models have been developed to augment histogram information with spatial regularization using Markov random fields (MRFs). In previous work, an approximate model was developed to allow adaptive determination of the parameter controlling the strength of spatial regularization. Inference was performed using Markov Chain Monte Carlo (MCMC) sampling. However, this approach is prohibitively slow for large datasets. In this work, a more efficient inference approach is presented. This combines a variational Bayes approximation with a second-order Taylor expansion of the components of the posterior distribution, which would otherwise be intractable to Variational Bayes. This provides inference on fully adaptive spatial mixture models an order of magnitude faster than MCMC. We examine the behavior of this approach when applied to artificial data with different spatial characteristics, and to functional magnetic resonance imaging statistical parametric maps