使用predictorEffects从线性模型可视化二次关系

Question

我试图将以下简单线性模型的二次(而不是线性)关系可视化：

Top_model <- lm(Response~Predictor^2,data=data)

我还试图可视化这种二次关系的残差和置信区间。有没有人能帮我提供一些建议，让我在predictorEffects中实现这些功能？下面的代码得到了模型的线性表示，但我寻求的是二次表示。我需要以某种方式变换x轴吗？

plot(predictorEffects(Top_model, ~ Predictor))

非常感谢。

Answer 1

您可以使用poly()以非线性方式拟合您的模型：

Top_model <- lm(Response ~ poly(Predictor,2), data=data)

也可以通过以下方式扩展模型：

Another_Model <- lm(Response ~ poly(Predictor_1,2) * poly(Predictor_2,2), data=data)

它为您提供了Predicator_1和Predictor_2的所有组合，即

Predictor_1
Predictor_1^2

Predictor_2
Predictor_2^2

Predictor_1 * Predictor_2
Predictor_1^2 * Predictor_2
Predictor_1 * Predictor_2^2
Predictor_1^2 * Predictor_2^2