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

【问题标题】：how to include all polynomial and interaction terms如何包含所有多项式和交互项
【发布时间】：2019-08-07 06:50:03
【问题描述】：

我的数据集中有 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)

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

【问题讨论】：

出于兴趣，您为什么要这样做？

标签： r interaction

【解决方案1】：

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

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

【讨论】：

嗨，非常感谢您提供的代码，但为什么这是一个坏主意？我正在考虑考虑交互作用和二次效应。您是否担心多重共线性问题？
@user11806155 它更多地与过度拟合、模型错误指定和使用多项式有关。