Source author record

Amanda F. Mejia

Amanda F. Mejia appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Methodology Applications Neurons and Cognition

Catalog footprint

What is connected

4works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

A spatial template independent component analysis model for subject-level brain network estimation and inference

Independent component analysis is commonly applied to functional magnetic resonance imaging (fMRI) data to extract independent components (ICs) representing functional brain networks. While ICA produces reliable group-level estimates, single-subject ICA often produces noisy results. Template ICA (tICA) is a hierarchical ICA model using empirical population priors to produce reliable subject-level IC estimates. However, this and other hierarchical ICA models assume unrealistically that subject effects are spatially independent. Here, we propose spatial template ICA (stICA), which incorporates spatial process priors into tICA. This results in greater estimation efficiency of ICs and subject effects. Additionally, the joint posterior distribution can be used to identify engaged areas using an excursions set approach. By leveraging spatial dependencies and avoiding massive multiple comparisons, stICA has high power to detect true effects. We derive an efficient expectation-maximization algorithm to obtain maximum likelihood estimates of the model parameters and posterior moments of the latent fields. Based on analysis of simulated data and fMRI data from the Human Connectome Project, we find that stICA produces estimates that are more accurate and reliable than benchmark approaches, and identifies larger and more reliable areas of engagement. The algorithm is quite tractable, achieving convergence within 7 hours in our fMRI analysis.

preprint2016arXiv

Effects of Scan Length and Shrinkage on Reliability of Resting-State Functional Connectivity in the Human Connectome Project

In this paper, we use data from the Human Connectome Project (N=461) to investigate the effect of scan length on reliability of resting-state functional connectivity (rsFC) estimates produced from resting-state functional magnetic resonance imaging (rsfMRI). Additionally, we study the benefits of empirical Bayes shrinkage, in which subject-level estimates borrow strength from the population average by trading a small increase in bias for a greater reduction in variance. For each subject, we compute raw and shrinkage estimates of rsFC between 300 regions identified through independent components analysis (ICA) based on rsfMRI scans varying from 3 to 30 minutes in length. The time course for each region is determined using dual regression, and rsFC is estimated as the Pearson correlation between each pair of time courses. Shrinkage estimates for each subject are computed as a weighted average between the raw subject-level estimate and the population average estimate, where the weight is determined for each connection by the relationship of within-subject variance to between-subject variance. We find that shrinkage estimates exhibit greater reliability than raw estimates for most connections, with 30-40% improvement using scans less than 10 minutes in length and 10-20% improvement using scans of 20-30 minutes. We also observe significant spatial variability in reliability of both raw and shrinkage estimates, with connections within the default mode and motor networks exhibiting the greatest reliability and between-network connections exhibiting the poorest reliability. We conclude that the scan length required for reliable estimation of rsFC depends on the specific connections of interest, and shrinkage can be used to increase reliability of rsFC, even when produced from long, high-quality rsfMRI scans.

preprint2016arXiv

PCA leverage: outlier detection for high-dimensional functional magnetic resonance imaging data

Outlier detection for high-dimensional (HD) data is a popular topic in modern statistical research. However, one source of HD data that has received relatively little attention is functional magnetic resonance images (fMRI), which consists of hundreds of thousands of measurements sampled at hundreds of time points. At a time when the availability of fMRI data is rapidly growing---primarily through large, publicly available grassroots datasets---automated quality control and outlier detection methods are greatly needed. We propose PCA leverage and demonstrate how it can be used to identify outlying time points in an fMRI run. Furthermore, PCA leverage is a measure of the influence of each observation on the estimation of principal components, which are often of interest in fMRI data. We also propose an alternative measure, PCA robust distance, which is less sensitive to outliers and has controllable statistical properties. The proposed methods are validated through simulation studies and are shown to be highly accurate. We also conduct a reliability study using resting-state fMRI data from the Autism Brain Imaging Data Exchange and find that removal of outliers using the proposed methods results in more reliable estimation of subject-level resting-state networks using ICA.

preprint2015arXiv

Improving Reliability of Subject-Level Resting-State fMRI Parcellation with Shrinkage Estimators

A recent interest in resting state functional magnetic resonance imaging (rsfMRI) lies in subdividing the human brain into anatomically and functionally distinct regions of interest. For example, brain parcellation is often used for defining the network nodes in connectivity studies. While inference has traditionally been performed on group-level data, there is a growing interest in parcellating single subject data. However, this is difficult due to the low signal-to-noise ratio of rsfMRI data, combined with typically short scan lengths. A large number of brain parcellation approaches employ clustering, which begins with a measure of similarity or distance between voxels. The goal of this work is to improve the reproducibility of single-subject parcellation using shrinkage estimators of such measures, allowing the noisy subject-specific estimator to "borrow strength" in a principled manner from a larger population of subjects. We present several empirical Bayes shrinkage estimators and outline methods for shrinkage when multiple scans are not available for each subject. We perform shrinkage on raw intervoxel correlation estimates and use both raw and shrinkage estimates to produce parcellations by performing clustering on the voxels. Our proposed method is agnostic to the choice of clustering method and can be used as a pre-processing step for any clustering algorithm. Using two datasets---a simulated dataset where the true parcellation is known and is subject-specific and a test-retest dataset consisting of two 7-minute rsfMRI scans from 20 subjects---we show that parcellations produced from shrinkage correlation estimates have higher reliability and validity than those produced from raw estimates. Application to test-retest data shows that using shrinkage estimators increases the reproducibility of subject-specific parcellations of the motor cortex by up to 30%.