A False Discovery Rate Control Method Using a Fully Connected Hidden Markov Random Field for Neuroimaging Data

TL;DR

This paper introduces fcHMRF-LIS, a novel spatial FDR control method for neuroimaging data that models complex spatial dependencies, improves power, stability, and scalability over existing methods, and effectively identifies relevant brain regions.

Contribution

It develops a fully connected hidden Markov random field integrated with LIS-based testing, employing efficient algorithms to enhance FDR control, stability, and computational scalability in neuroimaging analysis.

Findings

Achieves accurate FDR control and lower FNR in simulations

Identifies relevant brain regions in Alzheimer's PET data

Offers computational efficiency suitable for high-resolution data

Abstract

False discovery rate (FDR) control methods are essential for voxel-wise multiple testing in neuroimaging data analysis, where hundreds of thousands or even millions of tests are conducted to detect brain regions associated with disease-related changes. Classical FDR control methods (e.g., BH, q-value, and LocalFDR) assume independence among tests and often lead to high false non-discovery rates (FNR). Although various spatial FDR control methods have been developed to improve power, they still fall short of jointly addressing three major challenges in neuroimaging applications: capturing complex spatial dependencies, maintaining low variability in both false discovery proportion (FDP) and false non-discovery proportion (FNP) across replications, and achieving computational scalability for high-resolution data. To address these challenges, we propose fcHMRF-LIS, a powerful, stable, and scalable spatial FDR control method for voxel-wise multiple testing. It integrates the local index of significance (LIS)-based testing procedure with a novel fully connected hidden Markov random field (fcHMRF) designed to model complex spatial structures using a parsimonious parameterization. We develop an efficient expectation-maximization algorithm incorporating mean-field approximation, the Conditional Random Fields as Recurrent Neural Networks (CRF-RNN) technique, and permutohedral lattice filtering, reducing the time complexity from quadratic to linear in the number of tests. Extensive simulations demonstrate that fcHMRF-LIS achieves accurate FDR control, lower FNR, reduced variability in FDP and FNP, and a higher number of true positives compared to existing methods. Applied to an FDG-PET dataset from the Alzheimer's Disease Neuroimaging Initiative, fcHMRF-LIS identifies neurobiologically relevant brain regions and offers notable advantages in computational efficiency.