Rで判別分析（２）－分散が等しくない場合：マハラノビスの距離を求める－

内容

グループ１（Ｇ１）、グループ２（Ｇ２）とあって、それぞれのデータがある。
このとき、あるデータＡ＝７０があったとき、Ａは、グループ１に属するか、
グループ２に属するかを見たい。
ただし、前回と違い、このＧ１，Ｇ２の分散は異なる。

---

方法
１．データを作成・読み込む（前回と同じ）

＜＜グループ１、グループ２のデータが１項目しかないとき＞＞

例：ある試験　受験者１０人（うち、５人合格）
　　グループ１（合格者）：100,75,80,90,95
　　グループ２（不合格）：68,65,50,40,55

Rでの入力：
label＜-c(rep(1,5),rep(2,5))
g1＜-c(100,75,80,90,95)
g2＜-c(68,65,50,40,55)


２．各グループの平均と分散を求める

Rでの入力：
xbar1＜-mean(g1)
xbar2＜-mean(g2)
s1＜-var(g1)
s2＜-var(g2)

３．マハラノビスの距離の値を求める
a＜-70
mahalanobis(a,xbar1,s1)
mahalanobis(a,xbar2,s2)

距離が近いほう（数字が小さいほう）のグループに属すると考えられる。