Open Access
April 2017 Multiple testing of local maxima for detection of peaks in random fields
Dan Cheng, Armin Schwartzman
Ann. Statist. 45(2): 529-556 (April 2017). DOI: 10.1214/16-AOS1458


A topological multiple testing scheme is presented for detecting peaks in images under stationary ergodic Gaussian noise, where tests are performed at local maxima of the smoothed observed signals. The procedure generalizes the one-dimensional scheme of Schwartzman, Gavrilov and Adler [Ann. Statist. 39 (2011) 3290–3319] to Euclidean domains of arbitrary dimension. Two methods are developed according to two different ways of computing p-values: (i) using the exact distribution of the height of local maxima, available explicitly when the noise field is isotropic [Extremes 18 (2015) 213–240; Expected number and height distribution of critical points of smooth isotropic Gaussian random fields (2015) Preprint]; (ii) using an approximation to the overshoot distribution of local maxima above a pre-threshold, applicable when the exact distribution is unknown, such as when the stationary noise field is nonisotropic [Extremes 18 (2015) 213–240]. The algorithms, combined with the Benjamini–Hochberg procedure for thresholding p-values, provide asymptotic strong control of the False Discovery Rate (FDR) and power consistency, with specific rates, as the search space and signal strength get large. The optimal smoothing bandwidth and optimal pre-threshold are obtained to achieve maximum power. Simulations show that FDR levels are maintained in nonasymptotic conditions. The methods are illustrated in the analysis of functional magnetic resonance images of the brain.


Download Citation

Dan Cheng. Armin Schwartzman. "Multiple testing of local maxima for detection of peaks in random fields." Ann. Statist. 45 (2) 529 - 556, April 2017.


Received: 1 May 2014; Revised: 1 February 2016; Published: April 2017
First available in Project Euclid: 16 May 2017

zbMATH: 1369.62144
MathSciNet: MR3650392
Digital Object Identifier: 10.1214/16-AOS1458

Primary: 62H35
Secondary: 62H15

Keywords: False discovery rate , Gaussian random field , image analysis , kernel smoothing , overshoot distribution , selective inference , topological inference

Rights: Copyright © 2017 Institute of Mathematical Statistics

Vol.45 • No. 2 • April 2017
Back to Top