向线性模型添加多项式项

Question

我的模特是

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

我必须将 2 阶的多项式项 rm*rm 添加到上述模型中。我们如何在 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)，其中raw = TRUE 确保可以像往常一样解释系数。后者使用 as is operator I(...) 在公式上下文中使用它之前评估其中的内容。请注意，x^2 不会在公式上下文中转换为 R 中的正规多项式。