R中的卡尔曼滤波器采用Pykalman(Python)方法

时间:2017-08-23 12:54:53

标签: python r kalman-filter pykalman

下午好!

我在Python中使用Pykalman(https://pykalman.github.io/)创建了一个代码,我在其中创建并且没有为卡尔曼滤波器指定任何参数,只是我观察的维数。初始值自动启动(例如转换矩阵的标识),然后使用EM算法,优化卡尔曼滤波器的参数。完成后,我们最终计算出滤波器预测的值。

我很难用R做同样的事情。我使用“dle”包(https://cran.r-project.org/web/packages/dlmodeler/dlmodeler.pdf)。然后我启动我的模型并使用拟合(y为矩阵,84列和2行):

mod=dlmodeler.build(a0 = c(1.142857143,-0.142857143), P0 = diag(2), P0inf = diag(2), Tt = diag(2), Rt = diag(2), Qt = diag(2), Zt = diag(2), Ht = diag(2), name='test')
dlmodeler.fit(y, model=mod, method='MLE')

然后我收到此错误:

  

dlmodeler.build.function(model)出错:多变量案例尚未实施

您是否有任何想法,这不起作用,我应该在我的方法中修改什么?

感谢。

0 个答案:

没有答案