R 中的 GLM：具有 95% 置信区间的正态 Q-Q 图

Question

我想要一个函数或包来绘制具有 95% 置信区间的正态 Q-Q 图，但我找不到 GLM，只有 GAM 模型和包 car 中的响应变量。在我的例子中：

#Data set example
p <- read.csv("https://raw.githubusercontent.com/Leprechault/PEN-533/master/bradysia-greenhouse.csv")

#Quasi Poisson GLM
m1 <- glm(bradysia ~ area + mes, family="quasipoisson", data=p)

#Normal Q-Q Plot 
plot(m1, which = 2)

#Normal Q-Q Plot with the confidence interval

我想要像上面那样的 GLM 模型图，可能吗？有什么功能或者包吗？

Answer 1

有两种类型的置信区间：1.confidence 预测区间（您在图 2 中显示的内容）和 2. 回归置信区间。要拟合回归到 glm 模型残差的置信区间，您可以使用 car 包 like

library(car)
qqPlot(m1$resid, ylab="Residuals", xlab="Theoretical Quantiles")

R 中的 GLM：具有 95% 置信区间的正态 Q-Q 图

GLM in R: Normal Q-Q Plot with the 95% confidence interval

r

confidence-interval

glm