skip to main content


Title: Regularized joint estimation of related vector autoregressive models
In a number of applications, one has access to high-dimensional time series data on several related subjects. A motivating application area comes from the neuroimaging field, such as brain fMRI time series data, obtained from various groups of subjects (cases/controls) with a specific neurological disorder. The problem of regularized joint estimation of multiple related Vector Autoregressive (VAR) models is discussed, leveraging a group lasso penalty in addition to a regular lasso one, so as to increase statistical efficiency of the estimates by borrowing strength across the models. A modeling framework is developed that it allows for both group-level and subject-specific effects for related subjects, using a group lasso penalty to estimate the former. An estimation procedure is introduced, whose performance is illustrated on synthetic data and compared to other state-of-the-art methods. Moreover, the proposed approach is employed for the analysis of resting state fMRI data. In particular, a group-level descriptive analysis is conducted for brain inter-regional temporal effects of Attention Deficit Hyperactive Disorder (ADHD) patients as opposed to controls, with the data available from the ADHD-200 Global Competition repository.  more » « less
Award ID(s):
1821220
NSF-PAR ID:
10108540
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Computational statistics & data analysis
Volume:
139
ISSN:
1872-7352
Page Range / eLocation ID:
164-177
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Functional magnetic resonance imaging (fMRI) studies have shown altered brain dynamic functional connectivity (DFC) in mental disorders. Here, we aim to explore DFC across a spectrum of symptomatically‐related disorders including bipolar disorder with psychosis (BPP), schizoaffective disorder (SAD), and schizophrenia (SZ). We introduce a group information guided independent component analysis procedure to estimate both group‐level and subject‐specific connectivity states from DFC. Using resting‐state fMRI data of 238 healthy controls (HCs), 140 BPP, 132 SAD, and 113 SZ patients, we identified measures differentiating groups from the whole‐brain DFC and traditional static functional connectivity (SFC), separately. Results show that DFC provided more informative measures than SFC. Diagnosis‐related connectivity states were evident using DFC analysis. For the dominant state consistent across groups, we found 22 instances of hypoconnectivity (with decreasing trends from HC to BPP to SAD to SZ) mainly involving post‐central, frontal, and cerebellar cortices as well as 34 examples of hyperconnectivity (with increasing trends HC through SZ) primarily involving thalamus and temporal cortices. Hypoconnectivities/hyperconnectivities also showed negative/positive correlations, respectively, with clinical symptom scores. Specifically, hypoconnectivities linking postcentral and frontal gyri were significantly negatively correlated with the PANSS positive/negative scores. For frontal connectivities, BPP resembled HC while SAD and SZ were more similar. Three connectivities involving the left cerebellar crus differentiated SZ from other groups and one connection linking frontal and fusiform cortices showed a SAD‐unique change. In summary, our method is promising for assessing DFC and may yield imaging biomarkers for quantifying the dimension of psychosis.Hum Brain Mapp 38:2683–2708, 2017. ©2017 Wiley Periodicals, Inc.

     
    more » « less
  2. Background: Schizophrenia affects around 1% of the global population. Functional connectivity extracted from resting-state functional magnetic resonance imaging (rs-fMRI) has previously been used to study schizophrenia and has great potential to provide novel insights into the disorder. Some studies have shown abnormal functional connectivity in the default mode network (DMN) of individuals with schizophrenia, and more recent studies have shown abnormal dynamic functional connectivity (dFC) in individuals with schizophrenia. However, DMN dFC and the link between abnormal DMN dFC and symptom severity have not been well-characterized. Method: Resting-state fMRI data from subjects with schizophrenia (SZ) and healthy controls (HC) across two datasets were analyzed independently. We captured seven maximally independent subnodes in the DMN by applying group independent component analysis and estimated dFC between subnode time courses using a sliding window approach. A clustering method separated the dFCs into five reoccurring brain states. A feature selection method modeled the difference between SZs and HCs using the state-specific FC features. Finally, we used the transition probability of a hidden Markov model to characterize the link between symptom severity and dFC in SZ subjects. Results: We found decreases in the connectivity of the anterior cingulate cortex (ACC) and increases in the connectivity between the precuneus (PCu) and the posterior cingulate cortex (PCC) (i.e., PCu/PCC) of SZ subjects. In SZ, the transition probability from a state with weaker PCu/PCC and stronger ACC connectivity to a state with stronger PCu/PCC and weaker ACC connectivity increased with symptom severity. Conclusions: To our knowledge, this was the first study to investigate DMN dFC and its link to schizophrenia symptom severity. We identified reproducible neural states in a data-driven manner and demonstrated that the strength of connectivity within those states differed between SZs and HCs. Additionally, we identified a relationship between SZ symptom severity and the dynamics of DMN functional connectivity. We validated our results across two datasets. These results support the potential of dFC for use as a biomarker of schizophrenia and shed new light upon the relationship between schizophrenia and DMN dynamics. 
    more » « less
  3. Objectively differentiating patient mental states based on electrical activity, as opposed to overt behavior, is a fundamental neuroscience problem with medical applications, such as identifying patients in locked-in state vs. coma. Electroencephalography (EEG), which detects millisecond-level changes in brain activity across a range of frequencies, allows for assessment of external stimulus processing by the brain in a non-invasive manner. We applied machine learning methods to 26-channel EEG data of 24 fluent Deaf signers watching videos of sign language sentences (comprehension condition), and the same videos reversed in time (non-comprehension condition), to objectively separate vision-based high-level cognition states. While spectrotemporal parameters of the stimuli were identical in comprehension vs. non-comprehension conditions, the neural responses of participants varied based on their ability to linguistically decode visual data. We aimed to determine which subset of parameters (specific scalp regions or frequency ranges) would be necessary and sufficient for high classification accuracy of comprehension state. Optical flow, characterizing distribution of velocities of objects in an image, was calculated for each pixel of stimulus videos using MATLAB Vision toolbox. Coherence between optical flow in the stimulus and EEG neural response (per video, per participant) was then computed using canonical component analysis with NoiseTools toolbox. Peak correlations were extracted for each frequency for each electrode, participant, and video. A set of standard ML algorithms were applied to the entire dataset (26 channels, frequencies from .2 Hz to 12.4 Hz, binned in 1 Hz increments), with consistent out-of-sample 100% accuracy for frequencies in .2-1 Hz range for all regions, and above 80% accuracy for frequencies < 4 Hz. Sparse Optimal Scoring (SOS) was then applied to the EEG data to reduce the dimensionality of the features and improve model interpretability. SOS with elastic-net penalty resulted in out-of-sample classification accuracy of 98.89%. The sparsity pattern in the model indicated that frequencies between 0.2–4 Hz were primarily used in the classification, suggesting that underlying data may be group sparse. Further, SOS with group lasso penalty was applied to regional subsets of electrodes (anterior, posterior, left, right). All trials achieved greater than 97% out-of-sample classification accuracy. The sparsity patterns from the trials using 1 Hz bins over individual regions consistently indicated frequencies between 0.2–1 Hz were primarily used in the classification, with anterior and left regions performing the best with 98.89% and 99.17% classification accuracy, respectively. While the sparsity pattern may not be the unique optimal model for a given trial, the high classification accuracy indicates that these models have accurately identified common neural responses to visual linguistic stimuli. Cortical tracking of spectro-temporal change in the visual signal of sign language appears to rely on lower frequencies proportional to the N400/P600 time-domain evoked response potentials, indicating that visual language comprehension is grounded in predictive processing mechanisms. 
    more » « less
  4. Abstract Background and Hypothesis

    Risk-taking in specific contexts can be beneficial, leading to rewarding outcomes. Schizophrenia is associated with disadvantageous decision-making, as subjects pursue uncertain risky rewards less than controls. However, it is unclear whether this behavior is associated with more risk sensitivity or less reward incentivization. Matching on demographics and intelligence quotient (IQ), we determined whether risk-taking was more associated with brain activation in regions affiliated with risk evaluation or reward processing.

    Study Design

    Subjects (30 schizophrenia/schizoaffective disorder, 30 controls) completed a modified, fMRI Balloon Analogue Risk Task. Brain activation was modeled during decisions to pursue risky rewards and parametrically modeled according to risk level.

    Study Results

    The schizophrenia group exhibited less risky-reward pursuit despite previous adverse outcomes (Average Explosions; F(1,59) = 4.06, P = .048) but the comparable point at which risk-taking was volitionally discontinued (Adjusted Pumps; F(1,59) = 2.65, P = .11). Less activation was found in schizophrenia via whole brain and region of interest (ROI) analyses in the right (F(1,59) = 14.91, P < 0.001) and left (F(1,59) = 16.34, P < 0.001) nucleus accumbens (NAcc) during decisions to pursue rewards relative to riskiness. Risk-taking correlated with IQ in schizophrenia, but not controls. Path analyses of average ROI activation revealed less statistically determined influence of anterior insula upon dorsal anterior cingulate bilaterally (left: χ2 = 12.73, P < .001; right: χ2 = 9.54, P = .002) during risky reward pursuit in schizophrenia.

    Conclusions

    NAcc activation in schizophrenia varied less according to the relative riskiness of uncertain rewards compared to controls, suggesting aberrations in reward processing. The lack of activation differences in other regions suggests similar risk evaluation. Less insular influence on the anterior cingulate may relate to attenuated salience attribution or inability for risk-related brain region collaboration to sufficiently perceive situational risk.

     
    more » « less
  5. Attention deficit hyperactivity disorder (ADHD) is one of the most common brain disorders among children. This disorder is considered as a big threat for public health and causes attention, focus and organizing difficulties for children and even adults. Since the cause of ADHD is not known yet, data mining algorithms are being used to help discover patterns which discriminate healthy from ADHD subjects. Numerous efforts are underway with the goal of developing classification tools for ADHD diagnosis based on functional and structural magnetic resonance imaging data of the brain. In this paper, we used Eros, which is a technique for computing similarity between two multivariate time series along with k-Nearest-Neighbor classifier, to classify healthy vs ADHD children. We designed a model selection scheme called J-Eros which is able to pick the optimum value of k for k-Nearest-Neighbor from the training data. We applied this technique to the public data provided by ADHD-200 Consortium competition and our results show that J-Eros is capable of discriminating healthy from ADHD children such that we outperformed the best results reported by ADHD-200 competition about 20 percent for two datasets. The implemented code is available as GPL license on GitHub portal of our lab at https://github.com/pcdslab/J-Eros. 
    more » « less