1
我有一些数据有一个变量的缺失值,我想能够创建(随机)预测这些可能是什么。这是我的第一个念头:R中线性模型的随机预测
# miss indicates where the observations with missing response are
library(MASS)
model <- glm.nb(data[-miss,4] ~ ., data=data[-miss,-4])
predict(model, newdata=data[miss,-4])
但是,如果我重复的最后一行,一遍又一遍给出了同样的答案 - 这似乎给给出的答复预测平均数据和模型。我想要一个包含方差的随机预测,即从给定模型下的这种预测因子的观测响应分布的随机抽取。
它可能与pred.var参数有关,但我不确定如何使用它。