假设我适合使用x1,x2,... xn预测y的标准线性模型,该模型具有很多两个且可能是三向交互。然后我有一个单独的数据集,其中包含x1,x2,... xn,并且想要创建f(y)。是否有一种直接且简单的方法可以做到这一点,而不涉及在R中手动定义每个交互?
答案 0 :(得分: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可以包含任意数量的条目。在给出先前拟合模型的所有参数(包括相互作用)的情况下,结果将是新数据的预期结果。