TY - JOUR
T1 - Time-series filtering for replicated observations via a kernel approximate Bayesian computation
AU - Hasegawa, Takanori
AU - Kojima, Kaname
AU - Kawai, Yosuke
AU - Nagasaki, Masao
N1 - Funding Information:
Manuscript received July 2, 2017; revised March 19, 2018, June 4, 2018, and August 5, 2018; accepted September 14, 2018. Date of publication October 1, 2018; date of current version October 29, 2018. The associate editor coordinating the review of this manuscript and approving it for publication was Dr. Paolo Braca. This work was supported in part by Grant-in-Aid for Young Scientists (Start-up) Grant Number 15H06008 and by MEXT Tohoku Medical Megabank Project. The super-computing resource was provided by Human Genome Center, the Institute of Medical Science, the University of Tokyo. (Corresponding author: Takanori Hasegawa.) T. Hasegawa is with the Health Intelligence Center, The Institute of Medical Science, The University of Tokyo, Tokyo 108-8639, Japan (e-mail:,t-hasegw@ ims.u-tokyo.ac.jp).
Publisher Copyright:
© 2018 IEEE.
PY - 2018/12/1
Y1 - 2018/12/1
N2 - In time-series analysis, state-space models (SSMs) have been widely used to estimate the conditional probability distributions of hidden variables and parameter values, as well as to understand structures that can generate the data. For example, the Kalman filter is used to analytically calculate the conditional probability distribution on linear SSMs in terms of minimizing the variance, and several extensions, such as the unscented Kalman filter and particle filter, have been applied to calculate the approximate distribution on nonlinear SSMs. Recently, the approximate Bayesian computation (ABC) has been applied to such time-series filtering to handle intractable likelihoods; however, it remains problematic with respect to consistently achieving a reduction of the estimation bias, evaluating the validity of the models, and dealing with replicated observations. To address these problems, in this paper, we propose a novel method combined with the kernel ABC to perform filtering, parameter estimation, and the model evaluation in SSMs. Simulation studies show that the proposed method produces inference that is comparable to other ABC methods, with the advantage of not requiring a careful calibration of the ABC threshold. In addition, we evaluate the performance of the model-selection capability using true and competitive models on synthetic data from nonlinear SSMs. Finally, we apply the proposed method to real data in rat circadian oscillations, and demonstrated the usefulness in practical situations.
AB - In time-series analysis, state-space models (SSMs) have been widely used to estimate the conditional probability distributions of hidden variables and parameter values, as well as to understand structures that can generate the data. For example, the Kalman filter is used to analytically calculate the conditional probability distribution on linear SSMs in terms of minimizing the variance, and several extensions, such as the unscented Kalman filter and particle filter, have been applied to calculate the approximate distribution on nonlinear SSMs. Recently, the approximate Bayesian computation (ABC) has been applied to such time-series filtering to handle intractable likelihoods; however, it remains problematic with respect to consistently achieving a reduction of the estimation bias, evaluating the validity of the models, and dealing with replicated observations. To address these problems, in this paper, we propose a novel method combined with the kernel ABC to perform filtering, parameter estimation, and the model evaluation in SSMs. Simulation studies show that the proposed method produces inference that is comparable to other ABC methods, with the advantage of not requiring a careful calibration of the ABC threshold. In addition, we evaluate the performance of the model-selection capability using true and competitive models on synthetic data from nonlinear SSMs. Finally, we apply the proposed method to real data in rat circadian oscillations, and demonstrated the usefulness in practical situations.
KW - approximate Bayesian computation
KW - likelihood-free inference
KW - Nonlinear state-space model
KW - time-series analysis
UR - http://www.scopus.com/inward/record.url?scp=85054351441&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85054351441&partnerID=8YFLogxK
U2 - 10.1109/TSP.2018.2872864
DO - 10.1109/TSP.2018.2872864
M3 - Article
AN - SCOPUS:85054351441
SN - 1053-587X
VL - 66
SP - 6148
EP - 6161
JO - IEEE Transactions on Signal Processing
JF - IEEE Transactions on Signal Processing
IS - 23
M1 - 8478007
ER -