Multiple testing of local maxima for detection of peaks in random fields

DOI10.1214/16-AOS1458zbMATH Open1369.62144arXiv1405.1400OpenAlexW2963770316WikidataQ90135769 ScholiaQ90135769MaRDI QIDQ2012197FDOQ2012197

Armin Schwartzman, Dan Cheng

Publication date: 28 July 2017

Published in: The Annals of Statistics (Search for Journal in Brave)

Abstract: 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 et al. (2011) 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 (Cheng and Schwartzman, 2014), available explicitly when the noise field is isotropic; (ii) using an approximation to the overshoot distribution of local maxima above a pre-threshold (Cheng and Schwartzman, 2014), applicable when the exact distribution is unknown, such as when the stationary noise field is non-isotropic. 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 non-asymptotic conditions. The methods are illustrated in a nanoscopy image analysis problem of detecting fluorescent molecules against the image background.

Full work available at URL: https://arxiv.org/abs/1405.1400

Recommendations

• Multiple testing of local maxima for detection of peaks in 1D
• Multiple testing of local maxima for detection of peaks on the (celestial) sphere
• Testing for signals with unknown location and scale in a χ²random field, with an application to fMRI
• False Discovery Rate Smoothing
• False Discovery Rates for Spatial Signals

zbMATH Keywords

false discovery rateGaussian random fieldkernel smoothingselective inferenceimage analysistopological inferenceovershoot distribution

Mathematics Subject Classification ID

Random fields; image analysis (62M40) Hypothesis testing in multivariate analysis (62H15) Image analysis in multivariate analysis (62H35)

Cited In (18)

• Smooth Matérn Gaussian random fields: Euler characteristic, expected number and height distribution of critical points
• Local repulsion of planar Gaussian critical points
• On critical points of Gaussian random fields under diffeomorphic transformations
• Expected number and height distribution of critical points of smooth isotropic Gaussian random fields
• False discovery variance reduction in large scale simultaneous hypothesis tests
• False discovery rate for functional data
• Modelling persistence diagrams with planar point processes, and revealing topology with bagplots
• Distribution of the height of local maxima of Gaussian random fields
• Multiple testing of local maxima for detection of peaks on the (celestial) sphere
• Local maximin properties of tests in Gaussian shift experiments
• Mean number and correlation function of critical points of isotropic Gaussian fields and some results on GOE random matrices
• An approximation to peak detection power using Gaussian random field theory
• A multi-scale approach for testing and detecting peaks in time series
• Multiple Testing of Local Extrema for Detection of Change Points
• Semiparametric Bayesian Inference for Local Extrema of Functions in the Presence of Noise
• Detection of spurious maxima through random draw tests and specification tests
• Multidimensional multiscale scanning in exponential families: limit theory and statistical consequences
• dSTEM

This page was built for publication: Multiple testing of local maxima for detection of peaks in random fields

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2012197)