A Bayesian General Linear Modeling Approach to Cortical Surface fMRI Data Analysis

Amanda F. Mejia, Yu (Ryan) Yue, David Bolin, Finn Lindgren, Martin Lindquist

Research output: Contribution to journalArticle

Abstract

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.

Original languageEnglish (US)
JournalJournal of the American Statistical Association
DOIs
StatePublished - Jan 1 2019

Fingerprint

Functional Magnetic Resonance Imaging
Data analysis
Spatial Model
Bayesian Model
Modeling
Activation
Linear Model
Variational Bayes
Bayesian Computation
Spatial Modeling
Bayesian Modeling
Spatial Process
Excursion
Dimension Reduction
Marginal Distribution
Posterior distribution
Joint Distribution
Univariate
Functional magnetic resonance imaging
Alignment

Keywords

  • Bayesian smoothing
  • Brain imaging
  • Integrated nested Laplace approximation
  • Spatial statistics
  • Stochastic partial differential equation

ASJC Scopus subject areas

  • Statistics and Probability
  • Statistics, Probability and Uncertainty

Cite this

A Bayesian General Linear Modeling Approach to Cortical Surface fMRI Data Analysis. / Mejia, Amanda F.; Yue, Yu (Ryan); Bolin, David; Lindgren, Finn; Lindquist, Martin.

In: Journal of the American Statistical Association, 01.01.2019.

Research output: Contribution to journalArticle

@article{3514a31bf8684851bdf2f91c97a60580,
title = "A Bayesian General Linear Modeling Approach to Cortical Surface fMRI Data Analysis",
abstract = "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.",
keywords = "Bayesian smoothing, Brain imaging, Integrated nested Laplace approximation, Spatial statistics, Stochastic partial differential equation",
author = "Mejia, {Amanda F.} and Yue, {Yu (Ryan)} and David Bolin and Finn Lindgren and Martin Lindquist",
year = "2019",
month = "1",
day = "1",
doi = "10.1080/01621459.2019.1611582",
language = "English (US)",
journal = "Journal of the American Statistical Association",
issn = "0162-1459",
publisher = "Taylor and Francis Ltd.",

}

TY - JOUR

T1 - A Bayesian General Linear Modeling Approach to Cortical Surface fMRI Data Analysis

AU - Mejia, Amanda F.

AU - Yue, Yu (Ryan)

AU - Bolin, David

AU - Lindgren, Finn

AU - Lindquist, Martin

PY - 2019/1/1

Y1 - 2019/1/1

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

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

KW - Bayesian smoothing

KW - Brain imaging

KW - Integrated nested Laplace approximation

KW - Spatial statistics

KW - Stochastic partial differential equation

UR - http://www.scopus.com/inward/record.url?scp=85067548935&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85067548935&partnerID=8YFLogxK

U2 - 10.1080/01621459.2019.1611582

DO - 10.1080/01621459.2019.1611582

M3 - Article

JO - Journal of the American Statistical Association

JF - Journal of the American Statistical Association

SN - 0162-1459

ER -