绘制变量和估计变量之间的GAM关系

Question

我正在运行这样的GAM模型：

library(mgcv)
set.seed(2) ## simulate some data... 
dat <- gamSim(1,n=400,dist="normal",scale=2)
b <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat)

我知道我可以如此策划结果：

plot(b,pages=1)
plot(b, select = 1)
plot(b, select = 2)

但我正在寻找一个在x轴上显示x0，在y轴上显示y，所有数据点以及显示模型b的估计线的图。

为了得到这样的东西：