Effect of temporal autocorrelation due to physiological noise and stimulus paradigm on voxel‐level false‐positive rates in fMRI

Statistical mapping within a binary hypothesis testing framework is the most widely used analytical method in functional MRI of the brain. A common assumption in this kind of analysis is that the fMRI time series are independent and identically distributed in time, yet we know that fMRI data can have significant temporal correlation due to low‐frequency physiological fluctuation (Weisskoff et al. [1993]; Proc Soc Magn Reson Med 9:7; Biswal et al. [1995]: Mag Reson Med 34:537–541). Furthermore, since the signal‐to‐noise ratio will vary with imaging rate, we should expect that the degree of correlation will vary with imaging rate. In this paper, we investigate the effect of temporal correlation and experimental paradigm on false‐positive rates (type I error rates), using data synthesized through a simple autoregressive plus white‐noise model whose parameters were estimated from real data over a range of imaging rates. We demonstrate that actual false‐positive rates can be biased far above or below the assumed significance level α when temporal autocorrelation is ignored in a way that depends on both the degree of correlation as well as the paradigm frequency. Furthermore, we present a simple method, based on the noise model described above, for correcting such distortions, and relate this method to the extended general linear model of Worsley and Friston ([1995]: Neuroimage 2:173–181). Hum. Brain Mapping 6:239–249, 1998. © 1998 Wiley‐Liss, Inc.

[1]  John G. Proakis,et al.  Probability, random variables and stochastic processes , 1985, IEEE Trans. Acoust. Speech Signal Process..

[2]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[3]  Richard A. Davis,et al.  Time Series: Theory and Methods , 2013 .

[4]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[5]  Karl J. Friston,et al.  Analysis of functional MRI time‐series , 1994, Human Brain Mapping.

[6]  Karl J. Friston,et al.  Statistical parametric maps in functional imaging: A general linear approach , 1994 .

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

[8]  B. Biswal,et al.  Functional connectivity in the motor cortex of resting human brain using echo‐planar mri , 1995, Magnetic resonance in medicine.

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

[10]  J Xiong,et al.  Assessment and optimization of functional MRI analyses , 1996, Human brain mapping.

[11]  E. Bullmore,et al.  Statistical methods of estimation and inference for functional MR image analysis , 1996, Magnetic resonance in medicine.

[12]  P. J. Jennings,et al.  Time series analysis in the time domain and resampling methods for studies of functional magnetic resonance brain imaging , 1997, Human brain mapping.

[13]  M. D’Esposito,et al.  Empirical analyses of BOLD fMRI statistics. I. Spatially unsmoothed data collected under null-hypothesis conditions. , 1997, NeuroImage.

[14]  Scott L. Zeger,et al.  Non‐linear Fourier Time Series Analysis for Human Brain Mapping by Functional Magnetic Resonance Imaging , 1997 .