1
|
Bullmore E, Fadili J, Breakspear M, Salvador R, Suckling J, Brammer M. Wavelets and statistical analysis of functional magnetic resonance images of the human brain. Stat Methods Med Res 2016; 12:375-99. [PMID: 14599002 DOI: 10.1191/0962280203sm339ra] [Citation(s) in RCA: 99] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Wavelets provide an orthonormal basis for multiresolution analysis and decorrelation or ‘whitening’ of nonstationary time series and spatial processes. Wavelets are particularly well suited to analysis of biological signals and images, such as human brain imaging data, which often have fractal or scale-invariant properties. We briefly define some key properties of the discrete wavelet transform (DWT) and review its applications to statistical analysis of functional magnetic resonance imaging (fMRI) data. We focus on time series resampling by ‘wavestrapping’ of wavelet coefficients, methods for efficient linear model estimation in the wavelet domain, and wavelet-based methods for multiple hypothesis testing, all of which are somewhat simplified by the decorrelating property of the DWT.
Collapse
Affiliation(s)
- Ed Bullmore
- Brain Mapping Unit and Wolfson Brain Imaging Centre, University of Cambridge, Addenbrooke's Hospital, Cambridge, UK
| | | | | | | | | | | |
Collapse
|
2
|
Kapur K, Roy A, Bhaumik DK, Gibbons RD, Lazar NA, Sweeney JA, Aryal S, Patterson D. Estimation and Classification of BOLD Responses Over Multiple Trials. COMMUN STAT-THEOR M 2009; 38:3099-3113. [DOI: 10.1080/03610920902947576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
3
|
Neumann J, von Cramon DY, Lohmann G. Model-based clustering of meta-analytic functional imaging data. Hum Brain Mapp 2008; 29:177-92. [PMID: 17390315 PMCID: PMC2885605 DOI: 10.1002/hbm.20380] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
We present a method for the analysis of meta-analytic functional imaging data. It is based on Activation Likelihood Estimation (ALE) and subsequent model-based clustering using Gaussian mixture models, expectation-maximization (EM) for model fitting, and the Bayesian Information Criterion (BIC) for model selection. Our method facilitates the clustering of activation maxima from previously performed imaging experiments in a hierarchical fashion. Regions with a high concentration of activation coordinates are first identified using ALE. Activation coordinates within these regions are then subjected to model-based clustering for a more detailed cluster analysis. We demonstrate the usefulness of the method in a meta-analysis of 26 fMRI studies investigating the well-known Stroop paradigm.
Collapse
Affiliation(s)
- Jane Neumann
- Max-Planck-Institute for Human Cognitive and Brain Sciences, Stephanstrasse 1a, Leipzig, Germany.
| | | | | |
Collapse
|
4
|
Van De Ville D, Seghier ML, Lazeyras F, Blu T, Unser M. WSPM: Wavelet-based statistical parametric mapping. Neuroimage 2007; 37:1205-17. [PMID: 17689101 DOI: 10.1016/j.neuroimage.2007.06.011] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2007] [Revised: 05/23/2007] [Accepted: 06/03/2007] [Indexed: 10/23/2022] Open
Abstract
Recently, we have introduced an integrated framework that combines wavelet-based processing with statistical testing in the spatial domain. In this paper, we propose two important enhancements of the framework. First, we revisit the underlying paradigm; i.e., that the effect of the wavelet processing can be considered as an adaptive denoising step to "improve" the parameter map, followed by a statistical detection procedure that takes into account the non-linear processing of the data. With an appropriate modification of the framework, we show that it is possible to reduce the spatial bias of the method with respect to the best linear estimate, providing conservative results that are closer to the original data. Second, we propose an extension of our earlier technique that compensates for the lack of shift-invariance of the wavelet transform. We demonstrate experimentally that both enhancements have a positive effect on performance. In particular, we present a reproducibility study for multi-session data that compares WSPM against SPM with different amounts of smoothing. The full approach is available as a toolbox, named WSPM, for the SPM2 software; it takes advantage of multiple options and features of SPM such as the general linear model.
Collapse
Affiliation(s)
- Dimitri Van De Ville
- Biomedical Imaging Group, Ecole Polytechnique Fédérale de Lausanne (EPFL), and Department of Radiology and Medical Informatics, University Hospital Geneva, Switzerland.
| | | | | | | | | |
Collapse
|
5
|
Aston JAD, Turkheimer FE, Brett M. HBM functional imaging analysis contest data analysis in wavelet space. Hum Brain Mapp 2006; 27:372-9. [PMID: 16565952 PMCID: PMC6871402 DOI: 10.1002/hbm.20244] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
An analysis of the Functional Imaging Analysis Contest (FIAC) data is presented using spatial wavelet processing. This technique allows the image to be filtered adaptively according to the data itself, rather than relying on a predetermined filter. This adaptive filtering leads to better estimation of the parameters and contrasts in terms of mean squared error. It will be shown that by introducing a slight bias into the estimation, a large reduction in the variance can be achieved, leading to better overall mean squared error estimates. As no single filter needs to be preselected, results containing many scales of information can be found. In the FIAC data, it is shown that both small-scale and large-scale (smoother, more dispersed) effects occur. The combination of small- and large-scale effects detected in the FIAC data would be easy to miss using conventional single filter analysis.
Collapse
Affiliation(s)
- John A D Aston
- Institute of Statistical Science, Academia Sinica, Taiwan.
| | | | | |
Collapse
|
6
|
Neumann J, von Cramon DY, Forstmann BU, Zysset S, Lohmann G. The parcellation of cortical areas using replicator dynamics in fMRI. Neuroimage 2006; 32:208-19. [PMID: 16647272 DOI: 10.1016/j.neuroimage.2006.02.039] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2005] [Revised: 02/24/2006] [Accepted: 02/27/2006] [Indexed: 11/19/2022] Open
Abstract
In this paper, we show that replicator dynamics can be used as an exploratory analysis tool to detect subregions of cortical areas on the basis of the similarity between fMRI time series. As similarity measure, we propose to use canonical correlation, a multivariate extension to the typically employed Pearson's correlation coefficient. We applied the replicator process to data obtained from two different experimental paradigms in the search for subregions within the left lateral frontal cortex (LFC). In both cases, the replicator process resulted in a parcellation that corresponds to a recently suggested subdivision of the LFC in anterior-posterior direction. Most notably, these results were very consistent when compared across different measurements of a single subject and across a group of subjects.
Collapse
Affiliation(s)
- Jane Neumann
- Max-Planck-Institute for Human Cognitive and Brain Sciences, Stephanstrasse 1a, D-04103 Leipzig, Germany.
| | | | | | | | | |
Collapse
|
7
|
Affiliation(s)
- Dimitri Van De Ville
- Biomedical Imaging Group, Ecole Polytechnique Fédérale de Lausanne, Biomedical Imaging Group, Switzerland.
| | | | | |
Collapse
|
8
|
Anderson CM, Lowen SB, Renshaw PF. Emotional task-dependent low-frequency fluctuations and methylphenidate: Wavelet scaling analysis of 1/f-type fluctuations in fMRI of the cerebellar vermis. J Neurosci Methods 2006; 151:52-61. [PMID: 16427128 DOI: 10.1016/j.jneumeth.2005.09.020] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2004] [Revised: 09/07/2005] [Accepted: 09/07/2005] [Indexed: 10/25/2022]
Abstract
UNLABELLED Ion channel currents, neural firing patterns, and brain BOLD signals display 1/f-type fluctuations or fractal properties in time. By design, fMRI methods attempt to minimize the contribution of variance from low-frequency physiological 1/f-noise. New fMRI methods are described to visualize and measure 1/f-type BOLD fluctuations in volunteers recalling affectively neutral or emotional memories or meditating (i.e., attending to breathing) then retrospectively rating emotional content. A wavelet scaling exponent (alpha) was used to characterize signals from 0.015625 to 0.5Hz in cerebellar lobules VIII to X of the vermis (posterior inferior vermis; PIV), a region coordinating balance, eye tracking, locomotion, and vascular tone, and a possible site of pathology in attention deficit hyperactivity disorder (ADHD). RESULTS Changes in alpha and emotional measures were correlated in PIV voxels (r = 0.622, d.f .= 14, P < 0.0005), but not other regions examined. In contrast, conventional means and standard deviations of PIV voxels were unchanged. Methylphenidate, shown to decrease slow oscillations in rodent basal ganglia [Ruskin DN, Bergstrom DA, Shenker A, Freeman LE, Baek D, Walters JR. Drugs used in the treatment of attention-deficit/hyperactivity disorder affect postsynaptic firing rate and oscillation without preferential dopamine autoreceptor action. Biol Psychiatry 2001;49:340-50.], abolished task-dependent alpha changes in the PIV of an adult with ADHD. Wavelet analysis of long BOLD time series appears well suited to fractal physiology and studies of pharmacologically modulated cerebellar-thalamic-cortical function in ADHD or other psychiatric disorders.
Collapse
Affiliation(s)
- Carl M Anderson
- Department of Psychiatry, Harvard Medical School, Boston, MA, USA.
| | | | | |
Collapse
|
9
|
Hinrichs H, Scholz M, Noesselt T, Heinze HJ. Quantile estimation to derive optimized test thresholds for random field statistics. Neuroimage 2005; 27:116-29. [PMID: 15955713 DOI: 10.1016/j.neuroimage.2005.03.033] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2004] [Revised: 03/18/2005] [Accepted: 03/28/2005] [Indexed: 11/28/2022] Open
Abstract
We present a numerical method to estimate the true threshold values in random fields needed to determine the significance of apparent signals observed in noisy images. To accomplish this, a quantile estimation algorithm is applied to derive the threshold with a predefined confidence interval from a large number of simulated random fields. Also, a computationally efficient method for generating a random field simulation is presented using resampling techniques. Applying these techniques, thresholds have been determined for a large variety of parameter settings (smoothness, voxel size, brain shape, type of statistics). By means of interpolation techniques, thresholds for additional arbitrary settings can be quickly derived without the need to run individual simulations. Compared to the parametric approach of Worsley et al. (1996) (Worsley, K.J., Marrett, S., Neelin P., Vandal, A.C., Friston, K.J., Evans, A.C., 1996. A unified statistical approach for determining significant signals in images of cerebral activation. Hum. Brain Mapp. 4, 58-73) and Friston et al. (1991) (Friston, K.J., Frith, C.D., Liddle, P.F., Frackowiak, R.S. 1991. Comparing functional (PET) images: the assessment of significant change. J. Cereb. Blood Flow Metab. 11(4), 690-699), and to the Bonferroni approach, these optimized thresholds lead to higher levels of significance (i.e., lower p values) with a specific amount of activation especially with fields of moderate smoothness (i.e., with a relative full width half maximum between 2 and 6). Alternatively, the threshold for a specified level of significance can be lowered. This improved statistical sensitivity is illustrated by the analysis of an actual event related functional magnetic resonance data set, and its limitations are tested by determining the false positive rate with experimental MR noise data. The grid of estimated threshold values as well as the interpolation algorithm to derive thresholds for arbitrary parameter settings are made available over the internet (http://neuro2.med.uni-magdeburg.de/quantile_estimation).
Collapse
Affiliation(s)
- H Hinrichs
- Department of Neurology II, University of Magdeburg, Leipziger Street 44, D-39120 Magdeburg, Germany.
| | | | | | | |
Collapse
|
10
|
Bullmore E, Fadili J, Maxim V, Sendur L, Whitcher B, Suckling J, Brammer M, Breakspear M. Wavelets and functional magnetic resonance imaging of the human brain. Neuroimage 2005; 23 Suppl 1:S234-49. [PMID: 15501094 DOI: 10.1016/j.neuroimage.2004.07.012] [Citation(s) in RCA: 156] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2004] [Accepted: 07/01/2004] [Indexed: 02/08/2023] Open
Abstract
The discrete wavelet transform (DWT) is widely used for multiresolution analysis and decorrelation or "whitening" of nonstationary time series and spatial processes. Wavelets are naturally appropriate for analysis of biological data, such as functional magnetic resonance images of the human brain, which often demonstrate scale invariant or fractal properties. We provide a brief formal introduction to key properties of the DWT and review the growing literature on its application to fMRI. We focus on three applications in particular: (i) wavelet coefficient resampling or "wavestrapping" of 1-D time series, 2- to 3-D spatial maps and 4-D spatiotemporal processes; (ii) wavelet-based estimators for signal and noise parameters of time series regression models assuming the errors are fractional Gaussian noise (fGn); and (iii) wavelet shrinkage in frequentist and Bayesian frameworks to support multiresolution hypothesis testing on spatially extended statistic maps. We conclude that the wavelet domain is a rich source of new concepts and techniques to enhance the power of statistical analysis of human fMRI data.
Collapse
Affiliation(s)
- Ed Bullmore
- Brain Mapping Unit and Wolfson Brain Imaging Centre, University of Cambridge, Addenbrooke's Hospital, Cambridge, UK.
| | | | | | | | | | | | | | | |
Collapse
|
11
|
Aston JAD, Gunn RN, Hinz R, Turkheimer FE. Wavelet variance components in image space for spatiotemporal neuroimaging data. Neuroimage 2005; 25:159-68. [PMID: 15734352 DOI: 10.1016/j.neuroimage.2004.10.037] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2004] [Revised: 10/20/2004] [Accepted: 10/26/2004] [Indexed: 11/25/2022] Open
Abstract
Neuroimaging studies place great emphasis on not only the estimation but also the standard error estimates of underlying parameters derived from a temporal model. This allows inferences to be made about the signal estimates and resulting conclusions to be drawn about the underlying data. It can often be advantageous to interrogate temporal models after spatial transformation of the data into the wavelet domain. Wavelet bases provide a multiresolution decomposition of the spatial data dimension and an ensuing reduction in spatial correlation. However, widespread acceptance of these wavelet techniques has been hampered by the limited ability to reconstruct both parametric and error estimates into the image domain after analysis of temporal models in the wavelet domain. This paper introduces a derivation and a fast implementation of a method for the calculation of the variance of the parametric images obtained from wavelet filters. The technique is proposed for a class of estimators that have been shown to be useful in neuroimaging studies. The techniques are demonstrated for both functional magnetic resonance imaging (fMRI) and positron emission tomography (PET) data sets.
Collapse
Affiliation(s)
- John A D Aston
- Institute of Statistical Science, Academia Sinica, 128 Academia Road, Sec 2, Taipei 11529, Taiwan.
| | | | | | | |
Collapse
|
12
|
Van De Ville D, Blu T, Unser M. Integrated wavelet processing and spatial statistical testing of fMRI data. Neuroimage 2004; 23:1472-85. [PMID: 15589111 DOI: 10.1016/j.neuroimage.2004.07.056] [Citation(s) in RCA: 58] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2004] [Revised: 07/07/2004] [Accepted: 07/12/2004] [Indexed: 11/17/2022] Open
Abstract
We introduce an integrated framework for detecting brain activity from fMRI data, which is based on a spatial discrete wavelet transform. Unlike the standard wavelet-based approach for fMRI analysis, we apply the suitable statistical test procedure in the spatial domain. For a desired significance level, this scheme has one remaining degree of freedom, characterizing the wavelet processing, which is optimized according to the principle of minimal approximation error. This allows us to determine the threshold values in a way that does not depend on data. While developing our framework, we make only conservative assumptions. Consequently, the detection of activation is based on strong evidence. We have implemented this framework as a toolbox (WSPM) for the SPM2 software, taking advantage of multiple options and functions of SPM such as the setup of the linear model and the use of the hemodynamic response function. We show by experimental results that our method is able to detect activation patterns; the results are comparable to those obtained by SPM even though statistical assumptions are more conservative.
Collapse
Affiliation(s)
- Dimitri Van De Ville
- Biomedical Imaging Group, Swiss Federal Institute of Technology Lausanne (EPFL), CH-1015 Lausanne, Switzerland.
| | | | | |
Collapse
|
13
|
Shimizu Y, Barth M, Windischberger C, Moser E, Thurner S. Wavelet-based multifractal analysis of fMRI time series. Neuroimage 2004; 22:1195-202. [PMID: 15219591 DOI: 10.1016/j.neuroimage.2004.03.007] [Citation(s) in RCA: 79] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2003] [Revised: 02/27/2004] [Accepted: 03/01/2004] [Indexed: 12/28/2022] Open
Abstract
Functional magnetic resonance imaging (fMRI) time series are investigated with a multifractal method based on the Wavelet Modulus Maxima (WTMM) method to extract local singularity ("fractal") exponents. The spectrum of singularity exponents of each fMRI time series is quantified by spectral characteristics including its maximum and the corresponding dimension. We found that the range of Hölder exponents in voxels with activation is close to 1, whereas exponents are close to 0.5 in white matter voxels without activation. The maximum dimension decreases going from white matter to gray matter, and is lower still for activated time series. The full-width-at-half-maximum of the spectra is higher in activated areas. The proposed method becomes particularly effective when combining these spectral characteristics into a single parameter. Using these multifractal parameters, it is possible to identify activated areas in the human brain in both hybrid and in vivo fMRI data sets without knowledge of the stimulation paradigm applied.
Collapse
Affiliation(s)
- Yu Shimizu
- MR Centre of Excellence, Medical University of Vienna, Austria
| | | | | | | | | |
Collapse
|