假设我拟合了一个使用 x1,x2,..xn 预测 y 的标准线性模型,该模型具有很多二向交互作用,也可能是三向交互作用。然后我有一个单独的数据集,它有 x1,x2,...xn,并且想要创建 f(y)。是否有一种直接且简单的方法可以做到不涉及在 R 中手动定义每个交互?
问问题
1516 次
1 回答
3
mnel 是正确的。
dat<-data.frame(matrix(rnorm(30), 10, 3))
names(dat)<-c("y", "x1", "x2")
ans<-lm(y~x1*x2,dat)
summary(ans) # note three parameters + one intercept
predict(ans, newdata=list(x1=0,x2=0)) #only need x1 and x2 for newdata, not interaction
Newdata 可以包含任意数量的条目。给定先前拟合模型的所有参数,包括交互作用,结果将是新数据的预期结果。
于 2013-03-07T01:52:42.313 回答