Marc Lavielle | Search PubMed

PubMed2006

Are fMRI event-related response constant in time? A model selection answer.

An accurate estimation of the hemodynamic response function (HRF) in functional magnetic resonance imaging (fMRI) is crucial for a precise spatial and temporal estimate of the underlying neuronal processes. Recent works have proposed non-parametric estimation of the HRF under the hypotheses of linearity and stationarity in time. Biological literature suggests, however, that response magnitude may vary with attention or ongoing activity. We therefore test a more flexible model that allows for the variation of the magnitude of the HRF with time in a maximum likelihood framework. Under this model, the magnitude of the HRF evoked by a single event may vary across occurrences of the same type of event. This model is tested against a simpler model with a fixed magnitude using information theory. We develop a standard EM algorithm to identify the event magnitudes and the HRF. We test this hypothesis on a series of 32 regions (4 ROIS on eight subjects) of interest and find that the more flexible model is better than the usual model in most cases. The important implications for the analysis of fMRI time series for event-related neuroimaging experiments are discussed.

Algorithms↗

PubMed2005

A statistical approach for array CGH data analysis.

BACKGROUND: Microarray-CGH experiments are used to detect and map chromosomal imbalances, by hybridizing targets of genomic DNA from a test and a reference sample to sequences immobilized on a slide. These probes are genomic DNA sequences (BACs) that are mapped on the genome. The signal has a spatial coherence that can be handled by specific statistical tools. Segmentation methods seem to be a natural framework for this purpose. A CGH profile can be viewed as a succession of segments that represent homogeneous regions in the genome whose BACs share the same relative copy number on average. We model a CGH profile by a random Gaussian process whose distribution parameters are affected by abrupt changes at unknown coordinates. Two major problems arise: to determine which parameters are affected by the abrupt changes (the mean and the variance, or the mean only), and the selection of the number of segments in the profile. RESULTS: We demonstrate that existing methods for estimating the number of segments are not well adapted in the case of array CGH data, and we propose an adaptive criterion that detects previously mapped chromosomal aberrations. The performances of this method are discussed based on simulations and publicly available data sets. Then we discuss the choice of modeling for array CGH data and show that the model with a homogeneous variance is adapted to this context. CONCLUSIONS: Array CGH data analysis is an emerging field that needs appropriate statistical tools. Process segmentation and model selection provide a theoretical framework that allows precise biological interpretations. Adaptive methods for model selection give promising results concerning the estimation of the number of altered regions on the genome.

Algorithms↗