3

カイ二乗分布に対してプロットされた多変量データのマハラノビス距離の多変量qqplotを実行すると、付随するqqlineは切片0と傾き1の線になると予想されます。ただし、次のコードを実行すると、次のようになります。

scores<-matrix(rnorm(100*3),nc=3)
mah_scores = mahalanobis(scores,center=colMeans(scores),cov=cov(scores))
chi_scores = qchisq(ppoints(nrow(scores)),df=3)
qqplot(x=chi_scores,y=mah_scores,asp=1)
abline(a=0,b=1)
qqline(mah_scores, distribution = function(p) qchisq(p,df = 3),col="red")

次の図が表示されます。

ablineとqqlineを使用したqqplot

qqline（赤）は切片0と傾き1（黒）の線と同じであると予想しました。これらの2行が一致しない理由を誰かが私に説明できますか？

（私はRバージョン2.15.3（2013-03-01）を実行しています）

4

1 に答える 1

5

デフォルトでqqlineは、1 番目と 3 番目の四分位数を通る線を描画します。ヘルプ ファイル ( ?qqline) を参照してください。

#highlight first and third quartiles
points(quantile(chi_scores,c(.25,.75)),
       quantile(mah_scores,c(.25,.75)),col="blue",cex=2,bg="blue",pch=21)

ここに画像の説明を入力

于 2013-03-16T13:31:18.840 に答える