Cross-region dynamic connectivity, which describes spatio-temporal dependence of neural activity among multiple brain regions of interest (ROIs), can provide important information for understanding cognition. For estimating such connectivity, magnetoencephalography (MEG) and electroencephalography (EEG) are well-suited tools because of their millisecond temporal resolution. However, localizing source activity in the brain requires solving an under-determined linear problem. In typical two-step approaches, researchers first solve the linear problem with general priors assuming independence across ROIs, and secondly quantify cross-region connectivity. In this work, we propose a one-step state-space model to improve estimation of dynamic connectivity. The model treats the mean activity in individual ROIs as the state variable, and describes non-stationary dynamic dependence across ROIs using time-varying auto-regression. Compared with a two-step method, which first obtains the commonly used minimum-norm estimates of source activity, and then fits the auto-regressive model, our state-space model yielded smaller estimation errors on simulated data where the model assumptions held. When applied on empirical MEG data from one participant in a scene-processing experiment, our state-space model also demonstrated intriguing preliminary results, indicating leading and lagged linear dependence between the early visual cortex and a higher-level scene-sensitive region, which could reflect feed-forward and feedback information flow within the visual cortex during scene processing.

The anatomical connectivity of the human brain supports diverse patterns of correlated neural activity that are thought to underlie cognitive function. In a manner sensitive to underlying structural brain architecture, we examine the extent to which such patterns of correlated activity systematically vary across cognitive states. Anatomical white matter connectivity is compared with functional correlations in neural activity measured via blood oxygen level dependent (BOLD) signals. Functional connectivity is separately measured at rest, during an attention task, and during a memory task. We assess these structural and functional measures within previously-identified resting-state functional networks, denoted task-positive and task-negative networks, that have been independently shown to be strongly anticorrelated at rest but also involve regions of the brain that routinely increase and decrease in activity during task-driven processes. We find that the density of anatomical connections within and between task-positive and task-negative networks is differentially related to strong, task-dependent correlations in neural activity. The space mapped out by the observed structure-function relationships is used to define a quantitative measure of separation between resting, attention, and memory states. We find that the degree of separation between states is related to both general measures of behavioral performance and relative differences in task-specific measures of attention versus memory performance. These findings suggest that the observed separation between cognitive states reflects underlying organizational principles of human brain structure and function.

In order to construct spatial activation plots from functional magnetic resonance imaging (fMRI) data, a complex spatio-temporal modeling problem must be solved. A crucial part of this process is the estimation of the hemodynamic response (HR) function, an impulse response relating the stimulus signal to the measured noisy response. The estimation of the HR is complicated by the presence of low frequency colored noise. The standard approach to modeling the HR is to use simple parametric models, although FIR models have been used. We offer two contributions. First, we pursue a nonparametric approach using orthonormal causal Laguerre polynomials which have become popular in the system identification literature. It also happens that the shape of the basis elements is similar to that of a typical HR. We thus expect to achieve a compact, and so bias reduced, and low noise representation of the HR. Additionally, we develop a procedure for providing confidence intervals for the whole HR function. This feature is completely lacking in all previous work.

In order to construct spatial activation plots from functional magnetic resonance imaging (fMRI) data, a complex spatio-temporal modeling problem must be solved. A crucial part of this process is the estimation of the hemodynamic response (HR) function, an impulse response relating the stimulus signal to the measured noisy response. The estimation of the HR is complicated by the presence of low frequency colored noise. The standard approach to modeling the HR is to use simple parametric models, although FIR models have been used. We pursue a nonparametric approach using orthonormal causal Laguerre polynomials which have become popular in the system identification literature. It also happens that the shape of the basis elements is similar to that of a typical HR. We thus expect to achieve a compact and so bias reduced and low noise representation of the HR. This is not the case in FIR modeling, because a low FIR order is unable to cover the whole length of the HR over its region of support while a high FIR order results in overestimation of signal and underestimation of noise leading to misleading interpretations.