9

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.2 81 4] 1.2 7.4

95% の信頼区間をプロットして、外側にあるデータ ポイントを確認したいと思います。これを行う方法に関するヒントはありますか?

4

1 に答える 1

15

ライブラリは、デフォルトで通常の qq-plot に点ごとの信頼度エンベロープを追加するcar関数を提供します。qqPlot(...)

library(car)
qqPlot(pca1$x)
于 2013-10-11T13:23:48.230 に答える