我试图可视化以下简单线性模型的二次(而不是线性(关系：
Top_model <- lm(Response~Predictor^2,data=data)
我还试图可视化这种二次关系的残差和置信区间。任何人都可以帮助提供代码建议，允许我在 predictorEffects 中执行此操作？下面的代码导致模型的线性表示，但我寻求二次表示。我需要以某种方式转换 x 轴吗？
plot(predictorEffects(Top_model, ~ Predictor))
非常感谢。
您可以使用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