Concordance between the fast and efficient mixed-effects algorithm (FEMA) and conventional mixed-effects modeling for whole-brain connectivity in longitudinal chronic pain
- 1Chronic Pain and Fatigue Research Center, Department of Anesthesiology, University of Michigan Medical School
Abstract
Researchers increasingly have access to neuroimaging data with repeated measures within subjects. Linear mixed-effects modeling offers a valuable way to analyze such data, but its application remains limited in the neuroimaging literature due to high computational cost and a lack of mainstream analysis packages. Here we benchmark a recently published algorithm, the Fast and Efficient Mixed-Effects Algorithm (FEMA), against a conventional implementation (MATLAB's fitlme, restricted maximum likelihood) across 94,830 functional connectivity edges. In contrast to the Adolescent Brain Cognitive Development (ABCD) Study Dataset used in the development and testing of FEMA, the dataset used in the current analysis has a smaller sample size (n = 378) and comprises individuals with chronic pelvic pain scanned at up to four visits over three years. Despite these differences, the two implementations produced highly concordant results: fixed-effect estimates and test statistics correlated at r [≥] 0.999 (Lin's concordance correlation coefficient [≥] 0.999), and the two approaches reached the same statistical conclusion for 99.98% of edges. FEMA completed each analysis roughly 20 times faster. Agreement was weaker for the estimated variance components, where FEMA attributed systematically less variance to the participant random effect, yielding a slightly lower intraclass correlation in approximately 82% of edges. We further show that most residual disagreement reflected FEMA's default variance-grid resolution rather than the estimator itself, and could be eliminated at negligible computational cost. These findings support the use of FEMA for connectome-wide analysis of repeated-measures neuroimaging data.
Related
- Case report: Progressive organ-threatening IgG4-related disease despite conventional immunosuppression: lessons from a 5-year longitudinal follow-up
- ITHUIM: An efficient algorithm for mining inter-transaction high utility itemsets
- Efficient Deep Learning Models for Predicting Individualized Task Activation From Resting‐State Functional Connectivity
- GCOM-C/SGLI ground-area-based chlorophyll content estimation: a computationally efficient algorithm free from LAI and vegetation classification
- Rising contraceptive use, falling couple concordance: A mixed-methods analysis of the concordance paradox in Nigeria
BibTeX
@misc{mclain2026concordance,
title = {{Concordance between the fast and efficient mixed-effects algorithm (FEMA) and conventional mixed-effects modeling for whole-brain connectivity in longitudinal chronic pain}},
author = {N. McLain and C. Kaplan and S. Harte and A. Schrepf},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.20.753017},
url = {https://www.biorxiv.org/content/10.64898/2026.09.20.753017v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Concordance between the fast and efficient mixed-effects algorithm (FEMA) and conventional mixed-effects modeling for whole-brain connectivity in longitudinal chronic pain
AU - N. McLain
AU - C. Kaplan
AU - S. Harte
AU - A. Schrepf
PY - 2026
DA - 2026/09/20
PB - bioRxiv
DO - 10.64898/2026.09.20.753017
UR - https://www.biorxiv.org/content/10.64898/2026.09.20.753017v1
AB - Researchers increasingly have access to neuroimaging data with repeated measures within subjects. Linear mixed-effects modeling offers a valuable way to analyze such data, but its application remains limited in the neuroimaging literature due to high computational cost and a lack of mainstream analysis packages. Here we benchmark a recently published algorithm, the Fast and Efficient Mixed-Effects Algorithm (FEMA), against a conventional implementation (MATLAB's fitlme, restricted maximum likelihood) across 94,830 functional connectivity edges. In contrast to the Adolescent Brain Cognitive Development (ABCD) Study Dataset used in the development and testing of FEMA, the dataset used in the current analysis has a smaller sample size (n = 378) and comprises individuals with chronic pelvic pain scanned at up to four visits over three years. Despite these differences, the two implementations produced highly concordant results: fixed-effect estimates and test statistics correlated at r [≥] 0.999 (Lin's concordance correlation coefficient [≥] 0.999), and the two approaches reached the same statistical conclusion for 99.98% of edges. FEMA completed each analysis roughly 20 times faster. Agreement was weaker for the estimated variance components, where FEMA attributed systematically less variance to the participant random effect, yielding a slightly lower intraclass correlation in approximately 82% of edges. We further show that most residual disagreement reflected FEMA's default variance-grid resolution rather than the estimator itself, and could be eliminated at negligible computational cost. These findings support the use of FEMA for connectome-wide analysis of repeated-measures neuroimaging data.
ER -