下午好!
我在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)出错:多变量案例尚未实施
您是否有任何想法,这不起作用,我应该在我的方法中修改什么?
感谢。