Peak Detection as Multiple Testing

This paper considers the problem of detecting equal-shaped non-overlapping unimodal peaks in the presence of Gaussian ergodic stationary noise, where the number, location and heights of the peaks are unknown. A multiple testing approach is proposed in which, after kernel smoothing, the presence of a peak is tested at each observed local maximum. The procedure provides strong control of the family wise error rate and the false discovery rate asymptotically as both the signal-to-noise ratio (SNR) and the search space get large, where the search space may grow exponentially as a function of SNR. Simulations assuming a Gaussian peak shape and a Gaussian autocorrelation function show that desired error levels are achieved for relatively low SNR and are robust to partial peak overlap. Simulations also show that detection power is maximized when the smoothing bandwidth is close to the bandwidth of the signal peaks, akin to the well-known matched filter theorem in signal processing. The procedure is illustrated in an analysis of electrical recordings of neuronal cell activity.

[1]  Hui Zhang,et al.  Cluster mass inference via random field theory , 2009, NeuroImage.

[2]  Karl J. Friston,et al.  False discovery rate revisited: FDR and topological inference using Gaussian random fields , 2009, NeuroImage.

[3]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[4]  Chi-Sang Poon,et al.  Analysis of First-Derivative Based QRS Detection Algorithms , 2008, IEEE Transactions on Biomedical Engineering.

[5]  Marvin K. Simon,et al.  Digital Communication Techniques: Signal Design and Detection , 2008 .

[6]  David C Christiani,et al.  Biomarker discovery for arsenic exposure using functional data. Analysis and feature learning of mass spectrometry proteomic data. , 2008, Journal of proteome research.

[7]  Y. Benjamini,et al.  False Discovery Rates for Spatial Signals , 2007 .

[8]  Christopher R. Genovese,et al.  Scan clustering: A false discovery approach , 2007 .

[9]  R. Adler,et al.  Random Fields and Geometry , 2007 .

[10]  Jeffrey S. Morris,et al.  Feature extraction and quantification for mass spectrometry in biomedical applications using the mean spectrum , 2005, Bioinform..

[11]  R. Tibshirani,et al.  Sparsity and smoothness via the fused lasso , 2005 .

[12]  Christopher J. Miller,et al.  Spike Hunting in Galaxy Spectra , 2005 .

[13]  K. Worsley,et al.  Unified univariate and multivariate random field theory , 2004, NeuroImage.

[14]  I. Verdinelli,et al.  False Discovery Control for Random Fields , 2004 .

[15]  T. Speed,et al.  Deconvolution of Sparse Positive Spikes , 2004 .

[16]  J. Potter,et al.  A data-analytic strategy for protein biomarker discovery: profiling of high-dimensional proteomic data for cancer detection. , 2003, Biostatistics.

[17]  Noel A Cressie,et al.  Nonparametric hypothesis testing for a spatial signal , 2002, IEEE Workshop on Statistical Signal Processing, 2003.

[18]  M. Meister,et al.  Fast and Slow Contrast Adaptation in Retinal Circuitry , 2002, Neuron.

[19]  Thomas E. Nichols,et al.  Thresholding of Statistical Maps in Functional Neuroimaging Using the False Discovery Rate , 2002, NeuroImage.

[20]  I. Johnstone,et al.  Adapting to Unknown Smoothness via Wavelet Shrinkage , 1995 .

[21]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[22]  Anthony N. Sinclair,et al.  Recovery of a sparse spike time series by L1 norm deconvolution , 1994, IEEE Trans. Signal Process..

[23]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[24]  Georg Lindgren,et al.  Spectral moment estimation by means of level crossings , 1974 .

[25]  S. Berman.  Stationary and Related Stochastic Processes , 1967 .

[26]  harald Cramer,et al.  Stationary And Related Stochastic Processes , 1967 .