20 research outputs found
Computer model calibration with large non-stationary spatial outputs: application to the calibration of a climate model
Bayesian calibration of computer models tunes unknown input parameters by
comparing outputs with observations. For model outputs that are distributed
over space, this becomes computationally expensive because of the output size.
To overcome this challenge, we employ a basis representation of the model
outputs and observations: we match these decompositions to carry out the
calibration efficiently. In the second step, we incorporate the non-stationary
behaviour, in terms of spatial variations of both variance and correlations, in
the calibration. We insert two integrated nested Laplace
approximation-stochastic partial differential equation parameters into the
calibration. A synthetic example and a climate model illustration highlight the
benefits of our approach
A Bayesian General Linear Modeling Approach to Cortical Surface fMRI Data Analysis
Cortical surface functional magnetic resonance imaging (cs-fMRI) has recently grown in popularity versus traditional volumetric fMRI. In addition to offering better whole-brain visualization, dimension reduction, removal of extraneous tissue types, and improved alignment of cortical areas across subjects, it is also more compatible with common assumptions of Bayesian spatial models. However, as no spatial Bayesian model has been proposed for cs-fMRI data, most analyses continue to employ the classical general linear model (GLM), a âmassive univariateâ approach. Here, we propose a spatial Bayesian GLM for cs-fMRI, which employs a class of sophisticated spatial processes to model latent activation fields. We make several advances compared with existing spatial Bayesian models for volumetric fMRI. First, we use integrated nested Laplacian approximations, a highly accurate and efficient Bayesian computation technique, rather than variational Bayes. To identify regions of activation, we utilize an excursions set method based on the joint posterior distribution of the latent fields, rather than the marginal distribution at each location. Finally, we propose the first multi-subject spatial Bayesian modeling approach, which addresses a major gap in the existing literature. The methods are very computationally advantageous and are validated through simulation studies and two task fMRI studies from the Human Connectome Project. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement