A Bayesian time-course model for functional magnetic resonance imaging data

Functional magnetic resonance imaging (fMRI) is a new technique for studying the workings of the active human brain. During an fMRI experiment, a sequence of magnetic resonance images is acquired while the subject performs specific behavioral tasks. Changes in the measured signal can be used to identify and characterize the brain activity resulting from task performance. The data obtained from an fMRI experiment are a realization of a complex spatiotemporal process with many sources of variation, both biological and technological. This article describes a nonlinear Bayesian hierarchical model for fMRI data and presents inferential methods that enable investigators to directly target their scientific questions of interest, many of which are inaccessible to current methods. The article describes optimization and posterior sampling techniques to fit the model, both of which must be applied many thousands of times for a single dataset. The model is used to analyze data from a psychological experiment and to test a specific prediction of a cognitive theory.