Statistical techniques are needed to analyze data structures with complex dependencies such that clinically useful information can be extracted. Individual‐specific networks, which capture ...dependencies in complex biological systems, are often summarized by graph‐theoretical features. These features, which lend themselves to outcome modeling, can be subject to high variability due to arbitrary decisions in network inference and noise. Correlation‐based adjacency matrices often need to be sparsified before meaningful graph‐theoretical features can be extracted, requiring the data analysts to determine an optimal threshold. To address this issue, we propose to incorporate a flexible weighting function over the full range of possible thresholds to capture the variability of graph‐theoretical features over the threshold domain. The potential of this approach, which extends concepts from functional data analysis to a graph‐theoretical setting, is explored in a plasmode simulation study using real functional magnetic resonance imaging (fMRI) data from the Autism Brain Imaging Data Exchange (ABIDE) Preprocessed initiative. The simulations show that our modeling approach yields accurate estimates of the functional form of the weight function, improves inference efficiency, and achieves a comparable or reduced root mean square prediction error compared to competitor modeling approaches. This assertion holds true in settings where both complex functional forms underlie the outcome‐generating process and a universal threshold value is employed. We demonstrate the practical utility of our approach by using resting‐state fMRI data to predict biological age in children. Our study establishes the flexible modeling approach as a statistically principled, serious competitor to ad‐hoc methods with superior performance.
Working memory is a complex psychological construct referring to the temporary storage and active processing of information. We used functional connectivity brain network metrics quantifying local ...and global efficiency of information transfer for predicting individual variability in working memory performance on an n-back task in both young (n = 14) and older (n = 15) adults. Individual differences in both local and global efficiency during the working memory task were significant predictors of working memory performance in addition to age (and an interaction between age and global efficiency). Decreases in local efficiency during the working memory task were associated with better working memory performance in both age cohorts. In contrast, increases in global efficiency were associated with much better working performance for young participants; however, increases in global efficiency were associated with a slight decrease in working memory performance for older participants. Individual differences in local and global efficiency during resting-state sessions were not significant predictors of working memory performance. Significant group whole-brain functional network decreases in local efficiency also were observed during the working memory task compared to rest, whereas no significant differences were observed in network global efficiency. These results are discussed in relation to recently developed models of age-related differences in working memory.
Brain network analyses have exploded in recent years and hold great potential in helping us understand normal and abnormal brain function. Network science approaches have facilitated these analyses ...and our understanding of how the brain is structurally and functionally organized. However, the development of statistical methods that allow relating this organization to health outcomes has lagged behind. We have attempted to address this need by developing mixed modeling frameworks that allow relating system-level properties of brain networks to outcomes of interest. These frameworks serve as a synergistic fusion of multivariate statistical approaches with network science, providing a needed analytic (modeling and inferential) foundation for whole-brain network data. In this chapter we delineate these approaches that have been developed for single-task and multitask (longitudinal) brain network data, illustrate their utility with data applications, detail their implementation with a user-friendly Matlab toolbox, and discuss ongoing work to adapt the methods to (within-task) dynamic network analysis.
Exponential random graph models (ERGMs), also known as p* models, have been utilized extensively in the social science literature to study complex networks and how their global structure depends on ...underlying structural components. However, the literature on their use in biological networks (especially brain networks) has remained sparse. Descriptive models based on a specific feature of the graph (clustering coefficient, degree distribution, etc.) have dominated connectivity research in neuroscience. Corresponding generative models have been developed to reproduce one of these features. However, the complexity inherent in whole-brain network data necessitates the development and use of tools that allow the systematic exploration of several features simultaneously and how they interact to form the global network architecture. ERGMs provide a statistically principled approach to the assessment of how a set of interacting local brain network features gives rise to the global structure. We illustrate the utility of ERGMs for modeling, analyzing, and simulating complex whole-brain networks with network data from normal subjects. We also provide a foundation for the selection of important local features through the implementation and assessment of three selection approaches: a traditional p-value based backward selection approach, an information criterion approach (AIC), and a graphical goodness of fit (GOF) approach. The graphical GOF approach serves as the best method given the scientific interest in being able to capture and reproduce the structure of fitted brain networks.
Objective
The goal of this study was to determine whether the degree of weight loss after 6 months of a behavior‐based intervention is related to baseline connectivity within two functional networks ...(FNs) of interest, FN1 and FN2, in a group of older adults with obesity.
Methods
Baseline functional magnetic resonance imaging data were collected following an overnight fast in 71 older adults with obesity involved in a weight‐loss intervention. Functional brain networks in a resting state and during a food‐cue task were analyzed using a mixed‐regression framework to examine the relationships between baseline networks and 6‐month change in weight.
Results
During the resting condition, the relationship of baseline brain functional connectivity and network clustering in FN1, which includes the visual cortex and sensorimotor areas, was significantly associated with 6‐month weight loss. During the food‐cue condition, 6‐month weight loss was significantly associated with the relationship between baseline brain connectivity and network global efficiency in FN2, which includes executive control, attention, and limbic regions.
Conclusion
These findings provide further insight into complex functional circuits in the brain related to successful weight loss and may ultimately aid in developing tailored behavior‐based treatment regimens that target specific brain circuitry.
Statistical models to predict incident diabetes are often based on limited variables. Here we pursued two main goals: 1) investigate the relative performance of a machine learning method such as ...Random Forests (RF) for detecting incident diabetes in a high-dimensional setting defined by a large set of observational data, and 2) uncover potential predictors of diabetes. The Jackson Heart Study collected data at baseline and in two follow-up visits from 5,301 African Americans. We excluded those with baseline diabetes and no follow-up, leaving 3,633 individuals for analyses. Over a mean 8-year follow-up, 584 participants developed diabetes. The full RF model evaluated 93 variables including demographic, anthropometric, blood biomarker, medical history, and echocardiogram data. We also used RF metrics of variable importance to rank variables according to their contribution to diabetes prediction. We implemented other models based on logistic regression and RF where features were preselected. The RF full model performance was similar (AUC = 0.82) to those more parsimonious models. The top-ranked variables according to RF included hemoglobin A1C, fasting plasma glucose, waist circumference, adiponectin, c-reactive protein, triglycerides, leptin, left ventricular mass, high-density lipoprotein cholesterol, and aldosterone. This work shows the potential of RF for incident diabetes prediction while dealing with high-dimensional data.
Group-based brain connectivity networks have great appeal for researchers interested in gaining further insight into complex brain function and how it changes across different mental states and ...disease conditions. Accurately constructing these networks presents a daunting challenge given the difficulties associated with accounting for inter-subject topological variability. Viable approaches to this task must engender networks that capture the constitutive topological properties of the group of subjects' networks that it is aiming to represent. The conventional approach has been to use a mean or median correlation network (Achard et al., 2006; Song et al., 2009; Zuo et al., 2011) to embody a group of networks. However, the degree to which their topological properties conform with those of the groups that they are purported to represent has yet to be explored. Here we investigate the performance of these mean and median correlation networks. We also propose an alternative approach based on an exponential random graph modeling framework and compare its performance to that of the aforementioned conventional approach. Simpson et al. (2011) illustrated the utility of exponential random graph models (ERGMs) for creating brain networks that capture the topological characteristics of a single subject's brain network. However, their advantageousness in the context of producing a brain network that “represents” a group of brain networks has yet to be examined. Here we show that our proposed ERGM approach outperforms the conventional mean and median correlation based approaches and provides an accurate and flexible method for constructing group-based representative brain networks.
► Averaging connectivity values across subjects leads to a poor group network. ► Exponential random graph models are applied to construct group-based networks. ► ERGM group-based networks capture important properties of individual networks. ► ERGM group-based approach outperforms the mean/median approach.► These group networks can be used for many purposes including modularity analyses.
Standard methods for quantifying IncuCyte ZOOM(™) assays involve measurements that quantify how rapidly the initially-vacant area becomes re-colonised with cells as a function of time. Unfortunately, ...these measurements give no insight into the details of the cellular-level mechanisms acting to close the initially-vacant area. We provide an alternative method enabling us to quantify the role of cell motility and cell proliferation separately. To achieve this we calibrate standard data available from IncuCyte ZOOM(™) images to the solution of the Fisher-Kolmogorov model.
The Fisher-Kolmogorov model is a reaction-diffusion equation that has been used to describe collective cell spreading driven by cell migration, characterised by a cell diffusivity, D, and carrying capacity limited proliferation with proliferation rate, λ, and carrying capacity density, K. By analysing temporal changes in cell density in several subregions located well-behind the initial position of the leading edge we estimate λ and K. Given these estimates, we then apply automatic leading edge detection algorithms to the images produced by the IncuCyte ZOOM(™) assay and match this data with a numerical solution of the Fisher-Kolmogorov equation to provide an estimate of D. We demonstrate this method by applying it to interpret a suite of IncuCyte ZOOM(™) assays using PC-3 prostate cancer cells and obtain estimates of D, λ and K. Comparing estimates of D, λ and K for a control assay with estimates of D, λ and K for assays where epidermal growth factor (EGF) is applied in varying concentrations confirms that EGF enhances the rate of scratch closure and that this stimulation is driven by an increase in D and λ, whereas K is relatively unaffected by EGF.
Our approach for estimating D, λ and K from an IncuCyte ZOOM(™) assay provides more detail about cellular-level behaviour than standard methods for analysing these assays. In particular, our approach can be used to quantify the balance of cell migration and cell proliferation and, as we demonstrate, allow us to quantify how the addition of growth factors affects these processes individually.
The sliding window correlation (SWC) analysis is a straightforward and common approach for evaluating dynamic functional connectivity. Despite the fact that sliding window analyses have been long ...used, there are still considerable technical issues associated with the approach. A great effort has recently been dedicated to investigate the window setting effects on dynamic connectivity estimation. In this direction, tapered windows have been proposed to alleviate the effect of sudden changes associated with the edges of rectangular windows. Nevertheless, the majority of the windows exploited to estimate brain connectivity tend to suppress dynamic correlations, especially those with faster variations over time. Here, we introduced a window named modulated rectangular (mRect) to address the suppressing effect associated with the conventional windows. We provided a frequency domain analysis using simulated time series to investigate how sliding window analysis (using the regular window functions, e.g. rectangular and tapered windows) may lead to unwanted spectral modulations, and then we showed how this issue can be alleviated through the mRect window. Moreover, we created simulated dynamic network data with altering states over time using simulated fMRI time series, to examine the performance of different windows in tracking network states. We quantified the state identification rate of different window functions through the Jaccard index, and observed superior performance of the mRect window compared to the conventional window functions. Overall, the proposed window function provides an approach that improves SWC estimations, and thus the subsequent inferences and interpretations based on the connectivity network analyses.
•Sliding window analysis has been widely used for quantifying dynamic connectivity.•Sliding window operation is characterized by a non-uniform frequency spectrum.•Sliding window operation thus suppresses fast dynamic correlations.•We modulated rectangular window to make a new window with uniform spectrum.•The proposed window outperformed regular windows in network connectivity analysis.