PubMed Health⌕ Search

Biomedical subjects

M J Fadili

Publications and source records attributed to M J Fadili.

4 recordsLinked to original sources

A comparative evaluation of wavelet-based methods for hypothesis testing of brain activation maps.

Wavelet-based methods for hypothesis testing are described and their potential for activation mapping of human functional magnetic resonance imaging (fMRI) data is investigated. In this approach, we emphasise convergence between methods of wavelet thresholding or shrinkage and the problem of hypothesis testing in both classical and Bayesian contexts. Specifically, our interest will be focused on the trade-off between type I probability error control and power dissipation, estimated by the area under the ROC curve. We describe a technique for controlling the false discovery rate at an arbitrary level of error in testing multiple wavelet coefficients generated by a 2D discrete wavelet transform (DWT) of spatial maps of fMRI time series statistics. We also describe and apply change-point detection with recursive hypothesis testing methods that can be used to define a threshold unique to each level and orientation of the 2D-DWT, and Bayesian methods, incorporating a formal model for the anticipated sparseness of wavelet coefficients representing the signal or true image. The sensitivity and type I error control of these algorithms are comparatively evaluated by analysis of "null" images (acquired with the subject at rest) and an experimental data set acquired from five normal volunteers during an event-related finger movement task. We show that all three wavelet-based algorithms have good type I error control (the FDR method being most conservative) and generate plausible brain activation maps (the Bayesian method being most powerful). We also generalise the formal connection between wavelet-based methods for simultaneous multiresolution denoising/hypothesis testing and methods based on monoresolution Gaussian smoothing followed by statistical testing of brain activation maps.

Algorithms↗

Wavelet-generalized least squares: a new BLU estimator of linear regression models with 1/f errors.

Long-memory noise is common to many areas of signal processing and can seriously confound estimation of linear regression model parameters and their standard errors. Classical autoregressive moving average (ARMA) methods can adequately address the problem of linear time invariant, short-memory errors but may be inefficient and/or insufficient to secure type 1 error control in the context of fractal or scale invariant noise with a more slowly decaying autocorrelation function. Here we introduce a novel method, called wavelet-generalized least squares (WLS), which is (to a good approximation) the best linear unbiased (BLU) estimator of regression model parameters in the context of long-memory errors. The method also provides maximum likelihood (ML) estimates of the Hurst exponent (which can be readily translated to the fractal dimension or spectral exponent) characterizing the correlational structure of the errors, and the error variance. The algorithm exploits the whitening or Karhunen-Loéve-type property of the discrete wavelet transform to diagonalize the covariance matrix of the errors generated by an iterative fitting procedure after both data and design matrix have been transformed to the wavelet domain. Properties of this estimator, including its Cramèr-Rao bounds, are derived theoretically and compared to its empirical performance on a range of simulated data. Compared to ordinary least squares and ARMA-based estimators, WLS is shown to be more efficient and to give excellent type 1 error control. The method is also applied to some real (neurophysiological) data acquired by functional magnetic resonance imaging (fMRI) of the human brain. We conclude that wavelet-generalized least squares may be a generally useful estimator of regression models in data complicated by long-memory or fractal noise.

Adult↗

On the number of clusters and the fuzziness index for unsupervised FCA application to BOLD fMRI time series.

The aim of this paper is to present an exploratory data-driven strategy based on Unsupervised Fuzzy Clustering Analysis (UFCA) and its potential for fMRI data analysis in the temporal domain. The a priori definition of the number of clusters is addressed and solved using heuristics. An original validity criterion is proposed taking into account data geometry and the partition Membership Functions (MFs). From our simulations, this criterion is shown to outperform other indices used in the literature. The influence of the fuzziness index was studied using simulated activation combined with real life noise data acquired from subjects under a resting state. Receiver Operating Characteristics (ROC) methodology is implemented to assess the performance of the proposed UFCA with respect to the fuzziness index. An interval of choice around 2, a value widely used in FCA, is shown to yield the best performance.

Algorithms↗

A multistep unsupervised fuzzy clustering analysis of fMRI time series.

A paradigm independent multistage strategy based on the Unsupervised Fuzzy Clustering Analysis (UFCA) and its potential for fMRI data analysis are presented. The influence of the fuzziness index is studied using Receiver Operating Characteristics (ROC) methodology and an interval of choice, around the widely used value 2, is shown to yield the best performance. The ill-balanced data problem is also overcome using a pre-processing step to reduce the number of voxels presented to the method. Statistical and anatomical criteria are proposed to exclude some voxels and enhance the UFCA sensitivity. An original postprocessing step aiming at statistically characterizing the obtained clusters is also developed. Two similarity criteria are used: the correlation coefficient on temporal profiles and a novel fuzzy overlap coefficient on membership degree maps. This final step provides a useful analysis tool to study intra-individual reproducibility of the classes across series (stimulation vs. stimulation, noise vs. noise or stimulation vs. noise). Finally, a comparison between this technique and some existing or locally developed postprocessing algorithms is presented using ROC methods. Its sensitivity and robustness is compared to the classical FCA or other techniques as a function of several parameters such as Contrast-to-Noise Ratio (CNR) and noise amplitude. Even without knowledge about the paradigm, the hemodynamic response function and the number of clusters, the performances of the proposed strategy are comparable to those of the classical approaches where extensive prior knowledge has to be added.

Algorithms↗