2013-10-11 37 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]

Iがあろう95%の信頼区間をプロットして、どのデータポイントが外部にあるかを確認するのが好きです。これを行う方法に関するヒント?

+0

ですから、理論的な正規分布からあなたのサンプル配布を減算したいですか?あなたがしたいようなサウンドは、通常のdist関数にデータをフィットさせ、 'nls'の出力から信頼性データを取得する' 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) 
関連する問題