The invention described herein details a protocol to improve analysis and peak identification in spectroscopic data. Bayesian methods are used to automatically identify peaks in data sets. After identifying peak shapes, the method tests the hypothesis that a given number of peaks is found within any given data window. If a peak is identified within a given window, then the likelihood function is maximized in order to estimate peak position and amplitude. This process yields a spectrum with high resolution and minimal artifacts. The method described herein is particularly useful for identifying peaks in data sets obtained from spectroscopy.