2013-10-11 72 views
6

有沒有辦法將置信區間添加到qqplot?爲qq圖添加置信區間?

我有基因表達值的數據集,我一直在使用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%的置信區間來檢查哪些數據點在外面。有關如何做到這一點的任何提示?

+0

所以,你想從理論正態分佈減去你的樣本分佈?聽起來你想要做的就是使用'nls'來使你的數據適合正常的dist函數,並從'nls'的輸出中獲取置信度數據。 –

+0

如果您提供[最小,可重現的數據集],您更可能收到有用答案(http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example/5963610#5963610)以及您嘗試過的代碼。謝謝! – Henrik

+0

我用一些數據編輯了我的初始文章。 是否可以從qqnorm的輸出中獲取置信度數據? – user2846211

回答

11

car提供了默認情況下增加了一個逐點置信區間,以正常的QQ積函數qqPlot(...)

library(car) 
qqPlot(pca1$x)