Predicting brain activity using a Bayesian spatial model.
Ontology highlight
ABSTRACT:
SUBMITTER: Derado G
PROVIDER: S-EPMC4175991 | biostudies-literature | 2013 Aug
REPOSITORIES: biostudies-literature
Publications
Predicting brain activity using a Bayesian spatial model.
Derado Gordana G Bowman F Dubois FD Zhang Lijun L
Statistical methods in medical research 20120628 4
Increasing the clinical applicability of functional neuroimaging technology is an emerging objective, e.g. for diagnostic and treatment purposes. We propose a novel Bayesian spatial hierarchical framework for predicting follow-up neural activity based on an individual's baseline functional neuroimaging data. Our approach attempts to overcome some shortcomings of the modeling methods used in other neuroimaging settings, by borrowing strength from the spatial correlations present in the data. Our ...[more]
Similar Datasets
Project description:In vivo functional neuroimaging, including functional magnetic resonance imaging (fMRI) and positron emission tomography (PET), is becoming increasingly important in defining the pathophysiology of psychiatric disorders such as schizophrenia, major depression, and Alzheimer's disease. Furthermore, recent studies have begun to investigate the possibility of using functional neuroimaging to guide treatment selection for individual patients. By studying the changes between a patient's pre- and post-treatment brain activity, investigators are gaining insights into the impact of treatment on behavior-related neural processing traits associated with particular psychiatric disorders. Furthermore, these studies may shed light on the neural basis of response and nonresponse to specific treatments. The practical limitation of such studies is that the post-treatment scans offer little guidance to treatment selection in clinical settings, since treatment decisions precede the availability of post-treatment brain scans. This shortcoming represents the impetus for developing statistical methodology that would provide clinicians with predictive information concerning the effect of treatment on brain activity and, ultimately, symptom-related behaviors. We present a prediction algorithm that uses a patient's pretreatment scans, coupled with relevant patient characteristics, to forecast the patient's brain activity following a specified treatment regimen. We derive our predictive method from a Bayesian hierarchical model constructed on the pre- and post-treatment scans of designated training data. We perform estimation using the expectation-maximization algorithm. We evaluate the accuracy of our proposed prediction method using K-fold cross-validation, quantifying the error using two new measures that we propose for neuroimaging data. The proposed method is applicable to both PET and fMRI studies. We illustrate its use with a PET study of working memory in patients with schizophrenia and an fMRI data example is also provided.
| S-EPMC2757106 | biostudies-literature
Project description:Environmental health exposures to airborne chemicals often originate from chemical mixtures. Environmental health professionals may be interested in assessing exposure to one or more of the chemicals in these mixtures, but often exposure measurement data are not available, either because measurements were not collected/assessed for all exposure scenarios of interest or because some of the measurements were below the analytical methods' limits of detection (i.e. censored). In some cases, based on chemical laws, two or more components may have linear relationships with one another, whether in a single or in multiple mixtures. Although bivariate analyses can be used if the correlation is high, often correlations are low. To serve this need, this paper develops a multivariate framework for assessing exposure using relationships of the chemicals present in these mixtures. This framework accounts for censored measurements in all chemicals, allowing us to develop unbiased exposure estimates. We assessed our model's performance against simpler models at a variety of censoring levels and assessed our model's 95% coverage. We applied our model to assess vapor exposure from measurements of three chemicals in crude oil taken on the Ocean Intervention III during the Deepwater Horizon oil spill response and clean-up.
| S-EPMC6241297 | biostudies-literature
Project description:Streamlined Quantitative BOLD (sqBOLD) is an MR technique that can non-invasively measure physiological parameters including Oxygen Extraction Fraction (OEF) and deoxygenated blood volume (DBV) in the brain. Current sqBOLD methodology rely on fitting a linear model to log-transformed data acquired using an Asymmetric Spin Echo (ASE) pulse sequence. In this paper, a non-linear model implemented in a Bayesian framework was used to fit physiological parameters to ASE data. This model makes use of the full range of available ASE data, and incorporates the signal contribution from venous blood, which was ignored in previous analyses. Simulated data are used to demonstrate the intrinsic difficulty in estimating OEF and DBV simultaneously, and the benefits of the proposed non-linear model are shown. In vivo data are used to show that this model improves parameter estimation when compared with literature values. The model and analysis framework can be extended in a number of ways, and can incorporate prior information from external sources, so it has the potential to further improve OEF estimation using sqBOLD.
| S-EPMC7334042 | biostudies-literature
Project description:Relating disease status to imaging data stands to increase the clinical significance of neuroimaging studies. Many neurological and psychiatric disorders involve complex, systems-level alterations that manifest in functional and structural properties of the brain and possibly other clinical and biologic measures. We propose a Bayesian hierarchical model to predict disease status, which is able to incorporate information from both functional and structural brain imaging scans. We consider a two-stage whole brain parcellation, partitioning the brain into 282 subregions, and our model accounts for correlations between voxels from different brain regions defined by the parcellations. Our approach models the imaging data and uses posterior predictive probabilities to perform prediction. The estimates of our model parameters are based on samples drawn from the joint posterior distribution using Markov Chain Monte Carlo (MCMC) methods. We evaluate our method by examining the prediction accuracy rates based on leave-one-out cross validation, and we employ an importance sampling strategy to reduce the computation time. We conduct both whole-brain and voxel-level prediction and identify the brain regions that are highly associated with the disease based on the voxel-level prediction results. We apply our model to multimodal brain imaging data from a study of Parkinson's disease. We achieve extremely high accuracy, in general, and our model identifies key regions contributing to accurate prediction including caudate, putamen, and fusiform gyrus as well as several sensory system regions.
| S-EPMC5879954 | biostudies-literature
Project description:In this paper, we propose a spatial lag panel smoothing transition regression (SLPSTR) model ty considering spatial correlation of dependent variable in panel smooth transition regression model. This model combines advantages of both smooth transition model and spatial econometric model and can be used to deal with panel data with wide range of heterogeneity and cross-section correlation simultaneously. We also propose a Bayesian estimation approach in which the Metropolis-Hastings algorithm and the method of Gibbs are used for sampling design for SLPSTR model. A simulation study and a real data study are conducted to investigate the performance of the proposed model and the Bayesian estimation approach in practice. The results indicate that our theoretical method is applicable to spatial data with a wide range of spatial structures under finite sample.
| S-EPMC6398831 | biostudies-literature
Project description:A vast amount of ecological knowledge generated over the past two decades has hinged upon the ability of model selection methods to discriminate among various ecological hypotheses. The last decade has seen the rise of Bayesian hierarchical models in ecology. Consequently, commonly used tools, such as the AIC, become largely inapplicable and there appears to be no consensus about a particular model selection tool that can be universally applied. We focus on a specific class of competing Bayesian spatial capture-recapture (SCR) models and apply and evaluate some of the recommended Bayesian model selection tools: (1) Bayes Factor-using (a) Gelfand-Dey and (b) harmonic mean methods, (2) Deviance Information Criterion (DIC), (3) Watanabe-Akaike's Information Criterion (WAIC) and (4) posterior predictive loss criterion. In all, we evaluate 25 variants of model selection tools in our study. We evaluate these model selection tools from the standpoint of selecting the "true" model and parameter estimation. In all, we generate 120 simulated data sets using the true model and assess the frequency with which the true model is selected and how well the tool estimates N (population size), a parameter of much importance to ecologists. We find that when information content is low in the data, no particular model selection tool can be recommended to help realize, simultaneously, both the goals of model selection and parameter estimation. But, in general (when we consider both the objectives together), we recommend the use of our application of the Bayes Factor (Gelfand-Dey with MAP approximation) for Bayesian SCR models. Our study highlights the point that although new model selection tools are emerging (e.g., WAIC) in the applied statistics literature, those tools based on sound theory even under approximation may still perform much better.
| S-EPMC6822056 | biostudies-literature
Project description:BACKGROUND:Maps of influenza activity are important tools to monitor influenza epidemics and inform policymakers. In France, the availability of a high-quality data set from the Oscour® surveillance network, covering 92% of hospital emergency department (ED) visits, offers new opportunities for disease mapping. Traditional geostatistical mapping methods such as Kriging ignore underlying population sizes, are not suited to non-Gaussian data and do not account for uncertainty in parameter estimates. OBJECTIVE:Our objective was to create reliable weekly interpolated maps of influenza activity in the ED setting, to inform Santé publique France (the French national public health agency) and local healthcare authorities. METHODS:We used Oscour® data of ED visits covering the 2016-2017 influenza season. We developed a Bayesian model-based geostatistical approach, a class of generalized linear mixed models, with a multivariate normal random field as a spatially autocorrelated random effect. Using R-INLA, we developed an algorithm to create maps of the proportion of influenza-coded cases among all coded visits. We compared our results with maps obtained by Kriging. RESULTS:Over the study period, 45 565 (0.82%) visits were coded as influenza cases. Maps resulting from the model are presented for each week, displaying the posterior mean of the influenza proportion and its associated uncertainty. Our model performed better than Kriging. CONCLUSIONS:Our model allows producing smoothed maps where the random noise has been properly removed to reveal the spatial risk surface. The algorithm was incorporated into the national surveillance system to produce maps in real time and could be applied to other diseases.
| S-EPMC6185885 | biostudies-other
Project description:Daily flight activity patterns of forest insects are influenced by temporal and meteorological conditions. Temperature and time of day are frequently cited as key drivers of activity; however, complex interactions between multiple contributing factors have also been proposed. Here, we report individual Bayesian network models to assess the probability of flight activity of three exotic insects, Hylurgus ligniperda, Hylastes ater, and Arhopalus ferus in a managed plantation forest context. Models were built from 7,144 individual hours of insect sampling, temperature, wind speed, relative humidity, photon flux density, and temporal data. Discretized meteorological and temporal variables were used to build naïve Bayes tree augmented networks. Calibration results suggested that the H. ater and A. ferus Bayesian network models had the best fit for low Type I and overall errors, and H. ligniperda had the best fit for low Type II errors. Maximum hourly temperature and time since sunrise had the largest influence on H. ligniperda flight activity predictions, whereas time of day and year had the greatest influence on H. ater and A. ferus activity. Type II model errors for the prediction of no flight activity is improved by increasing the model's predictive threshold. Improvements in model performance can be made by further sampling, increasing the sensitivity of the flight intercept traps, and replicating sampling in other regions. Predicting insect flight informs an assessment of the potential phytosanitary risks of wood exports. Quantifying this risk allows mitigation treatments to be targeted to prevent the spread of invasive species via international trade pathways.
| S-EPMC5617153 | biostudies-literature
Project description:High throughput spatial transcriptomics (HST) is a rapidly emerging class of experimental technologies that allow for profiling gene expression in tissue samples at or near single-cell resolution while retaining the spatial location of each sequencing unit within the tissue sample. Through analyzing HST data, we seek to identify sub-populations of cells within a tissue sample that may inform biological phenomena. Existing computational methods either ignore the spatial heterogeneity in gene expression profiles, fail to account for important statistical features such as skewness, or are heuristic-based network clustering methods that lack the inferential benefits of statistical modeling. To address this gap, we develop SPRUCE: a Bayesian spatial multivariate finite mixture model based on multivariate skew-normal distributions, which is capable of identifying distinct cellular sub-populations in HST data. We further implement a novel combination of Pólya-Gamma data augmentation and spatial random effects to infer spatially correlated mixture component membership probabilities without relying on approximate inference techniques. Via a simulation study, we demonstrate the detrimental inferential effects of ignoring skewness or spatial correlation in HST data. Using publicly available human brain HST data, SPRUCE outperforms existing methods in recovering expertly annotated brain layers. Finally, our application of SPRUCE to human breast cancer HST data indicates that SPRUCE can distinguish distinct cell populations within the tumor microenvironment. An R package spruce for fitting the proposed models is available through The Comprehensive R Archive Network.
| S-EPMC10134739 | biostudies-literature
Project description:BackgroundAfrican trypanosomiasis is a tsetse-borne parasitic infection that affects humans, wildlife, and domesticated animals. Tsetse flies are endemic to much of Sub-Saharan Africa and a spatial and temporal understanding of tsetse habitat can aid surveillance and support disease risk management. Problematically, current fine spatial resolution remote sensing data are delivered with a temporal lag and are relatively coarse temporal resolution (e.g., 16 days), which results in disease control models often targeting incorrect places. The goal of this study was to devise a heuristic for identifying tsetse habitat (at a fine spatial resolution) into the future and in the temporal gaps where remote sensing and proximal data fail to supply information.MethodsThis paper introduces a generalizable and scalable open-access version of the tsetse ecological distribution (TED) model used to predict tsetse distributions across space and time, and contributes a geospatial Bayesian Maximum Entropy (BME) prediction model trained by TED output data to forecast where, herein the Morsitans group of tsetse, persist in Kenya, a method that mitigates the temporal lag problem. This model facilitates identification of tsetse habitat and provides critical information to control tsetse, mitigate the impact of trypanosomiasis on vulnerable human and animal populations, and guide disease minimization in places with ephemeral tsetse. Moreover, this BME analysis is one of the first to utilize cluster and parallel computing along with a Monte Carlo analysis to optimize BME computations. This allows for the analysis of an exceptionally large dataset (over 2 billion data points) at a finer resolution and larger spatiotemporal scale than what had previously been possible.ResultsUnder the most conservative assessment for Kenya, the BME kriging analysis showed an overall prediction accuracy of 74.8% (limited to the maximum suitability extent). In predicting tsetse distribution outcomes for the entire country the BME kriging analysis was 97% accurate in its forecasts.ConclusionsThis work offers a solution to the persistent temporal data gap in accurate and spatially precise rainfall predictions and the delayed processing of remotely sensed data collectively in the - 45 days past to + 180 days future temporal window. As is shown here, the BME model is a reliable alternative for forecasting future tsetse distributions to allow preplanning for tsetse control. Furthermore, this model provides guidance on disease control that would otherwise not be available. These 'big data' BME methods are particularly useful for large domain studies. Considering that past BME studies required reduction of the spatiotemporal grid to facilitate analysis. Both the GEE-TED and the BME libraries have been made open source to enable reproducibility and offer continual updates into the future as new remotely sensed data become available.
| S-EPMC10655428 | biostudies-literature