0

我想测试我已经适合一组问题数据的动态线性模型的适用性。我已经使用 R 中 dse 包中的 SS() 函数完成了这项工作。有没有任何方法可以使用可能性和信息测试来测试 R 中模型的拟合度？

出于说明目的，假设我的模型是随机游走。随机游走的理论形式为 X(t) = X(t-1) + e(t)~N(0,1) 状态演化 Y(t) = X(t) + w(t)~N (0,1)。R中的代码定义为：

kalman.filter=dse::SS(F = matrix(1,1,1), 
              Q = matrix(1,1,1),
              H = matrix(1,1,1),
              R = matrix(1,1,1),
              z0 = matrix(0,1,1),
              P0 = matrix(0,1,1)
              )

假设当时的实际观察结果是：

simulate.kalman.filter=simulate(kalman.filter, start = 1, freq = 1, sampleT = 100)

然后假设我们拟合一个名为“test”的模型：

test=l(kalman.filter, simulate.kalman.filter)

如何测试数据（simulate.kalman.filter）与 R 中的模型理论模型的拟合度？我正在寻找诸如可能性和贝叶斯信息准则之类的功能。

4

1 回答 1

0

我已经想出了问题的答案。

执行此操作的函数在 dse 的同一包中称为 informationTests()。它将拟合模型的 AIC、BIC 和负对数似然返回数据。在上面的示例中，这是通过以下方式完成的：

informationTests(test)

请记住，具有较低 BIC 的模型被认为更好。您还可以通过添加第二个模型作为参数来比较两个模型（假设您有第二个模型适合名为 test2 的数据）：

informationTests(test, test2)

这将 AIC、BIC 和可能性相互制表。

于 2017-10-13T15:06:01.743 回答