Simultaneous Estimation of Noise Variance and Number of Peaks in Bayesian Spectral Deconvolution
This addresses a domain-specific issue for researchers analyzing physical and chemical spectra, offering an incremental improvement over existing methods.
The paper tackles the problem of heuristic peak identification in noisy complex spectra by proposing a Bayesian framework that simultaneously estimates noise variance and number of peaks, preventing overfitting and improving parameter estimation precision.
The heuristic identification of peaks from noisy complex spectra often leads to misunderstanding of the physical and chemical properties of matter. In this paper, we propose a framework based on Bayesian inference, which enables us to separate multipeak spectra into single peaks statistically and consists of two steps. The first step is estimating both the noise variance and the number of peaks as hyperparameters based on Bayes free energy, which generally is not analytically tractable. The second step is fitting the parameters of each peak function to the given spectrum by calculating the posterior density, which has a problem of local minima and saddles since multipeak models are nonlinear and hierarchical. Our framework enables the escape from local minima or saddles by using the exchange Monte Carlo method and calculates Bayes free energy via the multiple histogram method. We discuss a simulation demonstrating how efficient our framework is and show that estimating both the noise variance and the number of peaks prevents overfitting, overpenalizing, and misunderstanding the precision of parameter estimation.