让我们考虑以下模型:
set.seed(42)
y <- rnorm(100)
x <- data.frame('Exp' = rexp(100),'Poiss' = rpois(100,1))
model_1 <- lm(y~., data = x)
我想预测(Exp=1,Poiss=0.5(的model_1
的值。
为此,我使用了代码:fitted(model_1, data.frame(Exp = 1, Poiss =0.5))
。
然而,我得到了100个值,而不是一个预测值。我做错了什么?
您可以使用:
#Code
predict(model_1, data.frame(Exp = 1, Poiss =0.5))
因为fitted()
是这样做的:
fitted是一个从对象中提取拟合值的通用函数由建模函数返回。fitted.values是它的别名。
这就是您获得100个值的原因,因为这是从原始model_1
中提取拟合的数据。