Finding electrophysiological sources of aging-related processes using penalized least squares with Modified Newton-Raphson algorithm

In this work, we evaluate the flexibility of a modified Newton-Raphson (MNR) algorithm for finding electrophysiological sources in both simulated and real data, and then apply it to different penalized models in order to compare the sources of the EEG theta rhythm in two groups of elderly subjects with different levels of declined physical performance. As a first goal, we propose the MNR algorithm for estimating general multiple penalized least squares (MPLS) models and show that it is capable to find solutions that are simultaneously sparse and smooth. This algorithm allowed to address known and novel models such as the Smooth Non-negative Garrote and the Non-negative Smooth LASSO. We test its ability to solve the EEG inverse problem with multiple penalties -using simulated data- in terms of localization error, blurring and visibility, as compared with traditional algorithms. As a second goal, we explore the electrophysiological sources of the theta activity extracted from resting-state EEG recorded in two groups of older adults, which belong to a longitudinal study to assess the relationship between measures of physical performance (gait speed) decline and normal cognition. The groups contained subjects with good and bad physical performance in the two evaluations (6 years apart). In accordance to clinical studies, we found differences in EEG theta sources for the two groups, specifically, subjects with declined physical performance presented decreased temporal sources while increased prefrontal sources that seem to reflect compensating mechanisms to ensure a stable walking

Revista de logopedia, foniatría y audiología

Resumen en inglés. Resumen tomado de la publicaciónSe realiza un estudio de los diferentes métodos de cribado de pérdidas auditivas en recién nacidos que se han difundido en los últimos años, con el propósito de lograr averiguar cuál de ellos logra mayor precisión a la hora de detectar los problemas auditivos. En el proceso se utilizaron diferentes materiales y métodos como: muestras, estímulos, calibraciones acústicas, procedimientos de registro e indicadores estadísticos que serán los que proporcionen la información sobre cuál es al final el método más preciso para el cribado temprano de pérdidas auditivas.CataluñaConsejería de Educación. Dirección General de Política Educativa; Calle Delgado Valencia, 6; 06800 Mérida (Badajoz); Tel. +34924006714; Fax +34924006716; [email protected]

Spatio Temporal EEG Source Imaging with the Hierarchical Bayesian Elastic Net and Elitist Lasso Models

The estimation of EEG generating sources constitutes an Inverse Problem (IP) in Neuroscience. This is an ill-posed problem due to the non-uniqueness of the solution and regularization or prior information is needed to undertake Electrophysiology Source Imaging. Structured Sparsity priors can be attained through combinations of (L1 norm-based) and (L2 norm-based) constraints such as the Elastic Net (ENET) and Elitist Lasso (ELASSO) models. The former model is used to find solutions with a small number of smooth nonzero patches, while the latter imposes different degrees of sparsity simultaneously along different dimensions of the spatio-temporal matrix solutions. Both models have been addressed within the penalized regression approach, where the regularization parameters are selected heuristically, leading usually to non-optimal and computationally expensive solutions. The existing Bayesian formulation of ENET allows hyperparameter learning, but using the computationally intensive Monte Carlo/Expectation Maximization methods, which makes impractical its application to the EEG IP. While the ELASSO have not been considered before into the Bayesian context. In this work, we attempt to solve the EEG IP using a Bayesian framework for ENET and ELASSO models. We propose a Structured Sparse Bayesian Learning algorithm based on combining the Empirical Bayes and the iterative coordinate descent procedures to estimate both the parameters and hyperparameters. Using realistic simulations and avoiding the inverse crime we illustrate that our methods are able to recover complicated source setups more accurately and with a more robust estimation of the hyperparameters and behavior under different sparsity scenarios than classical LORETA, ENET and LASSO Fusion solutions. We also solve the EEG IP using data from a visual attention experiment, finding more interpretable neurophysiological patterns with our methods. The Matlab codes used in this work, including Simulations, Methods, Quality Measures and Visualization Routines are freely available in a public website

Estimating brain functional connectivity with sparse multivariate autoregression

There is much current interest in identifying the anatomical and functional circuits that are the basis of the brain's computations, with hope that functional neuroimaging techniques will allow the in vivo study of these neural processes through the statistical analysis of the time-series they produce. Ideally, the use of techniques such as multivariate autoregressive (MAR) modelling should allow the identification of effective connectivity by combining graphical modelling methods with the concept of Granger causality. Unfortunately, current time-series methods perform well only for the case that the length of the time-series Nt is much larger than p, the number of brain sites studied, which is exactly the reverse of the situation in neuroimaging for which relatively short time-series are measured over thousands of voxels. Methods are introduced for dealing with this situation by using sparse MAR models. These can be estimated in a two-stage process involving (i) penalized regression and (ii) pruning of unlikely connections by means of the local false discovery rate developed by Efron. Extensive simulations were performed with idealized cortical networks having small world topologies and stable dynamics. These show that the detection efficiency of connections of the proposed procedure is quite high. Application of the method to real data was illustrated by the identification of neural circuitry related to emotional processing as measured by BOLD