-
Sparse electrophysiological source imaging predicts aging-related gait speed slowing
Authors:
Vega-Hernandez,
Mayrim,
Galan-Garcia,
Lidice,
Perez-Hidalgo-Gato,
Jhoanna,
Ontivero-Ortega,
Marlis,
Garcia-Agustin,
Daysi,
Garcia-Reyes,
Ronaldo,
Bosch-Bayard,
Jorge,
Marinazzo,
Daniele,
Martinez-Montes,
Eduardo,
Valdes-Sosa,
Pedro A
Abstract:
Objective: We seek stable Electrophysiological Source Imaging (ESI) biomarkers associated with Gait Speed (GS) as a measure of functional decline. Towards this end we determine the predictive value of ESI activation and connectivity patterns of resting-state EEG Theta rhythm on physical performance decline measured by a slowing GS in aging individuals. Methods: As potential biomarkers related to G…
▽ More
Objective: We seek stable Electrophysiological Source Imaging (ESI) biomarkers associated with Gait Speed (GS) as a measure of functional decline. Towards this end we determine the predictive value of ESI activation and connectivity patterns of resting-state EEG Theta rhythm on physical performance decline measured by a slowing GS in aging individuals. Methods: As potential biomarkers related to GS changes, we estimate ESI using flexible sparse/smooth/non-negative models (NN-SLASSO), from which activation ESI (aESI) and connectivity ESI (cESI) features are selected using the Stable Sparse Classifier method. Results and Conclusions: Novel sparse aESI models outperformed traditional methods such as the LORETA family. The models combining aESI and cESI features improved the predictability of GS changes. Selected biomarkers from activation/connectivity patterns were localized to orbitofrontal and temporal cortical regions. Significance: The proposed methodology contributes to understanding the activation and connectivity of ESI complex patterns related to GS, providing potential biomarker features for GS slowing. Given the known relationship between GS decline and cognitive impairment, this preliminary work suggests it might be applied to other, more complex measures of healthy and pathological aging. Importantly, it might allow an ESI-based evaluation of rehabilitation programs.
△ Less
Submitted 22 August, 2023; v1 submitted 20 July, 2023;
originally announced July 2023.
-
A Multi-objective Evolutionary Algorithm for EEG Inverse Problem
Authors:
José Enrique Alvarez Iglesias,
Mayrim Vega-Hernández,
Eduardo Martínez-Montes
Abstract:
In this paper, we proposed a multi-objective approach for the EEG Inverse Problem. This formulation does not need unknown parameters that involve empirical procedures. Due to the combinatorial characteristics of the problem, this alternative included evolutionary strategies to resolve it. The result is a Multi-objective Evolutionary Algorithm based on Anatomical Restrictions (MOEAAR) to estimate d…
▽ More
In this paper, we proposed a multi-objective approach for the EEG Inverse Problem. This formulation does not need unknown parameters that involve empirical procedures. Due to the combinatorial characteristics of the problem, this alternative included evolutionary strategies to resolve it. The result is a Multi-objective Evolutionary Algorithm based on Anatomical Restrictions (MOEAAR) to estimate distributed solutions. The comparative tests were between this approach and 3 classic methods of regularization: LASSO, Ridge-L and ENET-L. In the experimental phase, regression models were selected to obtain sparse and distributed solutions. The analysis involved simulated data with different signal-to-noise ratio (SNR). The indicators for quality control were Localization Error, Spatial Resolution and Visibility. The MOEAAR evidenced better stability than the classic methods in the reconstruction and localization of the maximum activation. The norm L0 was used to estimate sparse solutions with the evolutionary approach and its results were relevant.
△ Less
Submitted 26 December, 2021; v1 submitted 21 July, 2021;
originally announced July 2021.
-
Multiple penalized least squares and sign constraints with modified Newton-Raphson algorithms: application to EEG source imaging
Authors:
Mayrim Vega-Hernández,
José M. Sánchez-Bornot,
Agustín Lage-Castellanos,
Jhoanna Pérez-Hidalgo-Gato,
Darío Palmero-Ledón,
José E. Alvarez-Iglesias,
Daysi García-Agustin,
Eduardo Martínez-Montes,
Pedro A. Valdés-Sosa
Abstract:
Multiple penalized least squares (MPLS) models are a flexible approach to find adaptive least squares solutions required to be simultaneously sparse and smooth. This is particularly important when addressing real-life inverse problems where there is no ground truth available, such as electrophysiological source imaging. In this work we formalize a modified Newton-Raphson (MNR) algorithm to estimat…
▽ More
Multiple penalized least squares (MPLS) models are a flexible approach to find adaptive least squares solutions required to be simultaneously sparse and smooth. This is particularly important when addressing real-life inverse problems where there is no ground truth available, such as electrophysiological source imaging. In this work we formalize a modified Newton-Raphson (MNR) algorithm to estimate general MPLS models and propose its extension to perform efficient optimization over the active set of selected features (AMNR). This algorithm can be used to minimize continuously differentiable objective functions with multiple restrictions, including sign constraints. We show that these algorithms provide solutions with acceptable reconstruction in simulated scenarios that do not cope with model assumptions, and for low n/p ratios. We then use both algorithms for estimating different electroencephalography (EEG) inverse models with multiple penalties. We also show how the AMNR allows us to estimate new models in the EEG inverse problem context, such as nonnegative versions of Smooth Garrote and Smooth LASSO. Synthetic data were used for a comparison between solutions obtained using the two algorithms proposed here and the least angle regression (LARS) algorithm, according to well-known quality measures. A visual event-related EEG from healthy young subjects and a resting-state EEG study on the relationship between cognitive ageing and walking speed decline in active elders, were used to illustrate the usefulness of the proposed methods in the analysis of real experimental data.
△ Less
Submitted 24 August, 2021; v1 submitted 5 November, 2019;
originally announced November 2019.
-
Modeling functional resting-state brain networks through neural message passing on the human connectome
Authors:
Julio A. Peraza-Goicolea,
Eduardo Martínez-Montes,
Eduardo Aubert,
Pedro A. Valdés-Hernández,
Roberto Mulet
Abstract:
Understanding the relationship between the structure and function of the human brain is one of the most important open questions in Neurosciences. In particular, Resting State Networks (RSN) and more specifically the Default Mode Network (DMN) of the brain, which are defined from the analysis of functional data lack a definitive justification consistent with the anatomical structure of the brain.…
▽ More
Understanding the relationship between the structure and function of the human brain is one of the most important open questions in Neurosciences. In particular, Resting State Networks (RSN) and more specifically the Default Mode Network (DMN) of the brain, which are defined from the analysis of functional data lack a definitive justification consistent with the anatomical structure of the brain. In this work, we show that a possible connection may naturally rest on the idea that information flows in the brain through a neural message-passing dynamics between macroscopic structures, like those defined by the human connectome (HC). In our model, each brain region in the HC is assumed to have a binary behavior (active or not), the strength of interactions among them is encoded in the anatomical connectivity matrix defined by the HC, and the dynamics of the system is defined by a neural message-passing algorithm, Belief Propagation (BP), working near the critical point of the human connectome. We show that in the absence of direct external stimuli the BP algorithm converges to a spatial map of activations that is similar to the DMN. Moreover, we computed, using Susceptibility Propagation (SP), the matrix of correlations between the different regions and show that the modules defined by a clustering of this matrix resemble several Resting States Networks determined experimentally. Both results suggest that the functional DMN and RSNs can be seen as simple consequences of the anatomical structure of the brain and a neural message-passing dynamics between macroscopic regions. We then show preliminary results indicating our predictions on how functional DMN maps change when the anatomical brain network suffers structural anomalies, like in Alzheimers Disease and in lesions of the Corpus Callosum.
△ Less
Submitted 12 June, 2019;
originally announced June 2019.
-
Identification of oscillatory brain networks with Hidden Gaussian Graphical Spectral models of EEG/MEG
Authors:
Deirel Paz-Linares,
Eduardo Gonzalez-Moreira,
Ariosky Areces-Gonzalez,
Ying Wang,
Min Li,
Eduardo Martinez-Montes,
Jorge Bosch-Bayard,
Maria L. Bringas-Vega,
Mitchel J. Valdes-Sosa,
Pedro A. Valdes-Sosa
Abstract:
Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflect functional brain network activity with a temporal resolution that no other in vi…
▽ More
Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflect functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our inverse-solution with 1/3 more congruence according to Riemannian distances than state-of-the-art methods.
△ Less
Submitted 25 October, 2022; v1 submitted 2 October, 2018;
originally announced October 2018.
-
Empirical bayes formulation of the elastic net and mixed-norm models: application to the eeg inverse problem
Authors:
Deirel Paz-Linares,
Mayrim Vega-Hernández,
Pedro A. Rojas-López,
Pedro A. Valdés-Sosa,
Eduardo Martínez-Montes
Abstract:
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 many kinds of prior information have been used to constrain it. A combination of smoothness (L2 norm-based) and sparseness (L1 norm-based) constraints is a flexible approach that have been pursued by important examples such as th…
▽ More
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 many kinds of prior information have been used to constrain it. A combination of smoothness (L2 norm-based) and sparseness (L1 norm-based) constraints is a flexible approach that have been pursued by important examples such as the Elastic Net (ENET) and mixed-norm (MXN) models. The former is used to find solutions with a small number of smooth non-zero patches, while the latter imposes sparseness and smoothness 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 solutions. The existing Bayesian formulation of ENET allows hyperparameter learning, but using computationally intensive Monte Carlo/Expectation Maximization methods. In this work we attempt to solve the EEG IP using a Bayesian framework for models based on mixtures of L1/L2 norms penalization functions (Laplace/Normal priors) such as ENET and MXN. We propose a Sparse Bayesian Learning algorithm based on combining the Empirical Bayes and the iterative coordinate descent procedures to estimate both the parameters and hyperparameters. Using simple but realistic simulations we found that our methods are able to recover complicated source setups more accurately and with a more robust variable selection than the ENET and LASSO solutions using classical algorithms. We also solve the EEG IP using data coming from a visual attention experiment, finding more interpretable neurophysiological patterns with our methods, as compared with other known methods such as LORETA, ENET and LASSO FUSION using the classical regularization approach.
△ Less
Submitted 6 February, 2016; v1 submitted 25 January, 2016;
originally announced January 2016.