如何自动知道 R 中 stepAIC 中选定的变量？

Question

运行“stepAIC”后，我得到以下结果。

 fit1=lm(y~ x1+x2+x3+x4+x5)
 fit2=stepAIC(fit1)
 coef=fit2$coefficients

 >coef
 >intercept, x1,  x3,  x5
       5      1,    3,   5

我有另一个预测向量 z=(z1,...,z5)。因为它是 for 循环的一部分，所以我想使用 'coef %*% z' 自动计算预测 y。

我可以想出两种方法来实现它：
1.将stepAIC的未选择系数设为0；所以而不是

   coef=c(5,1,3,5)

我有

   coef=c(5,1,0,3,0,5)

1. 找到选定的系数并找到对应的“z_i”。

我不知道如何实现这一点。任何帮助，将不胜感激。提前致谢。

Answer 1

step 方法/设施只返回另一个lm 对象，因此您可以应用任何通用函数，包括predict。

predict(fit2, newdata = a.data.frame)

---

如果最终目标不是预测，而是如您的问题标题所述，请使用attr(terms(fit2), 'term.labels')。