Spatio-temporal correlations from fMRI time series based on the NN-ARx model.

For the purpose of statistical characterization of the spatio-temporal correlation structure of brain functioning from high-dimensional fMRI time series, we introduce an innovation approach. This is based on whitening the data by the Nearest-Neighbors AutoRegressive model with external inputs (NN-ARx). Correlations between the resulting innovations are an extension of the usual correlations, in which mean-correction is carried out by the dynamic NN-ARx model instead of the static, standard linear model for fMRI time series. Measures of dependencies between regions are defined by summarizing correlations among innovations at several time lags over pairs of voxels. Such summarization does not involve averaging the data over each region, which prevents loss of information in case of non-homogeneous regions. Statistical tests based on these measures are elaborated, which allow for assessing the correlation structure in search of connectivity. Results of application of the NN-ARx approach to fMRI data recorded in visual stimuli experiments are shown. Finally, a number of issues related with its potential and limitations are commented.

3D statistical parametric mapping of EEG source spectra by means of variable resolution electromagnetic tomography (VARETA).

This article describes a new method for 3D QEEG tomography in the frequency domain. A variant of Statistical Parametric Mapping is presented for source log spectra. Sources are estimated by means of a Discrete Spline EEG inverse solution known as Variable Resolution Electromagnetic Tomography (VARETA). Anatomical constraints are incorporated by the use of the Montreal Neurological Institute (MNI) probabilistic brain atlas. Efficient methods are developed for frequency domain VARETA in order to estimate the source spectra for the set of 10(3)-10(5) voxels that comprise an EEG/MEG inverse solution. High resolution source Z spectra are then defined with respect to the age dependent mean and standard deviations of each voxel, which are summarized as regression equations calculated from the Cuban EEG normative database. The statistical issues involved are addressed by the use of extreme value statistics. Examples are shown that illustrate the potential clinical utility of the methods herein developed.

The estimation of event related potentials affected by random shifts and scalings.

There is considerable evidence for trial to trial variability of the event related potentials (ERPs) within a given subject's recording. This variability influences the outcome of usual procedures in ERP analysis. Better results may be obtained if the sources of variability are explicitly taken into account in an appropriate model. This paper considers a probabilistic model, the random shift and scaling (RSS) model, where the response is modified by a random time shift and a random scale factor. In addition to this, an additional random scale factor which affects both the response and the background noise is taken into account. This time shift and these scale factors are handled as nuisance parameters. Maximum likelihood and least squares estimators of these parameters and the waveform of response are derived for the RSS model. It is shown that the Woody estimate of the ERP reported in earlier work can be derived by restricting the assumptions for the RSS model. Test statistics for hypotheses on means are obtained for the RSS model and a new type of discriminant function. The usefulness of the method is illustrated by means of simulation studies. Receiver operating characteristic (ROC) curves are used to demonstrate that the new type of discriminant performs better than the usual Fisher's Linear Discriminant.

Model parameters estimation when the evoked potential recordings are affected by a random scale factor.

In many situations an important source of the average evoked potentials (EPs) variability is a random scale factor affecting each recording. As a result, the outcome of any EP detection method may be greatly affected. However, using an appropriate probabilistic model these scale factor can be estimated, and the performance of any available detection index improved by data rescaling. In this paper the Maximum Likelihood Estimators of the waveform of the response and the scale factor affecting both background noise and this waveform are obtained. Also, an iterative algorithm for model parameters estimation is presented and its convergence is examined in a simulation study. The Linear Discriminant function is computed using simulated test data in both situations, before and after rescaling of recordings. The performance of these statistics is evaluated by mean of ROC curves.