我正在尝试使用以下状态空间模型。
(1)Kt = K(t-1)* + ε1t
(2)Yt = Kt + βZt + ε2t
其中,t是时间,Yt是可观察变量(在t),Kt
是不可观察的趋势,Zt
是可以解释Yt
的可观察变量矩阵。 ε1
和ε2
是状态空间模型中的常用错误术语。
我需要通过MLE估计以下过程,得到趋势Kt和系数β矩阵。但是我没有找到使用dlm和dlmModReg的方法。
我知道如何在Matlab中估算这个模型(参见下面的链接)。但我没有看到用dlm包指定像这样的模型的方法。是否可以使用dlm包的功能来实现它?
(我看到有一篇未答复的相关帖子标题为“dlm包中的外生变量”)
真诚地感谢任何帮助(即使它意味着使用另一个包裹)!