An MCMC algorithm for parameter estimation in signals with hidden intermittent instability. (English) Zbl 1350.60071
Summary: Prediction of extreme events is a highly important and challenging problem in science, engineering, finance, and many other areas. The observed extreme events in these areas are often associated with complex nonlinear dynamics with intermittent instability. However, due to lack of resolution or incomplete knowledge of the dynamics of nature, these instabilities are typically hidden. To describe nature with hidden instability, a stochastic parameterized model is used as the low-order reduced model. Bayesian inference incorporating data augmentation, regarding the missing path of the hidden processes as the augmented variables, is adopted in a Markov chain Monte Carlo (MCMC) algorithm to estimate the parameters in this reduced model from the partially observed signal. Howerver, direct application of this algorithm leads to an extremely low acceptance rate of the missing path. To overcome this shortcoming, an efficient MCMC algorithm which includes a pre-estimation of hidden processes is developed. This algorithm greatly increases the acceptance rate and provides the low-order reduced model with a high skill in capturing the extreme events due to intermittency.
MSC:
60J22 | Computational methods in Markov chains |
60J10 | Markov chains (discrete-time Markov processes on discrete state spaces) |
60G25 | Prediction theory (aspects of stochastic processes) |
62F15 | Bayesian inference |
62M20 | Inference from stochastic processes and prediction |
65C05 | Monte Carlo methods |
65C40 | Numerical analysis or methods applied to Markov chains |
62P12 | Applications of statistics to environmental and related topics |
94A15 | Information theory (general) |