Testing for signals with unknown location and scale in a χ2 random field, with an application to fMRI

Siegmund and Worsley (1995) considered the problem of testing for signals with unknown location and scale in a Gaussian random field defined on ℝ N . The test statistic was the maximum of a Gaussian random field in an N+1 dimensional ‘scale space’, N dimensions for location and 1 dimension for the scale of a smoothing filter. Scale space is identical to a continuous wavelet transform with a kernel smoother as the wavelet, though the emphasis here is on signal detection rather than image compression or enhancement. Two methods were used to derive an approximate null distribution for N=2 and N=3: one based on the method of volumes of tubes, the other based on the expected Euler characteristic of the excursion set. The purpose of this paper is two-fold: to show how the latter method can be extended to higher dimensions, and to apply this more general result to χ2 fields. The result of Siegmund and Worsley (1995) then follows as a special case. In this paper the results are applied to the problem of searching for activation in brain images obtained by functional magnetic resonance imaging (fMRI).

[1]  Azriel Rosenfeld,et al.  Digital Picture Processing , 1976 .

[2]  K. Worsley.  Confidence regions and tests for a change-point in a sequence of exponential family random variables , 1986 .

[3]  R. Turner,et al.  Dynamic magnetic resonance imaging of human brain activity during primary sensory stimulation. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[4]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[5]  Alan C. Evans,et al.  A Three-Dimensional Statistical Analysis for CBF Activation Studies in Human Brain , 1992, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[6]  Donald E. Knuth.  Two notes on notation , 1992 .

[7]  Ravi S. Menon,et al.  Intrinsic signal changes accompanying sensory stimulation: functional brain mapping with magnetic resonance imaging. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Adrian T. Lee,et al.  fMRI of human visual cortex , 1994, Nature.

[9]  Jean-Baptiste Poline,et al.  Analysis of individual brain activation maps using hierarchical description and multiscale detection , 1994, IEEE Trans. Medical Imaging.

[10]  Sergio N. Torres.  Topological Analysis of COBE-DMR Cosmic Microwave Background Maps , 1994 .

[11]  J. Hogg.  Magnetic resonance imaging. , 1994, Journal of the Royal Naval Medical Service.

[12]  Changbom Park,et al.  Topological analysis of the CfA redshift survey , 1994 .

[13]  K. Worsley,et al.  Local Maxima and the Expected Euler Characteristic of Excursion Sets of χ 2, F and t Fields , 1994, Advances in Applied Probability.

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

[15]  J B Poline,et al.  Enhanced Detection in Brain Activation Maps Using a Multifiltering Approach , 1994, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

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

[17]  D. Siegmund,et al.  Testing for a Signal with Unknown Location and Scale in a Stationary Gaussian Random Field , 1995 .

[18]  K. Worsley,et al.  Boundary corrections for the expected Euler characteristic of excursion sets of random fields, with an application to astrophysics , 1995, Advances in Applied Probability.

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

[20]  J. B. Poline,et al.  Intersubject Variability in Functional Neuroanatomy of Silent Verb Generation: Assessment by a New Activation Detection Algorithm Based on Amplitude and Size Information , 1995, NeuroImage.

[21]  K. Worsley.  Estimating the number of peaks in a random field using the Hadwiger characteristic of excursion sets, with applications to medical images , 1995 .

[22]  Karl J. Friston,et al.  Non-linear regression in parametric activation studies , 1996, NeuroImage.

[23]  Karl J. Friston,et al.  Nonlinear Regression in Parametric Activation Studies , 1996, NeuroImage.

[24]  Keith J. Worsley,et al.  The Geometry of Random Images , 1996 .

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

[26]  Karl J. Friston,et al.  Is Multivariate Analysis of PET Data More Revealing Than the Univariate Approach? Evidence from a Study of Episodic Memory Retrieval , 1996, NeuroImage.

[27]  Alan C. Evans,et al.  Searching scale space for activation in PET images , 1996, Human brain mapping.

[28]  C. Mountford.  Non-linear Fourier time series analysis for human brain mapping by functional magnetic resonance imaging - Discussion , 1997 .

[29]  Alan C. Evans,et al.  Scale space searches for a periodic signal in fMRI data with spatially varying hemodynamic response , 1997 .

[30]  Tony Lindeberg,et al.  Scale-Space Theory in Computer Vision , 1993, Lecture Notes in Computer Science.

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

[32]  Keith Worsley,et al.  The geometry of the Hotelling's random field with applications to the detection of shape changes , 1998 .

[33]  P. McMullen.  INTRODUCTION TO GEOMETRIC PROBABILITY , 1999 .

[34]  K. Worsley,et al.  The geometry of correlation fields with an application to functional connectivity of the brain , 1999 .

[35]  K. Worsley,et al.  THE DETECTION OF LOCAL SHAPE CHANGES VIA THE GEOMETRY OF HOTELLING’S T 2 FIELDS 1 , 1999 .

[36]  R. Adler.  On excursion sets, tube formulas and maxima of random fields , 2000 .

[37]  Emery N. Brown,et al.  A signal estimation approach to functional MRI , 2001, IEEE Transactions on Medical Imaging.