r - 比较r中的两组数据时如何绘制置信区间？

Question

我正在使用 qqplot 比较两组生成的数据，它们都适合伽马。如何将 95% 的置信范围添加到图中？

x <- rgamma(100,shape=1,scale=1)
y <- rgamma(100,shape=1,scale=1.2)
library(car)
myqq <- function(x, y) {
rg <- range(x, y, na.rm=T)
qqplot(x, y, xlim=rg, ylim=rg)
}
myqq(x, y)
abline(0,1, col="red")