Accurate modeling of temporal correlations in rapidly sampled fMRI time series

Rapid imaging techniques are increasingly used in functional MRI studies because they allow a greater number of samples to be acquired per unit time, thereby increasing statistical power. However, temporal correlations limit the increase in functional sensitivity and must be accurately accounted for to control the false‐positive rate. A common approach to accounting for temporal correlations is to whiten the data prior to estimating fMRI model parameters. Models of white noise plus a first‐order autoregressive process have proven sufficient for conventional imaging studies, but more elaborate models are required for rapidly sampled data. Here we show that when the “FAST” model implemented in SPM is used with a well‐controlled number of parameters, it can successfully prewhiten 80% of grey matter voxels even with volume repetition times as short as 0.35 s. We further show that the temporal signal‐to‐noise ratio (tSNR), which has conventionally been used to assess the relative functional sensitivity of competing imaging approaches, can be augmented to account for the temporal correlations in the time series. This amounts to computing the t‐score testing for the mean signal. We show in a visual perception task that unlike the tSNR weighted by the number of samples, the t‐score measure is directly related to the t‐score testing for activation when the temporal correlations are correctly modeled. This score affords a more accurate means of evaluating the functional sensitivity of different data acquisition options.

[1]  G. Box,et al.  Distribution of Residual Autocorrelations in Autoregressive-Integrated Moving Average Time Series Models , 1970 .

[2]  Klaus Scheffler,et al.  Effect of temporal resolution and serial autocorrelations in event‐related functional MRI , 2016, Magnetic resonance in medicine.

[3]  Catie Chang,et al.  Effects of model-based physiological noise correction on default mode network anti-correlations and correlations , 2009, NeuroImage.

[4]  Steen Moeller,et al.  Functional Sensitivity of 2D Simultaneous Multi-Slice Echo-Planar Imaging: Effects of Acceleration on g-factor and Physiological Noise , 2017, Front. Neurosci..

[5]  Guy B. Williams,et al.  Autocorrelation bias still exists in fMRI results , 2017 .

[6]  Steen Moeller,et al.  Evaluation of 2D multiband EPI imaging for high-resolution, whole-brain, task-based fMRI studies at 3T: Sensitivity and slice leakage artifacts , 2016, NeuroImage.

[7]  J. Polimeni,et al.  Blipped‐controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g‐factor penalty , 2012, Magnetic resonance in medicine.

[8]  Andrea Bergmann,et al.  Statistical Parametric Mapping The Analysis Of Functional Brain Images , 2016 .

[9]  R. Weisskoff,et al.  Effect of temporal autocorrelation due to physiological noise and stimulus paradigm on voxel‐level false‐positive rates in fMRI , 1998, Human brain mapping.

[10]  Mayur Narsude,et al.  Three‐dimensional echo planar imaging with controlled aliasing: A sequence for high temporal resolution functional MRI , 2016, Magnetic resonance in medicine.

[11]  Lawrence L. Wald,et al.  Three dimensional echo-planar imaging at 7 Tesla , 2010, NeuroImage.

[12]  G H Glover,et al.  Image‐based method for retrospective correction of physiological motion effects in fMRI: RETROICOR , 2000, Magnetic resonance in medicine.

[13]  P. Boesiger,et al.  SENSE: Sensitivity encoding for fast MRI , 1999, Magnetic resonance in medicine.

[14]  Peter M. Jakob,et al.  Controlled aliasing in parallel imaging , 2004 .

[15]  Karl J. Friston,et al.  Classical and Bayesian Inference in Neuroimaging: Theory , 2002, NeuroImage.

[16]  An Zi EFFECTIVE NUMBER OF OBSERVATIONS AND UNBIASED ESTIMATORS OF VARIANCE FOR AUTOCORRELATED DATA − AN OVERVIEW , 2010 .

[17]  Peter A. Bandettini,et al.  Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI , 2006, NeuroImage.

[18]  Kawin Setsompop,et al.  Interslice leakage artifact reduction technique for simultaneous multislice acquisitions , 2014, Magnetic resonance in medicine.

[19]  Robin M Heidemann,et al.  Controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA) for multi‐slice imaging , 2005, Magnetic resonance in medicine.

[20]  Nikolaus Weiskopf,et al.  Quantitative multi-parameter mapping of R1, PD*, MT, and R2* at 3T: a multi-center validation , 2013, Front. Neurosci..

[21]  Karl J. Friston,et al.  Analysis of fMRI Time-Series Revisited—Again , 1995, NeuroImage.

[22]  M. Fukunaga,et al.  Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: a 7 T study. , 2009, Magnetic resonance imaging.

[23]  D. Larkman,et al.  Use of multicoil arrays for separation of signal from multiple slices simultaneously excited , 2001, Journal of magnetic resonance imaging : JMRI.

[24]  Peter A. Bandettini,et al.  The respiration response function: The temporal dynamics of fMRI signal fluctuations related to changes in respiration , 2008, NeuroImage.

[25]  Daniel B. Rowe,et al.  Impacts of simultaneous multislice acquisition on sensitivity and specificity in fMRI , 2018, NeuroImage.

[26]  Mark W. Woolrich,et al.  Resting-state fMRI in the Human Connectome Project , 2013, NeuroImage.

[27]  H. Abdi The General Linear Model , 2009 .

[28]  Karl J. Friston,et al.  Variational free energy and the Laplace approximation , 2007, NeuroImage.

[29]  Robin M Heidemann,et al.  Generalized autocalibrating partially parallel acquisitions (GRAPPA) , 2002, Magnetic resonance in medicine.

[30]  Markus Barth,et al.  Serial correlations in single-subject fMRI with sub-second TR , 2017, NeuroImage.

[31]  Oliver Speck,et al.  The impact of physiological noise correction on fMRI at 7 T , 2011, NeuroImage.

[32]  Karl J. Friston,et al.  Classical and Bayesian Inference in Neuroimaging: Applications , 2002, NeuroImage.

[33]  Karl J. Friston,et al.  Unified segmentation , 2005, NeuroImage.

[34]  Kevin Murphy,et al.  How long to scan? The relationship between fMRI temporal signal to noise ratio and necessary scan duration , 2007, NeuroImage.

[35]  Dirk Ostwald,et al.  Variational Bayesian Parameter Estimation Techniques for the General Linear Model , 2017, bioRxiv.

[36]  Yunjie Tong,et al.  Studying the Spatial Distribution of Physiological Effects on BOLD Signals Using Ultrafast fMRI , 2014, Front. Hum. Neurosci..

[37]  Karl J. Friston,et al.  Variational Bayesian inference for fMRI time series , 2003, NeuroImage.

[38]  Hans Knutsson,et al.  Does Parametric Fmri Analysis with Spm Yield Valid Results? -an Empirical Study of 1484 Rest Datasets Does Parametric Fmri Analysis with Spm Yield Valid Results? - an Empirical Study of 1484 Rest Datasets , 2022 .