The Annals of Applied Statistics
- Ann. Appl. Stat.
- Volume 5, Number 2B (2011), 1488-1511.
Bayesian nonparametric models for peak identification in MALDI-TOF mass spectroscopy
We present a novel nonparametric Bayesian approach based on Lévy Adaptive Regression Kernels (LARK) to model spectral data arising from MALDI-TOF (Matrix Assisted Laser Desorption Ionization Time-of-Flight) mass spectrometry. This model-based approach provides identification and quantification of proteins through model parameters that are directly interpretable as the number of proteins, mass and abundance of proteins and peak resolution, while having the ability to adapt to unknown smoothness as in wavelet based methods. Informative prior distributions on resolution are key to distinguishing true peaks from background noise and resolving broad peaks into individual peaks for multiple protein species. Posterior distributions are obtained using a reversible jump Markov chain Monte Carlo algorithm and provide inference about the number of peaks (proteins), their masses and abundance. We show through simulation studies that the procedure has desirable true-positive and false-discovery rates. Finally, we illustrate the method on five example spectra: a blank spectrum, a spectrum with only the matrix of a low-molecular-weight substance used to embed target proteins, a spectrum with known proteins, and a single spectrum and average of ten spectra from an individual lung cancer patient.
Ann. Appl. Stat. Volume 5, Number 2B (2011), 1488-1511.
First available in Project Euclid: 13 July 2011
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
House, Leanna L.; Clyde, Merlise A.; Wolpert, Robert L. Bayesian nonparametric models for peak identification in MALDI-TOF mass spectroscopy. Ann. Appl. Stat. 5 (2011), no. 2B, 1488--1511. doi:10.1214/10-AOAS450. https://projecteuclid.org/euclid.aoas/1310562730.
- Supplementary material: Additional results for the simulation study. True positive rates for LARK estimates from the simulation study broken down by peak prevalence and average intensity of peaks across samples.