将多项式项添加到线性模型

Question

我的模型是

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

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

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)whereraw = TRUE确保系数可以像平常一样解释。后者使用as is 运算符 I(...)，在公式上下文中使用它之前先评估其中的内容。请注意，在公式上下文中，x^2它不会转换为R中的正规多项式。