今回は、状態空間モデルの中でも基礎となる、ローカルレベルモデルをPyStan、PyMC3、Edwardで実装してみました。 コードは以下にも置いておきました。 ※Stanのログ出力も入ってしまって、とても長いので注意 GitHub: https://github.com/Gin04gh/datascience/blob/master/compare_state_space_model/notebook.ipynb Statsmodelsで状態空間モデル 適当に時系列データを作成します。 import numpy as np import matplotlib import matplotlib.pylab as plt %matplotlib inline from tqdm import tqdm y = np.cumsum(np.random.normal(size=100)) plt