Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2013 Nov 1:81:358-370.
doi: 10.1016/j.neuroimage.2013.05.049. Epub 2013 May 20.

Spatiotemporal linear mixed effects modeling for the mass-univariate analysis of longitudinal neuroimage data

Affiliations

Spatiotemporal linear mixed effects modeling for the mass-univariate analysis of longitudinal neuroimage data

Jorge L Bernal-Rusiel et al. Neuroimage. .

Erratum in

  • Neuroimage. 2015 Mar;108:123

Abstract

We present an extension of the Linear Mixed Effects (LME) modeling approach to be applied to the mass-univariate analysis of longitudinal neuroimaging (LNI) data. The proposed method, called spatiotemporal LME or ST-LME, builds on the flexible LME framework and exploits the spatial structure in image data. We instantiated ST-LME for the analysis of cortical surface measurements (e.g. thickness) computed by FreeSurfer, a widely-used brain Magnetic Resonance Image (MRI) analysis software package. We validate the proposed ST-LME method and provide a quantitative and objective empirical comparison with two popular alternative methods, using two brain MRI datasets obtained from the Alzheimer's disease neuroimaging initiative (ADNI) and Open Access Series of Imaging Studies (OASIS). Our experiments revealed that ST-LME offers a dramatic gain in statistical power and repeatability of findings, while providing good control of the false positive rate.

Keywords: Linear Mixed Effects models; Longitudinal studies; Mass-univariate analysis; Statistical analysis.

PubMed Disclaimer

Figures

Figure 1
Figure 1
Empirical sensitivity (statistical power) as a function of FDR q-value on AD+HC sub-samples with 2N = 30, randomly drawn from the complete ADNI data (800 random sub-samples). Sensitivity is quantified as the fraction of instances, where the corresponding statistical method detected some group difference at a given FDR q-value. X-Slope: vertex-wise cross-subject analysis of cortical thinning rates estimated by fitting a line to serial measurements; V-LME: vertex-wise application of the LME approach to longitudinal thickness data; ST-LME: the proposed spatiotemporal LME modeling method applied to longitudinal thickness data.
Figure 2
Figure 2
Statistical power as a function of sample size (2N) with FDR q-value = 0.05. See caption of Figure 1 for further details.
Figure 3
Figure 3
Repeatability quantified as the agreement (area overlap) of the detected regions between two independent samples (400 independent AD+HC sample pairs of size 2N=30) as a function of FDR q-value. Error bars show standard error of the mean. See caption of Figure 1 for legend.
Figure 4
Figure 4
Repeatability quantified as the agreement (area overlap) of the detected regions between two independent samples as a function of sample size (with FDR q-value = 0.05). See caption of Figure 3 for further details.
Figure 5
Figure 5
Uncorrected statistical significance –negative log10(p-value)- maps comparing longitudinal cortical thinning rates between HC (N=210) and AD (N=188) subjects (from the entire ADNI sample) visualized on the pial surface of the FreeSurfer template (fsaverage): (A) ST-LME method and (B) X-Slope method. The left hemisphere is shown on the left, and the right hemisphere is on the right. Vertices that have an uncorrected p-value less than 0.05 are shown in color. The odd-numbered rows show the lateral, superior, and anterior views. The even-numbered rows show the medial, inferior, and posterior views. Colorbar shows the corresponding significance value.
Figure 6
Figure 6
Cortical regions exhibiting a statistically significant difference in longitudinal thinning between HC and AD subjects (in red) on the entire ADNI sample. These maps were derived by thresholding the values shown in Figure 5 with an FDR correction at q=0.05. (A) ST-LME method and (B) X-Slope method. ST-LME reveals a much wider extent of significant thinning in AD.
Figure 7
Figure 7
Uncorrected statistical significance –negative log10(p-value)- maps comparing longitudinal cortical thinning rates between stable MCI (N=227) and converter MCI (N=166) subjects (from the entire ADNI sample) visualized on the pial surface of the FreeSurfer template (fsaverage): (A) ST-LME method and (B) X-Slope method. See caption of Figure 5 for further details.
Figure 8
Figure 8
Cortical regions exhibiting a statistically significant difference in longitudinal thinning between stable and converter MCI subjects (in red). These maps were derived by thresholding the values shown in Figure 7 with an FDR correction at q=0.05. (A) ST-LME method and (B) X-Slope method. ST-LME reveals a much more dramatic extent of significant thinning differences between two groups.

References

    1. Asami T, Bouix S, Whitford TJ, Shenton ME, Salisbury DF, McCarley RW. Longitudinal loss of gray matter volume in patients with first-episode schizophrenia: DARTEL automated analysis and ROI validation. Neuroimage 2011 - PMC - PubMed
    1. Benjamini Y, Krieger AM, Yekutieli D. Adaptive linear step-up procedures that control the false discovery rate. Biometrika. 2006;93:491–507.
    1. Bernal-Rusiel JL, Atienza M, Cantero JL. Determining the optimal level of smoothing in cortical thickness analysis: A hierarchical approach based on sequential statistical thresholding. Neuroimage. 2010;52:158–171. - PubMed
    1. Bernal-Rusiel JL, Greve DN, Reuter M, Fischl B, Sabuncu MR. Statistical analysis of longitudinal neuroimage data with Linear Mixed Effects models. Neuroimage. 2012:249–260. - PMC - PubMed
    1. Blockx I, Van Camp N, Verhoye M, Boisgard R, Dubois A, Jego B, Jonckers E, Raber K, Siquier K, Kuhnast B, Dolle F, Nguyen HP, Von Horsten S, Tavitian B, Van der Linden A. Genotype specific age related changes in a transgenic rat model of Huntington’s disease. Neuroimage. 2011;58:1006–1016. - PubMed

Publication types

LinkOut - more resources