Question

我正在使用KFAS来拟合表格的动态逻辑模型;

Y 1 =β_t* X +ε

其中y是长度为n的向量，β是长度为p的向量，x为n * p的矩阵。

β_t=β_（T-1）+η

因此，回归参数随时间变化，并作为潜在变量由过滤器估算。

如何在R中指定这样的模型？我不能使用MARSS包，因为我需要一个逻辑链接函数，但我一直在尝试使用的KFAS包（它接受二项式分布）的文档很少。

关键问题似乎是我每个时间段都有多个观察结果，这个包可能不支持。关于我的问题的一个例子，请看下面的代码 - 它应该展示十个时间段，每个时间段有2-3个，但是KFAS认为每一行都是一个单独的时间段，其中有22个而不是10个。

library(KFAS)
y = c(1,0,0,0,1,0,1,1,0,1,1,0,1,1,0,1,1,0,1,1,0,1)
i = seq.Date(from = as.Date("2014-01-01"), as.Date("2014-01-10"), length.out = 22)
x = rnorm(n = 22, mean = 1, sd = 2)

a =   model = SSModel(y ~ 
                    SSMregression(~x),
                  distribution = "binomial")

fit = fitSSM(a, inits = c(0,0))

如何在R（KFAS）

0 个答案: