在qq图中添加置信区间？

Question

有没有办法在qqplot中添加置信区间？

我有一个基因表达值的数据集，我用PCA可视化：

pca1 = prcomp(data, scale. = TRUE)

我现在通过以下方式检查数据与正态分布的分布来寻找异常值：

qqnorm(pca1$x,pch = 20, col = c(rep("red", 73), rep("blue", 33)))

qqline(pca1$x)

这是我的数据：

数据= [2.48 104 4.25 219 0.682 0.302 1.09 0.586 90.7 344 13.8 1.17 305 2.8 79.7 3.18 109 0.932 562 0.958 1.87 0.59 114 391 13.5 1.41 208 2.37 166 3.42]

我现在想绘制95％置信区间来检查哪些数据点在外面。关于如何做到这一点的任何提示？

Answer 1

图书馆car提供的功能qqPlot(...)默认情况下会为普通的qq图添加一个逐点置信度：

library(car)
qqPlot(pca1$x)