Source author record

Roberto D. Pascual-Marqui

Roberto D. Pascual-Marqui appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2024arXiv

Lagged coherence: explicit and testable definition

Measures of association between cortical regions based on activity signals provide useful information for studying brain functional connectivity. Difficulties occur with signals of electric neuronal activity, where an observed signal is a mixture, i.e. an instantaneous weighted average of the true, unobserved signals from all regions, due to volume conduction and low spatial resolution. This is why measures of lagged association are of interest, since at least theoretically, "lagged association" is of physiological origin. In contrast, the actual physiological instantaneous zero-lag association is masked and confounded by the mixing artifact. A minimum requirement for a measure of lagged association is that it must not tend to zero with an increase of strength of true instantaneous physiological association. Such biased measures cannot tell apart if a change in its value is due to a change in lagged or a change in instantaneous association. An explicit testable definition for frequency domain lagged connectivity between two multivariate time series is proposed. It is endowed with two important properties: it is invariant to non-singular linear transformations of each vector time series separately, and it is invariant to instantaneous association. As a first sanity check: in the case of two univariate time series, the new definition leads back to the bivariate lagged coherence of 2007 (eqs 25 and 26 in https://doi.org/10.48550/arXiv.0706.1776). As a second stronger sanity check: in the case of a univariate and multivariate vector time series, the new measure presented here leads back to the original multivariate lagged coherence of 2007 (eq 31 in https://doi.org/10.48550/arXiv.0711.1455), which again trivially includes the bivariate case.

preprint2022arXiv

Cortical Xi-Alpha model for resting state electric neuronal activity

Power spectra of awake resting state EEG recordings in humans typically have an Alpha peak at around 10 Hz riding a decreasing background "Xi process". Normal and pathological variations may have more than one peak or none. The single channel Xi-Alpha model (Pascual-Marqui et al 1988, https://doi.org/10.3109/00207458808985730) separated these two processes, providing a low-dimensional parametric description of EEG spectra. Currently lacking is a generative whole cortex model for activity spectra and intracortical functional connectivity. Here we introduce the "cortical Xi-Alpha model". The cross-spectral density matrices are modeled as additive components, each one consisting of a scalar spectrum that multiplies a frequency invariant Hermitian covariance matrix. This simple "separation of variables" form offers a very rich repertoire of spatio-spectral properties, as well as diverse whole cortex functional connectivity patterns. The scalp EEG model conserves the same form, allowing simple estimation from scalp to cortex. Two independent open-access, resting state eyes open and closed EEG data sets (203 participants with 61 electrodes, and 47 participants with 26 electrodes) were used to demonstrate, test, and validate the model. Results summary: - The average dimension of cross-spectra lies between 1.7 and 2.6, indicating two processes. - Non-negative matrix factorization of population power spectra sampled at 6239 cortical voxels with only two components (identified as Xi and Alpha) explains 99% of the variance. - The median value of explained variance was 95% for the "cortical Xi-Alpha model" across all datasets and conditions. - Alpha process generators more occipital, Xi more frontal. - Xi cortical lagged connectivities are isotropic with interdistance. - Laminar recordings suggest layers 2/3 pyramidal neurons generate Xi; layers 5/6 pyramidal neurons generate Alpha.

preprint2022arXiv

On the relation between EEG microstates and cross-spectra

Brain function as measured by multichannel EEG recordings can be described to a high level of accuracy by microstates, characterized as a sequence of time intervals within which the sign invariant normalized scalp electric potential field remains quasi-stable, concatenated by fast transitions. Filtering the EEG has a small effect on the spatial microstate scalp maps, but a large effect on the dynamics (e.g. duration, frequency of occurrence, and transition rates). In addition, spectral power has been found to be strongly correlated with microstate dynamics. And yet, the nature of the relation between spectra and microstates remains poorly understood. Here we show that the multivariate EEG cross-spectrum contains sufficient generative information for estimating the microstate scalp maps and their dynamics, demonstrating an underlying fundamental link between the microstate model and the multivariate cross-spectrum. Empirically, based on EEG recordings from 203 participants in eyes-closed resting state, their cross-spectral matrices were computed, from which stochastic EEG was generated. No significant differences were found for the microstate model (maps and dynamics) estimated from the actual EEG and from the stochastic EEG based solely on the cross-spectra. In addition, with the aim of quantifying the spatio-cross-spectral properties of the microstate model, we introduce here the topographic likelihood spectrum, based on the Watson distribution, which provides a frequency-by-frequency account of the contribution of a normalized microstate map to the normalized EEG cross-spectrum, independent of power. The topographic likelihood spectra are distinct for the different microstate maps. In a comparison between eyes-closed and eyes open conditions, they are shown to be significantly different in frequency specific patterns.

preprint2013arXiv

A measure of association between vectors based on "similarity covariance"

The "maximum similarity correlation" definition introduced in this study is motivated by the seminal work of Szekely et al on "distance covariance" (Ann. Statist. 2007, 35: 2769-2794; Ann. Appl. Stat. 2009, 3: 1236-1265). Instead of using Euclidean distances "d" as in Szekely et al, we use "similarity", which can be defined as "exp(-d/s)", where the scaling parameter s>0 controls how rapidly the similarity falls off with distance. Scale parameters are chosen by maximizing the similarity correlation. The motivation for using "similarity" originates in spectral clustering theory (see e.g. Ng et al 2001, Advances in Neural Information Processing Systems 14: 849-856). We show that a particular form of similarity correlation is asymptotically equivalent to distance correlation for large values of the scale parameter. Furthermore, we extend similarity correlation to coherence between complex valued vectors, including its partitioning into real and imaginary contributions. Several toy examples are used for comparing distance and similarity correlations. For instance, points on a noiseless straight line give distance and similarity correlation values equal to 1; but points on a noiseless circle produces near zero distance correlation (dCorr=0.02) while the similarity correlation is distinctly non zero (sCorr=0.36). In distinction to the distance approach, similarity gives more importance to small distances, which emphasizes the local properties of functional relations. This paper represents a preliminary empirical study, showing that the novel similarity association has some distinct practical advantages over distance based association.For the sake of reproducible research, the software code implementing all methods here (using lazarus free-pascal "www.lazarus.freepascal.org"), including all test data, are freely available at: "sites.google.com/site/pascualmarqui/home/similaritycovariance".

preprint2011arXiv

Cortical current source connectivity by means of partial coherence fields

An important field of research in functional neuroimaging is the discovery of integrated, distributed brain systems and networks, whose different regions need to work in unison for normal functioning. The EEG is a non-invasive technique that can provide information for massive connectivity analyses. Cortical signals of time varying electric neuronal activity can be estimated from the EEG. Although such techniques have very high time resolution, two cortical signals even at distant locations will appear to be highly similar due to the low spatial resolution nature of the EEG. In this study a method for eliminating the effect of common sources due to low spatial resolution is presented. It is based on an efficient estimation of the whole-cortex partial coherence matrix. Using as a starting point any linear EEG tomography that satisfies the EEG forward equation, it is shown that the generalized partial coherences for the cortical grey matter current density time series are invariant to the selected tomography. It is empirically shown with simulation experiments that the generalized partial coherences have higher spatial resolution than the classical coherences. The results demonstrate that with as little as 19 electrodes, lag-connected brain regions can often be missed and misplaced even with lagged coherence measures, while the new method detects and localizes correctly the connected regions using the lagged partial coherences.

preprint2011arXiv

Interaction patterns of brain activity across space, time and frequency. Part I: methods

We consider exploratory methods for the discovery of cortical functional connectivity. Typically, data for the i-th subject (i=1...NS) is represented as an NVxNT matrix Xi, corresponding to brain activity sampled at NT moments in time from NV cortical voxels. A widely used method of analysis first concatenates all subjects along the temporal dimension, and then performs an independent component analysis (ICA) for estimating the common cortical patterns of functional connectivity. There exist many other interesting variations of this technique, as reviewed in [Calhoun et al. 2009 Neuroimage 45: S163-172]. We present methods for the more general problem of discovering functional connectivity occurring at all possible time lags. For this purpose, brain activity is viewed as a function of space and time, which allows the use of the relatively new techniques of functional data analysis [Ramsay & Silverman 2005: Functional data analysis. New York: Springer]. In essence, our method first vectorizes the data from each subject, which constitutes the natural discrete representation of a function of several variables, followed by concatenation of all subjects. The singular value decomposition (SVD), as well as the ICA of this new matrix of dimension [rows=(NT*NV); columns=NS] will reveal spatio-temporal patterns of connectivity. As a further example, in the case of EEG neuroimaging, Xi of size NVxNW may represent spectral density for electric neuronal activity at NW discrete frequencies from NV cortical voxels, from the i-th EEG epoch. In this case our functional data analysis approach would reveal coupling of brain regions at possibly different frequencies.

preprint2010arXiv

Dynamic interactions in terms of senders, hubs, and receivers (SHR) using the singular value decomposition of time series: Theory and brain connectivity applications

Understanding of normal and pathological brain function requires the identification and localization of functional connections between specialized regions. The availability of high time resolution signals of electric neuronal activity at several regions offers information for quantifying the connections in terms of information flow. When the signals cover the whole cortex, the number of connections is very large, making visualization and interpretation very difficult. We introduce here the singular value decomposition of time-lagged multiple signals, which localizes the senders, hubs, and receivers (SHR) of information transmission. Unlike methods that operate on large connectivity matrices, such as correlation thresholding and graph-theoretic analyses, this method operates on the multiple time series directly, providing 3D brain images that assign a score to each location in terms of its sending, relaying, and receiving capacity. The scope of the method is general and encompasses other applications outside the field of brain connectivity.