Time-varying parametric modelling and time-dependent spectral characterisation with applications to EEG signals using multi-wavelets

A new time-varying autoregressive (TVAR) modelling approach is proposed for nonstationary signal processing and analysis, with application to EEG data modelling and power spectral estimation. In the new parametric modelling framework, the time-dependent coefficients of the TVAR model are represented using a novel multi-wavelet decomposition scheme. The time-varying modelling problem is then reduced to regression selection and parameter estimation, which can be effectively resolved by using a forward orthogonal regression algorithm. Two examples, one for an artificial signal and another for an EEG signal, are given to show the effectiveness and applicability of the new TVAR modelling method

Spectral Characterization of functional MRI data on voxel-resolution cortical graphs

The human cortical layer exhibits a convoluted morphology that is unique to each individual. Conventional volumetric fMRI processing schemes take for granted the rich information provided by the underlying anatomy. We present a method to study fMRI data on subject-specific cerebral hemisphere cortex (CHC) graphs, which encode the cortical morphology at the resolution of voxels in 3-D. We study graph spectral energy metrics associated to fMRI data of 100 subjects from the Human Connectome Project database, across seven tasks. Experimental results signify the strength of CHC graphs' Laplacian eigenvector bases in capturing subtle spatial patterns specific to different functional loads as well as experimental conditions within each task.Comment: Fixed two typos in the equations; (1) definition of L in section 2.1, paragraph 1. (2) signal de-meaning and normalization in section 2.4, paragraph

Activation detection in fMRI data via multi-scale singularity detection

Detection of active areas in the brain by functional magnetic resonance imaging (fMRI) is a challenging problem in medical imaging. Moreover, determining the onset and end of activation signal at specific locations in 3-space can determined networks of temporal relationships required for brain mapping. We introduce a method for activation detection in fMRI data via wavelet analysis of singular features. We pose the problem of determining activated areas as singularity detection in the temporal domain. Overcomplete wavelet expansion at integer scales are used to trace wavelet modulus maxima across scales to construct maxima lines. Form these maxima lines, singularities in the signal are located corresponding to the onset and end of an activation signal. We present result for simulated phantom waveforms and clinical fMRI dat from human finger tapping experiments. Different levels of noise were added to two waveforms of phantom data. No assumptions about specific frequency and amplitude of an activation signal were made prior to analysis. Detection was reliable for modest levels of random noise, but less precise at higher levels. For clinical fMRI data, activation maps were comparable to those of existing standard techniques

Automated analysis of quantitative image data using isomorphic functional mixed models, with application to proteomics data

Image data are increasingly encountered and are of growing importance in many areas of science. Much of these data are quantitative image data, which are characterized by intensities that represent some measurement of interest in the scanned images. The data typically consist of multiple images on the same domain and the goal of the research is to combine the quantitative information across images to make inference about populations or interventions. In this paper we present a unified analysis framework for the analysis of quantitative image data using a Bayesian functional mixed model approach. This framework is flexible enough to handle complex, irregular images with many local features, and can model the simultaneous effects of multiple factors on the image intensities and account for the correlation between images induced by the design. We introduce a general isomorphic modeling approach to fitting the functional mixed model, of which the wavelet-based functional mixed model is one special case. With suitable modeling choices, this approach leads to efficient calculations and can result in flexible modeling and adaptive smoothing of the salient features in the data. The proposed method has the following advantages: it can be run automatically, it produces inferential plots indicating which regions of the image are associated with each factor, it simultaneously considers the practical and statistical significance of findings, and it controls the false discovery rate.Comment: Published in at http://dx.doi.org/10.1214/10-AOAS407 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

ROC evaluation of statistical wavelet-based analysis of brain activation in [15O]-H2O PET scans

This paper presents and evaluates a wavelet-based statistical analysis of PET images for the detection of brain activation areas. Brain regions showing significant activations were obtained by performing Student's t tests in the wavelet domain, reconstructing the final image from only those wavelet coefficients that passed the statistical test at a given significance level, and discarding artifacts introduced during the reconstruction process. Using Receiver Operating Characteristic (ROC) curves, we have compared this statistical analysis in the wavelet domain to the conventional image-domain Statistical Parametric Mapping (SPM) method. For obtaining an accurate assessment of sensitivity and specificity, we have simulated realistic single subject [15O]-H2O PET studies with different hyperactivation levels of the thalamic region. The results obtained from an ROC analysis show that the wavelet approach outperforms conventional SPM in identifying brain activation patterns. Using the wavelet method, activation areas detected were closer in size and shape to the region actually activated in the reference image.Publicad

Diffusion-adapted spatial filtering of fMRI data for improved activation mapping in white matter

Brain activation mapping using fMRI data has been mostly focused on finding detections in gray matter. Activations in white matter are harder to detect due to anatomical differences between both tissue types, which are rarely acknowledged in experimental design. However, recent publications have started to show evidence for the possibility of detecting meaningful activations in white matter. The shape of the activations arising from the BOLD signal is fundamentally different between white matter and gray matter, a fact which is not taken into account when applying isotropic Gaussian filtering in the preprocessing of fMRI data. We explore a graph-based description of the white matter developed from diffusion MRI data, which is capable of encoding the anisotropic domain. Based on this representation, two approaches to white matter filtering are tested, and their performance is evaluated on both semi-synthetic phantoms and real fMRI data. The first approach relies on heat kernel filtering in the graph spectral domain, and produced a clear increase in both sensitivity and specificity over isotropic Gaussian filtering. The second approach is based on spectral decomposition for the denosing of the signal, and showed increased specificity at the cost of a lower sensitivity.Novel approach to white matter filtering We introduced new advanced methods for filtering brain scans. Using them, we managed to improve the detection of activity in the white matter of the brain

Multiscale Granger causality analysis by \`a trous wavelet transform

Since interactions in neural systems occur across multiple temporal scales, it is likely that information flow will exhibit a multiscale structure, thus requiring a multiscale generalization of classical temporal precedence causality analysis like Granger's approach. However, the computation of multiscale measures of information dynamics is complicated by theoretical and practical issues such as filtering and undersampling: to overcome these problems, we propose a wavelet-based approach for multiscale Granger causality (GC) analysis, which is characterized by the following properties: (i) only the candidate driver variable is wavelet transformed (ii) the decomposition is performed using the \`a trous wavelet transform with cubic B-spline filter. We measure GC, at a given scale, by including the wavelet coefficients of the driver times series, at that scale, in the regression model of the target. To validate our method, we apply it to publicly available scalp EEG signals, and we find that the condition of closed eyes, at rest, is characterized by an enhanced GC among channels at slow scales w.r.t. eye open condition, whilst the standard Granger causality is not significantly different in the two conditions.Comment: 4 pages, 3 figure

Time course of information processing in visual and haptic object classification

Peer reviewedPublisher PD