R中的GLM：正常的Q-Q图，置信区间为95％

问题描述

我想要一个函数或程序包以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模型图，这可能吗？有功能或包装吗？

解决方法

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

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