@article {819, title = {Analytic programming with FMRI data: a quick-start guide for statisticians using R.}, journal = {PLoS One}, volume = {9}, year = {2014}, month = {2014}, pages = {e89470}, abstract = {

Functional magnetic resonance imaging (fMRI) is a thriving field that plays an important role in medical imaging analysis, biological and neuroscience research and practice. This manuscript gives a didactic introduction to the statistical analysis of fMRI data using the R project, along with the relevant R code. The goal is to give statisticians who would like to pursue research in this area a quick tutorial for programming with fMRI data. References of relevant packages and papers are provided for those interested in more advanced analysis.

}, keywords = {Brain, brain mapping, Humans, Magnetic Resonance Imaging, Programming Languages}, issn = {1932-6203}, doi = {10.1371/journal.pone.0089470}, author = {Eloyan, Ani and Li, Shanshan and Muschelli, John and Pekar, Jim J and Mostofsky, Stewart H and Caffo, Brian S} } @article {833, title = {Cross-sectional and longitudinal association of body mass index and brain volume.}, journal = {Hum Brain Mapp}, volume = {35}, year = {2014}, month = {2014 Jan}, pages = {75-88}, abstract = {

Although a link between body mass index (BMI) and brain volume has been established in several cross-sectional studies, evidence of the association between change in BMI over time and changes in brain structure is limited. Using data from a cohort of 347 former lead workers and community controls with two magnetic resonance imaging scans over a period of ~5 years, we estimated cross-sectional and longitudinal associations of BMI and brain volume using both region of interest (ROI) and voxel-based morphometric (VBM) methods. We found that associations of BMI and brain volume were not significantly different in former lead workers when compared with community controls. In the cross-sectional analysis, higher BMIs were associated with smaller brain volumes in gray matter (GM) using both ROI and VBM approaches. No associations with white matter (WM) were observed. In the longitudinal analysis, higher baseline BMI was associated with greater decline in temporal and occipital GM ROI volumes. Change in BMI over the 5-year period was only associated with change in hippocampal volume and was not associated with change in any of the GM ROIs. Overall, higher BMI was associated with lower GM volume in several ROIs and with declines in volume in temporal and occipital GM over time. These results suggest that sustained high body mass may contribute to progressive temporal and occipital atrophy.

}, keywords = {Aged, Atrophy, Body Mass Index, Brain, Cross-Sectional Studies, Humans, Image Interpretation, Computer-Assisted, Longitudinal Studies, Magnetic Resonance Imaging, Male, Middle Aged, Obesity}, issn = {1097-0193}, doi = {10.1002/hbm.22159}, author = {Bobb, Jennifer F and Schwartz, Brian S and Davatzikos, Christos and Caffo, Brian} } @article {831, title = {Disruption of functional organization within the primary motor cortex in children with autism.}, journal = {Hum Brain Mapp}, volume = {35}, year = {2014}, month = {2014 Feb}, pages = {567-80}, abstract = {

Accumulating evidence suggests that motor impairments are prevalent in autism spectrum disorder (ASD), relate to the social and communicative deficits at the core of the diagnosis and may reflect abnormal connectivity within brain networks underlying motor control and learning. Parcellation of resting-state functional connectivity data using spectral clustering approaches has been shown to be an effective means of visualizing functional organization within the brain but has most commonly been applied to explorations of normal brain function. This article presents a parcellation of a key area of the motor network, the primary motor cortex (M1), a key area of the motor control network, in adults, typically developing (TD) children and children with ASD and introduces methods for selecting the number of parcels, matching parcels across groups and testing group differences. The parcellation is based solely on patterns of connectivity between individual M1 voxels and all voxels outside of M1, and within all groups, a gross dorsomedial to ventrolateral organization emerged within M1 which was left-right symmetric. Although this gross organizational scheme was present in both groups of children, statistically significant group differences in the size and segregation of M1 parcels within regions of the motor homunculus corresponding to the upper and lower limbs were observed. Qualitative comparison of the M1 parcellation for children with ASD with that of younger and older TD children suggests that these organizational differences, with a lack of differentiation between lower limb/trunk regions and upper limb/hand regions, may be due, at least in part, to a delay in functional specialization within the motor cortex.

}, keywords = {Adult, Autistic Disorder, brain mapping, Child, Developmental Disabilities, Female, Humans, Image Processing, Computer-Assisted, Magnetic Resonance Imaging, Male, Middle Aged, Motor Cortex, Oxygen, Young Adult}, issn = {1097-0193}, doi = {10.1002/hbm.22188}, author = {Nebel, Mary Beth and Joel, Suresh E and Muschelli, John and Barber, Anita D and Caffo, Brian S and Pekar, James J and Mostofsky, Stewart H} } @article {822, title = {An evaluation of independent component analyses with an application to resting-state fMRI.}, journal = {Biometrics}, volume = {70}, year = {2014}, month = {2014 Mar}, pages = {224-36}, abstract = {

We examine differences between independent component analyses (ICAs) arising from different assumptions, measures of dependence, and starting points of the algorithms. ICA is a popular method with diverse applications including artifact removal in electrophysiology data, feature extraction in microarray data, and identifying brain networks in functional magnetic resonance imaging (fMRI). ICA can be viewed as a generalization of principal component analysis (PCA) that takes into account higher-order cross-correlations. Whereas the PCA solution is unique, there are many ICA methods-whose solutions may differ. Infomax, FastICA, and JADE are commonly applied to fMRI studies, with FastICA being arguably the most popular. Hastie and Tibshirani (2003) demonstrated that ProDenICA outperformed FastICA in simulations with two components. We introduce the application of ProDenICA to simulations with more components and to fMRI data. ProDenICA was more accurate in simulations, and we identified differences between biologically meaningful ICs from ProDenICA versus other methods in the fMRI analysis. ICA methods require nonconvex optimization, yet current practices do not recognize the importance of, nor adequately address sensitivity to, initial values. We found that local optima led to dramatically different estimates in both simulations and group ICA of fMRI, and we provide evidence that the global optimum from ProDenICA is the best estimate. We applied a modification of the Hungarian (Kuhn-Munkres) algorithm to match ICs from multiple estimates, thereby gaining novel insights into how brain networks vary in their sensitivity to initial values and ICA method.

}, keywords = {Adolescent, Algorithms, Attention Deficit Disorder with Hyperactivity, brain mapping, Child, Child, Preschool, Computer Simulation, Humans, Magnetic Resonance Imaging, Models, Statistical, Principal Component Analysis}, issn = {1541-0420}, doi = {10.1111/biom.12111}, author = {Risk, Benjamin B and Matteson, David S and Ruppert, David and Eloyan, Ani and Caffo, Brian S} } @article {684, title = {Corrected confidence bands for functional data using principal components.}, journal = {Biometrics}, volume = {69}, year = {2013}, month = {2013 Mar}, pages = {41-51}, abstract = {

Functional principal components (FPC) analysis is widely used to decompose and express functional observations. Curve estimates implicitly condition on basis functions and other quantities derived from FPC decompositions; however these objects are unknown in practice. In this article, we propose a method for obtaining correct curve estimates by accounting for uncertainty in FPC decompositions. Additionally, pointwise and simultaneous confidence intervals that account for both model- and decomposition-based variability are constructed. Standard mixed model representations of functional expansions are used to construct curve estimates and variances conditional on a specific decomposition. Iterated expectation and variance formulas combine model-based conditional estimates across the distribution of decompositions. A bootstrap procedure is implemented to understand the uncertainty in principal component decomposition quantities. Our method compares favorably to competing approaches in simulation studies that include both densely and sparsely observed functions. We apply our method to sparse observations of CD4 cell counts and to dense white-matter tract profiles. Code for the analyses and simulations is publicly available, and our method is implemented in the R package refund on CRAN.

}, keywords = {Brain, CD4 Lymphocyte Count, Computer Simulation, Confidence Intervals, HIV, HIV Infections, Humans, Magnetic Resonance Imaging, Models, Statistical, Multiple Sclerosis, Principal Component Analysis}, issn = {1541-0420}, doi = {10.1111/j.1541-0420.2012.01808.x}, author = {Goldsmith, J and Greven, S and Crainiceanu, C} } @article {824, title = {A decision-theory approach to interpretable set analysis for high-dimensional data.}, journal = {Biometrics}, volume = {69}, year = {2013}, month = {2013 Sep}, pages = {614-23}, abstract = {

A key problem in high-dimensional significance analysis is to find pre-defined sets that show enrichment for a statistical signal of interest; the classic example is the enrichment of gene sets for differentially expressed genes. Here, we propose a new decision-theory approach to the analysis of gene sets which focuses on estimating the fraction of non-null variables in a set. We introduce the idea of "atoms," non-overlapping sets based on the original pre-defined set annotations. Our approach focuses on finding the union of atoms that minimizes a weighted average of the number of false discoveries and missed discoveries. We introduce a new false discovery rate for sets, called the atomic false discovery rate (afdr), and prove that the optimal estimator in our decision-theory framework is to threshold the afdr. These results provide a coherent and interpretable framework for the analysis of sets that addresses the key issues of overlapping annotations and difficulty in interpreting p values in both competitive and self-contained tests. We illustrate our method and compare it to a popular existing method using simulated examples, as well as gene-set and brain ROI data analyses.

}, keywords = {Algorithms, Bayes Theorem, Biometry, Brain, Computer Simulation, Data Interpretation, Statistical, Decision Theory, Functional Neuroimaging, Gene Expression Profiling, Genomics, Humans, Magnetic Resonance Imaging, Models, Statistical, Oligonucleotide Array Sequence Analysis}, issn = {1541-0420}, doi = {10.1111/biom.12060}, author = {Boca, Simina M and Bravo, H{\'e}ctor C{\'e}orrada and Caffo, Brian and Leek, Jeffrey T and Parmigiani, Giovanni} } @article {829, title = {Developmental changes in within- and between-network connectivity between late childhood and adulthood.}, journal = {Neuropsychologia}, volume = {51}, year = {2013}, month = {2013 Jan}, pages = {156-67}, abstract = {

A number of behavioral changes occur between late childhood and adulthood, including maturation of social cognition, reward receptivity, impulsiveness, risk-taking and cognitive control. Although some of these abilities show linear improvements with age, some abilities may temporarily worsen, reflecting both the restructuring and/or strengthening of connections within some brain systems. The current study uses resting state functional connectivity to examine developmental differences between late childhood and adulthood in task positive (TP) regions, which play a role in cognitive control functions, and task negative (TN) regions, which play a role in social cognition, self-referential, and internally-directed thought. Within the TP network, developmental differences in connectivity were found with the left dorsolateral prefrontal cortex. Within the TN network, developmental differences in connectivity were found with a broad area of the medial prefrontal cortex and the right parahippocampal gyrus. Connections between the two networks also showed significant developmental differences. Stronger anticorrelations were found in the TN maps of the adult group for the right anterior insula/inferior frontal gyrus, bilateral anterior inferior parietal lobule, bilateral superior parietal lobule and an anterior portion of the right posterior cingulate cortex. There was a significant brain-behavior relationship between the strength of anticorrelation in these regions and inhibitory control performance on two Go/No-go tasks suggesting that the development of anticorrelations between late childhood and adulthood supports mature inhibitory control. Overall, maturation of these networks occurred in specific regions which are associated with cognitive control of goal-directed behavior, including those involved in working memory, social cognition, and inhibitory control.

}, keywords = {Adolescent, Adult, Age Factors, Analysis of Variance, Brain, brain mapping, Child, Child Development, Cognition, Female, Humans, Image Processing, Computer-Assisted, Inhibition (Psychology), Magnetic Resonance Imaging, Male, Middle Aged, Neural Pathways, Neuropsychological Tests, Oxygen, Reaction Time, Young Adult}, issn = {1873-3514}, doi = {10.1016/j.neuropsychologia.2012.11.011}, author = {Barber, Anita D and Caffo, Brian S and Pekar, James J and Mostofsky, Stewart H} } @article {826, title = {Effects of working memory demand on neural mechanisms of motor response selection and control.}, journal = {J Cogn Neurosci}, volume = {25}, year = {2013}, month = {2013 Aug}, pages = {1235-48}, abstract = {

Inhibitory control commonly recruits a number of frontal regions: pre-supplementary motor area (pre-SMA), frontal eye fields (FEFs), and right-lateralized posterior inferior frontal gyrus (IFG), dorsal anterior insula (DAI), dorsolateral prefrontal cortex (DLPFC), and inferior frontal junction (IFJ). These regions may directly implement inhibitory motor control or may be more generally involved in executive control functions. Two go/no-go tasks were used to distinguish regions specifically recruited for inhibition from those that additionally show increased activity with working memory demand. The pre-SMA and IFG were recruited for inhibition in both tasks and did not have greater activation for working memory demand on no-go trials, consistent with a role in inhibitory control. Activation in pre-SMA also responded to response selection demand and was increased with working memory on go trials specifically. The bilateral FEF and right DAI were commonly active for no-go trials. The FEF was also recruited to a greater degree with working memory demand on go trials and may bias top-down information when stimulus-response mappings change. The DAI, additionally responded to increased working memory demand on both go and no-go trials and may be involved in accessing sustained task information, alerting, or autonomic changes when cognitive demands increase. DLPFC activation was consistent with a role in working memory retrieval on both go and no-go trials. The inferior frontal junction, on the other hand, had greater activation with working memory specifically for no-go trials and may detect salient stimuli when the task requires frequent updating of working memory representations.

}, keywords = {Adult, Brain, brain mapping, Choice Behavior, Female, Humans, Image Processing, Computer-Assisted, Inhibition (Psychology), Magnetic Resonance Imaging, Male, Memory, Short-Term, Motor Activity, Neuropsychological Tests, Oxygen, Reaction Time, Young Adult}, issn = {1530-8898}, doi = {10.1162/jocn_a_00394}, author = {Barber, Anita D and Caffo, Brian S and Pekar, James J and Mostofsky, Stewart H} } @article {839, title = {An fMRI-based neurologic signature of physical pain.}, journal = {N Engl J Med}, volume = {368}, year = {2013}, month = {2013 Apr 11}, pages = {1388-97}, abstract = {

BACKGROUND: Persistent pain is measured by means of self-report, the sole reliance on which hampers diagnosis and treatment. Functional magnetic resonance imaging (fMRI) holds promise for identifying objective measures of pain, but brain measures that are sensitive and specific to physical pain have not yet been identified.

METHODS: In four studies involving a total of 114 participants, we developed an fMRI-based measure that predicts pain intensity at the level of the individual person. In study 1, we used machine-learning analyses to identify a pattern of fMRI activity across brain regions--a neurologic signature--that was associated with heat-induced pain. The pattern included the thalamus, the posterior and anterior insulae, the secondary somatosensory cortex, the anterior cingulate cortex, the periaqueductal gray matter, and other regions. In study 2, we tested the sensitivity and specificity of the signature to pain versus warmth in a new sample. In study 3, we assessed specificity relative to social pain, which activates many of the same brain regions as physical pain. In study 4, we assessed the responsiveness of the measure to the analgesic agent remifentanil.

RESULTS: In study 1, the neurologic signature showed sensitivity and specificity of 94\% or more (95\% confidence interval [CI], 89 to 98) in discriminating painful heat from nonpainful warmth, pain anticipation, and pain recall. In study 2, the signature discriminated between painful heat and nonpainful warmth with 93\% sensitivity and specificity (95\% CI, 84 to 100). In study 3, it discriminated between physical pain and social pain with 85\% sensitivity (95\% CI, 76 to 94) and 73\% specificity (95\% CI, 61 to 84) and with 95\% sensitivity and specificity in a forced-choice test of which of two conditions was more painful. In study 4, the strength of the signature response was substantially reduced when remifentanil was administered.

CONCLUSIONS: It is possible to use fMRI to assess pain elicited by noxious heat in healthy persons. Future studies are needed to assess whether the signature predicts clinical pain. (Funded by the National Institute on Drug Abuse and others.).

}, keywords = {Adult, Analgesics, Opioid, Artificial Intelligence, Brain, brain mapping, Female, Hot Temperature, Humans, Magnetic Resonance Imaging, Male, Pain, Pain Measurement, Piperidines, ROC Curve, Sensitivity and Specificity, Young Adult}, issn = {1533-4406}, doi = {10.1056/NEJMoa1204471}, author = {Wager, Tor D and Atlas, Lauren Y and Lindquist, Martin A and Roy, Mathieu and Woo, Choong-Wan and Kross, Ethan} } @article {812, title = {Likelihood-based population independent component analysis.}, journal = {Biostatistics}, volume = {14}, year = {2013}, month = {2013 Jul}, pages = {514-27}, abstract = {

Independent component analysis (ICA) is a widely used technique for blind source separation, used heavily in several scientific research areas including acoustics, electrophysiology, and functional neuroimaging. We propose a scalable two-stage iterative true group ICA methodology for analyzing population level functional magnetic resonance imaging (fMRI) data where the number of subjects is very large. The method is based on likelihood estimators of the underlying source densities and the mixing matrix. As opposed to many commonly used group ICA algorithms, the proposed method does not require significant data reduction by a 2-fold singular value decomposition. In addition, the method can be applied to a large group of subjects since the memory requirements are not restrictive. The performance of our approach is compared with a commonly used group ICA algorithm via simulation studies. Furthermore, the proposed method is applied to a large collection of resting state fMRI datasets. The results show that established brain networks are well recovered by the proposed algorithm.

}, keywords = {Adult, Algorithms, Attention Deficit Disorder with Hyperactivity, Biostatistics, Brain, Child, Computer Simulation, Databases, Factual, Humans, Imaging, Three-Dimensional, Likelihood Functions, Magnetic Resonance Imaging, Principal Component Analysis, Statistics, Nonparametric}, issn = {1468-4357}, doi = {10.1093/biostatistics/kxs055}, author = {Eloyan, Ani and Crainiceanu, Ciprian M and Caffo, Brian S} } @article {811, title = {Relationships between retinal axonal and neuronal measures and global central nervous system pathology in multiple sclerosis.}, journal = {JAMA Neurol}, volume = {70}, year = {2013}, month = {2013 Jan}, pages = {34-43}, abstract = {

OBJECTIVE: To determine the relationships between conventional and segmentation-derived optical coherence tomography (OCT) retinal layer thickness measures with intracranial volume (a surrogate of head size) and brain substructure volumes in multiple sclerosis (MS).

DESIGN: Cross-sectional study.

SETTING: Johns Hopkins University, Baltimore, Maryland.

PARTICIPANTS: A total of 84 patients with MS and 24 healthy control subjects.

MAIN OUTCOME MEASURES: High-definition spectral-domain OCT conventional and automated segmentation-derived discrete retinal layer thicknesses and 3-T magnetic resonance imaging brain substructure volumes.

RESULTS: Peripapillary retinal nerve fiber layer as well as composite ganglion cell layer+inner plexiform layer thicknesses in the eyes of patients with MS without a history of optic neuritis were associated with cortical gray matter (P=.01 and P=.04, respectively) and caudate (P=.04 and P=.03, respectively) volumes. Inner nuclear layer thickness, also in eyes without a history of optic neuritis, was associated with fluid-attenuated inversion recovery lesion volume (P=.007) and inversely associated with normal-appearing white matter volume (P=.005) in relapsing-remitting MS. As intracranial volume was found to be related with several of the OCT measures in patients with MS and healthy control subjects and is already known to be associated with brain substructure volumes, all OCT-brain substructure relationships were adjusted for intracranial volume. CONCLUSIONS Retinal measures reflect global central nervous system pathology in multiple sclerosis, with thicknesses of discrete retinal layers each appearing to be associated with distinct central nervous system processes. Moreover, OCT measures appear to correlate with intracranial volume in patients with MS and healthy control subjects, an important unexpected factor unaccounted for in prior studies examining the relationships between peripapillary retinal nerve fiber layer thickness and brain substructure volumes.

}, keywords = {Adult, Axons, Caudate Nucleus, Central Nervous System, Cerebral Cortex, Cross-Sectional Studies, Female, Humans, Magnetic Resonance Imaging, Male, Middle Aged, Multiple Sclerosis, Retina, Retinal Neurons, Tomography, Optical Coherence}, issn = {2168-6157}, doi = {10.1001/jamaneurol.2013.573}, author = {Saidha, Shiv and Sotirchos, Elias S and Oh, Jiwon and Syc, Stephanie B and Seigo, Michaela A and Shiee, Navid and Eckstein, Chistopher and Durbin, Mary K and Oakley, Jonathan D and Meyer, Scott A and Frohman, Teresa C and Newsome, Scott and Ratchford, John N and Balcer, Laura J and Pham, Dzung L and Crainiceanu, Ciprian M and Frohman, Elliot M and Reich, Daniel S and Calabresi, Peter A} } @article {827, title = {Analysis of group ICA-based connectivity measures from fMRI: application to Alzheimer{\textquoteright}s disease.}, journal = {PLoS One}, volume = {7}, year = {2012}, month = {2012}, pages = {e49340}, abstract = {

Functional magnetic resonance imaging (fMRI) is a powerful tool for the in vivo study of the pathophysiology of brain disorders and disease. In this manuscript, we propose an analysis stream for fMRI functional connectivity data and apply it to a novel study of Alzheimer{\textquoteright}s disease. In the first stage, spatial independent component analysis is applied to group fMRI data to obtain common brain networks (spatial maps) and subject-specific mixing matrices (time courses). In the second stage, functional principal component analysis is utilized to decompose the mixing matrices into population-level eigenvectors and subject-specific loadings. Inference is performed using permutation-based exact logistic regression for matched pairs data. The method is applied to a novel fMRI study of Alzheimer{\textquoteright}s disease risk under a verbal paired associates task. We found empirical evidence of alternative ICA-based metrics of connectivity when comparing subjects evidencing mild cognitive impairment relative to carefully matched controls.

}, keywords = {Aged, Algorithms, Alzheimer Disease, Brain, brain mapping, Case-Control Studies, Female, Humans, Image Processing, Computer-Assisted, Magnetic Resonance Imaging, Male, Middle Aged, Mild Cognitive Impairment, Nerve Net, Principal Component Analysis, Task Performance and Analysis, Time Factors}, issn = {1932-6203}, doi = {10.1371/journal.pone.0049340}, author = {Li, Shanshan and Eloyan, Ani and Joel, Suresh and Mostofsky, Stewart and Pekar, James and Bassett, Susan Spear and Caffo, Brian} } @article {841, title = {Dissociable influences of opiates and expectations on pain.}, journal = {J Neurosci}, volume = {32}, year = {2012}, month = {2012 Jun 6}, pages = {8053-64}, abstract = {

Placebo treatments and opiate drugs are thought to have common effects on the opioid system and pain-related brain processes. This has created excitement about the potential for expectations to modulate drug effects themselves. If drug effects differ as a function of belief, this would challenge the assumptions underlying the standard clinical trial. We conducted two studies to directly examine the relationship between expectations and opioid analgesia. We administered the opioid agonist remifentanil to human subjects during experimental thermal pain and manipulated participants{\textquoteright} knowledge of drug delivery using an open-hidden design. This allowed us to test drug effects, expectancy (knowledge) effects, and their interactions on pain reports and pain-related responses in the brain. Remifentanil and expectancy both reduced pain, but drug effects on pain reports and fMRI activity did not interact with expectancy. Regions associated with pain processing showed drug-induced modulation during both Open and Hidden conditions, with no differences in drug effects as a function of expectation. Instead, expectancy modulated activity in frontal cortex, with a separable time course from drug effects. These findings reveal that opiates and placebo treatments both influence clinically relevant outcomes and operate without mutual interference.

}, keywords = {Analgesics, Opioid, Anticipation, Psychological, Behavior, brain mapping, Dose-Response Relationship, Drug, Female, Hemodynamics, Hot Temperature, Humans, Image Processing, Computer-Assisted, Injections, Intravenous, Linear Models, Magnetic Resonance Imaging, Male, Pain, Pain Measurement, Piperidines, Young Adult}, issn = {1529-2401}, doi = {10.1523/JNEUROSCI.0383-12.2012}, author = {Atlas, Lauren Y and Whittington, Robert A and Lindquist, Martin A and Wielgosz, Joe and Sonty, Nomita and Wager, Tor D} } @article {842, title = {Dynamic connectivity regression: determining state-related changes in brain connectivity.}, journal = {Neuroimage}, volume = {61}, year = {2012}, month = {2012 Jul 16}, pages = {907-20}, abstract = {

Most statistical analyses of fMRI data assume that the nature, timing and duration of the psychological processes being studied are known. However, often it is hard to specify this information a priori. In this work we introduce a data-driven technique for partitioning the experimental time course into distinct temporal intervals with different multivariate functional connectivity patterns between a set of regions of interest (ROIs). The technique, called Dynamic Connectivity Regression (DCR), detects temporal change points in functional connectivity and estimates a graph, or set of relationships between ROIs, for data in the temporal partition that falls between pairs of change points. Hence, DCR allows for estimation of both the time of change in connectivity and the connectivity graph for each partition, without requiring prior knowledge of the nature of the experimental design. Permutation and bootstrapping methods are used to perform inference on the change points. The method is applied to various simulated data sets as well as to an fMRI data set from a study (N=26) of a state anxiety induction using a socially evaluative threat challenge. The results illustrate the method{\textquoteright}s ability to observe how the networks between different brain regions changed with subjects{\textquoteright} emotional state.

}, keywords = {Algorithms, Brain, brain mapping, Computer Simulation, Emotions, Humans, Image Interpretation, Computer-Assisted, Magnetic Resonance Imaging, Neural Pathways}, issn = {1095-9572}, doi = {10.1016/j.neuroimage.2012.03.070}, author = {Cribben, Ivor and Haraldsdottir, Ragnheidur and Atlas, Lauren Y and Wager, Tor D and Lindquist, Martin A} } @article {844, title = {Estimating and testing variance components in a multi-level GLM.}, journal = {Neuroimage}, volume = {59}, year = {2012}, month = {2012 Jan 2}, pages = {490-501}, abstract = {

Most analysis of multi-subject fMRI data is concerned with determining whether there exists a significant population-wide {\textquoteright}activation{\textquoteright} in a comparison between two or more conditions. This is typically assessed by testing the average value of a contrast of parameter estimates (COPE) against zero in a general linear model (GLM) analysis. However, important information can also be obtained by testing whether there exist significant individual differences in effect magnitude between subjects, i.e. whether the variance of a COPE is significantly different from zero. Intuitively, such a test amounts to testing whether inter-individual differences are larger than would be expected given the within-subject error variance. We compare several methods for estimating variance components, including a) a na{\"\i}ve estimate using ordinary least squares (OLS); b) linear mixed effects in R (LMER); c) a novel Matlab implementation of iterative generalized least squares (IGLS) and its restricted maximum likelihood variant (RIGLS). All methods produced reasonable estimates of within- and between-subject variance components, with IGLS providing an attractive balance between sensitivity and appropriate control of false positives. Finally, we use the IGLS method to estimate inter-subject variance in a perfusion fMRI study (N=18) of social evaluative threat, and show evidence for significant inter-individual differences in ventromedial prefrontal cortex (VMPFC), amygdala, hippocampus and medial temporal lobes, insula, and brainstem, with predicted inverse coupling between VMPFC and the midbrain periaqueductal gray only when high inter-individual variance was used to define the seed for functional connectivity analyses. In sum, tests of variance provides a way of selecting regions that show significant inter-individual variability for subsequent analyses that attempt to explain those individual differences.

}, keywords = {Brain, brain mapping, Humans, Image Interpretation, Computer-Assisted, Linear Models, Magnetic Resonance Imaging, Models, Neurological, Sensitivity and Specificity}, issn = {1095-9572}, doi = {10.1016/j.neuroimage.2011.07.077}, author = {Lindquist, Martin A and Spicer, Julie and Asllani, Iris and Wager, Tor D} } @article {830, title = {Total white matter hyperintensity volume in bipolar disorder patients and their healthy relatives.}, journal = {Bipolar Disord}, volume = {14}, year = {2012}, month = {2012 Dec}, pages = {888-93}, abstract = {

OBJECTIVES: White matter hyperintensities (WMH) are more common in subjects with bipolar disorder (BP) than in healthy subjects (HS). Few studies have examined the effect of the diagnostic type of bipolar illness on WMH burden, and none have approached this question through a direct measurement of the volume of affected white matter in relationship to familiality. In this pilot study, we utilized a volumetric measurement of WMH to investigate the relationship between the total volume of WMH and the familiality and type of BP.

METHODS: Forty-five individuals with bipolar I disorder (BP-I) with psychotic features, BP-I without psychotic features, or bipolar II disorder (BP-II), seven of their unaffected relatives, and 32 HS were recruited for participation. T-2 weighted magnetic resonance imaging scans were obtained on all subjects, and the total volume of all WMH for each subject was measured in cubic centimeters. The significance of difference between groups was tested using ANOVA with post-hoc adjustment for multiple comparisons. Further, we used logistic regression to test for trends between symptom load and total WMH volume.

RESULTS: The mean total volume of WMH in BP-I patients with psychotic features was significantly higher (p < 0.05) than that of HS. Further, we observed a positive linear trend by familiality and type of affectedness when comparing mean total WMH volume of HS, unaffected family members, subjects with BP-II, and BP-I with and without a history of psychosis (p < 0.05).

CONCLUSIONS: Based on a quantitative technique, WMH burden appears to be associated with familiality and type of BP. The significance of these findings remains to be fully elucidated.

}, keywords = {Adult, Analysis of Variance, Bipolar Disorder, Brain, Female, Humans, Image Processing, Computer-Assisted, Leukoencephalopathies, Magnetic Resonance Imaging, Male, Nerve Fibers, Myelinated}, issn = {1399-5618}, doi = {10.1111/bdi.12019}, author = {Tighe, Sarah K and Reading, Sarah A and Rivkin, Paul and Caffo, Brian and Schweizer, Barbara and Pearlson, Godfrey and Potash, James B and Depaulo, J Raymond and Bassett, Susan S} }