在Rstudio中将2阶多项式项rm * rm添加到模型lm（公式= medv〜克里姆+印度斯+ rm + dis +克里姆*印度斯）

Question

我的模特是

lm(formula = medv ~ crim + indus + rm + dis + crim * indus)

我必须将2阶多项式项rm*rm添加到上述模型中。我们如何在RStudio中做到这一点？

Answer 1

这里有两个选择。

lm(formula = medv ~ crim + indus + poly(rm, 2, raw = TRUE) + dis + crim * indus)
lm(formula = medv ~ crim + indus + rm + dis + crim * indus + I(rm^2))

首先使用poly(..., raw = TRUE)，其中raw = TRUE确保可以像通常一样解释系数。后者使用as as operator I(...)，它在公式上下文中使用它之前先评估其中的内容。请注意，x^2不会转换为R中的普通多项式。