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

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)

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

Answer 1

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

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

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

how to include all polynomial and interaction terms

interaction

r