1

下午好!

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

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

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) 中的错误:尚未实现多变量案例

你有什么想法,这不起作用,我应该在我的方法中修改什么?

谢谢。

4

0 回答 0