qqplot を使用して、生成された 2 つのデータ セットを両方ともガンマに適合させて比較しています。プロットに 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")