绘制线性和非线性多元回归线的附加变量图

Question

我想绘制包含线性和非线性关系的多元回归模型的结果。我希望每个预测器有一个绘图，条件是其他预测变量。

以下是我的示例数据：

library(car)
library(broom)
library(jtools)


## fit model: 
model <- lm(mpg ~ poly(disp,3) + wt + carb + factor(gear), data = mtcars)
summary(model)

功能＆＃34; avPlots＆＃34;创建附加变量图，这是一个很好的起点：

avPlots(model)

有没有办法为每个预测器（类似于avPlots）生成单个附加变量图，它还绘制了曲线关系？可能是ggplot的解决方案吗？

我的解决方案就是这样，但这似乎不起作用：

sub <- augment(model, mtcars) # use augment from broom to get fitted values
plot(sub$wt, sub$.fitted) # 
abline(sub$.fitted, sub$wt, col="blue", lwd=3)

回归线不符合要点。此外，这不允许绘制非线性关系。

同样，＆＃34; effect_plot＆＃34;绘制单个添加变量图。但是，这只能在没有多项式关系的情况下起作用：

model <- lm(mpg ~  wt + carb + factor(gear), data = mtcars) # model without poly
summary(model)
effect_plot(model, pred=wt, plot.points = TRUE)

感谢您的帮助!!!