R 如何仅从gam.check获取绘图

在应用

mgcv

包中的

gam.check

时，

R

生成一些残差图和基准尺寸输出。是否有一种方法只生成绘图而不生成打印输出

library(mgcv)
set.seed(0)
dat <- gamSim(1,n=200)
b   <- gam(y~s(x0)+s(x1)+s(x2)+s(x3), data=dat)
plot(b, pages=1)
gam.check(b, pch=19, cex=.3)

库（mgcv）
种子集（0）
dat有四个图，从左上角开始，向下和穿过我们有：

残差的QQ图
残差直方图
残差图与线性预测图
观察值与拟合值的曲线图
在下面的代码中，我假设
b
包含您的已安装型号，如您的示例所示。首先，我们需要一些东西

type <- "deviance"  ## "pearson" & "response" are other valid choices
resid <- residuals(b, type = type)
linpred <- napredict(b$na.action, b$linear.predictors)
observed.y <- napredict(b$na.action, b$y)

在R提示下

每个地块的制作如下所示：

QQ图
这是通过
qq.gam（）
生成的：

残差直方图
这是使用

hist(resid, xlab = "Residuals", main = "Histogram of residuals")

plot(linpred, resid, main = "Resids vs. linear pred.", 
     xlab = "linear predictor", ylab = "residuals")

plot(fitted(b), observed.y, xlab = "Fitted Values", 
     ylab = "Response", main = "Response vs. Fitted Values")

残差与线性预测
这是使用

hist(resid, xlab = "Residuals", main = "Histogram of residuals")

plot(linpred, resid, main = "Resids vs. linear pred.", 
     xlab = "linear predictor", ylab = "residuals")

plot(fitted(b), observed.y, xlab = "Fitted Values", 
     ylab = "Response", main = "Response vs. Fitted Values")

观测值与拟合值
这是使用

hist(resid, xlab = "Residuals", main = "Histogram of residuals")

plot(linpred, resid, main = "Resids vs. linear pred.", 
     xlab = "linear predictor", ylab = "residuals")

plot(fitted(b), observed.y, xlab = "Fitted Values", 
     ylab = "Response", main = "Response vs. Fitted Values")

所有三个涉及残差的地块？还是所有的情节都跳过印刷品？绘图很容易制作，但如果你只想要一个，我不想展示如何获得所有的绘图。@GavinSimpson对不起，我想要的是所有的绘图（4），而不是印刷的材料，我怀疑这最终会被移植到，因为你的问题中没有任何内在的统计信息，这是正确的R咒语。是的，在问了这个问题后，我也这么想