Development of spectral decomposition based on Bayesian information criterion with estimation of confidence interval

ABSTRACT We develop an automatic peak fitting algorithm using the Bayesian information criterion (BIC) fitting method with confidence-interval estimation in spectral decomposition. First, spectral decomposition is carried out by adopting the Bayesian exchange Monte Carlo method for various artificial spectral data, and the confidence interval of fitting parameters is evaluated. From the results, an approximated model formula that expresses the confidence interval of parameters and the relationship between the peak-to-peak distance and the signal-to-noise ratio is derived. Next, for real spectral data, we compare the confidence interval of each peak parameter obtained using the Bayesian exchange Monte Carlo method with the confidence interval obtained from the BIC-fitting with the model selection function and the proposed approximated formula. We thus confirm that the parameter confidence intervals obtained using the two methods agree well. It is therefore possible to not only simply estimate the appropriate number of peaks by BIC-fitting but also obtain the confidence interval of fitting parameters. Graphical abstract

[1]  Hideki Yoshikawa,et al.  Automated information compression of XPS spectrum using information criteria , 2020 .

[2]  Kenji Nagata,et al.  Spectrum adapted expectation-maximization algorithm for high-throughput peak shift analysis , 2019, Science and technology of advanced materials.

[3]  Kazuhiro Yoshihara.  The Introduction of Common Data Processing System Version 12 , 2017 .

[4]  Masato Okada,et al.  Bayesian spectral deconvolution with the exchange Monte Carlo method , 2012, Neural Networks.

[5]  Robert Langer,et al.  High throughput methods applied in biomaterial development and discovery. , 2010, Biomaterials.

[6]  Sumio Watanabe,et al.  Asymptotic behavior of exchange ratio in exchange Monte Carlo method , 2008, Neural Networks.

[7]  R. Hesse,et al.  Product or sum: comparative tests of Voigt, and product or sum of Gaussian and Lorentzian functions in the fitting of synthetic Voigt‐based X‐ray photoelectron spectra , 2007 .

[8]  T. Ishikawa,et al.  Development of hard X-ray photoelectron spectroscopy at BL29XU in SPring-8 , 2005 .

[9]  T. Ishikawa,et al.  High resolution-high energy x-ray photoelectron spectroscopy using third-generation synchrotron radiation source, and its application to Si-high k insulator systems , 2003 .

[10]  Sumio Watanabe.  Algebraic geometrical methods for hierarchical learning machines , 2001, Neural Networks.

[11]  Sumio Watanabe,et al.  Algebraic Analysis for Nonidentifiable Learning Machines , 2001, Neural Computation.

[12]  Thomas F. Coleman,et al.  A Subspace, Interior, and Conjugate Gradient Method for Large-Scale Bound-Constrained Minimization Problems , 1999, SIAM J. Sci. Comput..

[13]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[14]  Y. Ogata.  A Monte Carlo method for an objective Bayesian procedure , 1990 .

[15]  R. Fletcher.  A modified Marquardt subroutine for non-linear least squares , 1971 .

[16]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[17]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .