Question

我正在尝试使用函数zeroinfl（）设置零膨胀模型的置信区间

如果我使用函数

从线性模型或GLM计算它们

predict(glm, newdata, type = "response", se.fit = TRUE)

它工作正常并返回$ fit，$ se.fit，$ df和$ residual.scale的值（我不知道最后两个输出代表什么 - 对此的澄清也非常受欢迎）。< / p>

然而，当我用零膨胀模型替换glm时，我只得到$ fit的值。

这是一个例子：（我不知道如何使用set.seed - 所以不要责怪我）

set.seed(123)
a <- data.frame(participant = c(1:10),
                activity = c(round(abs(rnorm(10)))), 
                METmin = c(round(abs(rnorm(10)))),
                Var4 = c(round(abs(rnorm(10)))),
                Var5 = c(round(abs(rnorm(10)))))

b <- data.frame(participant = c(1:10),
                     activity = c(round(abs(rnorm(10)))), 
                     METmin = c(round(abs(rnorm(10)))),
                     Var4 = c(round(abs(rnorm(10)))),
                     Var5 = c(round(abs(rnorm(10)))))
identical(a,b) #FALSE
model <- zeroinfl(METmin~activity + Var4 + Var5, data = a)
predict(model, newdata = b, type = "response", se.fit = T)

有谁知道这里出了什么问题以及如何计算这些标准错误？

非常感谢你。

Answer 1

不幸的是，predict()（和zeroinfl）对象的hurdle方法目前没有se.fit参数。如果你想探索这个，你要么必须使用蒙特卡罗方法（通过从系数估计的分布中提取），或者看看包lsmeans中的zeroinfl方法是否可以用于你的目的

R：预测（）为零膨胀模型不返回se.fit

1 个答案: