在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)