Functional magnetic resonance imaging (fMRI) aims to locate activated regions in human brains when specific tasks are performed. The conventional tool for analyzing fMRI data applies some variant of the linear model, which is restrictive in modeling assumptions. To yield more accurate prediction of the time-course behavior of neuronal responses, the semiparametric inference for the underlying hemodynamic response function is developed to identify significantly activated voxels. Under mild regularity conditions, we demonstrate that a class of the proposed semiparametric test statistics, based on the local linear estimation technique, follow χ2 distributions under null hypotheses for a number of useful hypotheses. Furthermore, the asymptotic power functions of the constructed tests are derived under the fixed and contiguous alternatives. Simulation evaluations and real fMRI data application suggest that the semiparametric inference procedure provides more efficient detection of activated brain areas than the popular imaging analysis tools AFNI and FSL.
References
Bickel, P. J. and Levina, E. (2008). Regularized estimation of large covariance matrices. Ann. Statist. 36 199–227.
Bickel, P. J. and Li, B. (2006). Regularization in statistics (with discussion). Test 15 271–344.
Benjamini, Y. and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. Roy. Statist. Soc. Ser. B 57 289–300.
Bosq, D. (1998). Nonparametric Statistics for Stochastic Processes: Estimation and Prediction, 2nd ed. Springer, Berlin.
Cox, R. W. (1996). AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. 29 162–73.
Demko, S., Moss, W. F. and Smith, P. W. (1984). Decay rates for inverses of band matrices. Math. Comp. 43 491–499.
Mathematical Reviews (MathSciNet):
MR758197
Fahrmeir, L. and Gössl, C. (2002). Semiparametric Bayesian models for human brain mapping. Statistical Modelling 2 235–250.
Fan, J. and Gijbels, I. (1996). Local Polynomial Modeling and Its Applications. Chapman and Hall, London.
Fan, J. and Yao, Q. (2003). Nonlinear Time Series: Nonparametric and Parametric Methods. Springer, New York.
Fan, J., Zhang, C. M. and Zhang, J. (2001). Generalized likelihood ratio statistics and Wilks phenomenon. Ann. Statist. 29 153–193.
Friston, K. J. et. al. (1997). SPM course notes. Available at http://www.fil.ion.ucl.ac.uk/spm/.
Genovese, C. R. (2000). A Bayesian time-course model for functional magnetic resonance imaging data (with discussion). J. Amer. Statist. Assoc. 95 691–703.
Glover, G. H. (1999). Deconvolution of impulse response in event-related BOLD fMRI. NeuroImage 9 416–429.
Golub, G. H. and Van Loan, C. F. (1996). Matrix Computations, 3rd ed. Johns Hopkins Univ. Press, Baltimore, MD.
Goutte, C., Nielsen, F. A. and Hansen, L. K. (2000). Modeling the haemodynamic response in fMRI using smooth FIR filters. IEEE Trans. Med. Imag. 19 1188–1201.
Josephs, O. and Henson, R. N. A. (1999). Event-related functional magnetic resonance imaging: Modelling, inference and optimization. Philos. Trans. Roy. Soc. 354 1215–1228.
Lange, N. (1996). Tutorial in biostatistics: Statistical approaches to human brain mapping by functional magnetic resonance imaging. Statistics in Medicine 15 389–428.
Lange, N. and Zeger, S. L. (1997). Non-linear Fourier times series analysis for human brain mapping by functional magnetic resonance imaging. Appl. Statist. 46 1–29.
Lazar, N. A, Eddy, W. F., Genovese, C. R. and Welling, J. (2001). Statistical issues in fMRI for brain imaging. Internat. Statist. Rev. 69 105–127.
Lu, Y. (2006). Contributions to functional data analysis with biological applications. Ph.D. dissertation, Dept. Statistics, Univ. Wisconsin, Madison.
Purdon, P. L., Solo, V., Weissko, R. M. and Brown, E. (2001). Locally regularized spatio temporal modeling and model comparison for functional MRI. NeuroImage 14 912–923.
Shao, J. (2003). Mathematical Statistics, 2nd ed. Springer, New York.
Silverman, B. W. (1986). Density Estimation For Statistics and Data Analysis. Chapman and Hall, London.
Mathematical Reviews (MathSciNet):
MR848134
Smith, S., Jenkinson, M., Woolrich, M., Beckmann, C. F., Behrens, T. E. J., Johansen-Berg, H., Bannister, P. R., De Luca, M., Drobnjak, I. Flitney, D. E., Niazy, R. K., Saunders, J., Vickers, J., Zhang, Y., De Stefano, N., Brady, J. M. and Matthews, P. M. (2004). Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage 23 208–219.
Storey, J. D. (2002). A direct approach to false discovery rates. J. Roy. Statist. Soc. Ser. B 64 479–498.
Ward, B. D. (2001). Deconvolution analysis of fMRI time series data. Technical report, Biophysics Research Institute, Medical College of Wisconsin.
Woolrich, M. W., Ripley, B. D., Brady, M. and Smith, S. M. (2001). Temporal autocorrelation in univariate linear modelling of fMRI data. NeuroImage 14 1370–1386.
Worsley, K. J. and Friston, K. J. (1995). Analysis of fMRI time-series revisited-again. NeuroImage 2 173–181.
Worsley, K. J., Liao, C. H., Aston, J., Petre, V., Duncan, G. H., Morales, F. and Evans, A. C. (2002). A general statistical analysis for fMRI data. NeuroImage 15 1–15.
Zhang, C. M. (2003). Calibrating the degrees of freedom for automatic data smoothing and effective curve checking. J. Amer. Statist. Assoc. 98 609–628.
Zhang, C. M. and Dette, H. (2004). A power comparison between nonparametric regression tests. Statist. Probab. Lett. 66 289–301.
Zhang, C. M., Lu, Y., Johnstone, T., Oakes, T. and Davidson, R. (2006). Efficient modeling and inference for event-related fMRI data. Technical report #1125, Dept. Statistics, Univ. Wisconsin-Madison.