r - 如何包含所有多项式和交互项

Question

我的数据集中有 50 个预测变量，我想获得这些预测变量的二阶项和二阶交互项。

以 R base 中的 mtcars 数据集为例：

> mtcars2=mtcars^2
> names(mtcars2)=paste0('sq_',names(mtcars2))
> mtcars3=cbind(mtcars,mtcars2)
> formula=paste0('mpg~',paste0(names(mtcars3[,2:22]),collapse='+'),collapse='+')
> fit=lm(formula,mtcars3)
> summary(fit)

我想获得单个模型中所有预测变量的双向交互和多项式项。非常感谢你。

标签： rinteraction

Answer 1

这几乎肯定是个坏主意。无论如何，在删除分类变量后，使用 mtcars 数据集

mtcars=mtcars[,1:7]
lm(as.formula(paste0("mpg~.^2+",paste0("poly(",colnames(mtcars),",2)",collapse="+")))
,data=mtcars)