アルゴリズム (判別分析)
判別分析は観測値のグループ所属が分かっているものをさらにグループに振り分け、トレーニングデータとするために使われます。
/math-687a245ac15fad69d39d2ac8207ff385.png "X_t\") をn個の観測値とp個の変数を持った/math-6857d71951ca137bcda4fe14dccae251.png "n_g") グループに所属しているトレーニングデータとします。/math-fbcc679a0035f6fef5099eb0106f996d.png "\bar{x}_j") はｊグループのサンプル平均の行ベクトルで、/math-f898efd8ac64b6621a2aeed350ce3164.png "n_j\") はｊグループの観測された個数です。ｊグループのグループ内共分散行列は次のように表示されます。
/math-42484c7b3cb9982bcabccb7cfda729e9.png "S_j=\frac{1}{n_j-1}\cdot (X_{t}-\bar{x}_j)^T(X_{t}-\bar{x}_j)")
プールされたグループ内共分散行列は次のようになります。
/math-25624e90423e449a29d89e9cb2ef1335.png "S=\frac{1}{n-n_g}\cdot\sum_{j=1}^{n_g} (X_{t}-\bar{x}_j)^T(X_{t}-\bar{x}_j)")
欠損値はこのリスト形式で除かれています。つまり、欠損値を一つ以上含んでいる観測は分析から除外されています。
グループ内共分散行列の均衡検定
もしトレーニングデータが多変量正規分布になると仮定されているとき、次の尤度比検定統計Gを使い、グループ内共分散行列の均衡検定を行えます。
/math-4f141902b982a67990e391daacc238dd.png "G=C{(n-n_g) \mathrm{log} |S|-\sum_{j=1}^{n_g} (n_j-1) \mathrm{log} |S_j|}")
ここで
/math-7b51445c6fda672aa6cfd4a759c53f7d.png "C=1-\frac{2p^2+3p-1}{6(p+1)(n_g-1)}\cdot(\sum_{j=1}^{n_g} \frac{1}{n_j-1} -\frac{1}{n-n_g})")
大きいｎの値は、Gがおおよそ/math-cfbb03435e8c5628a146d4d6c1697e6c.png "\chi^2\") 変数と自由度 /math-d9da675dbc50d16c00b81e4fdf268f83.png "\frac{1}{2}\cdot p(p+1)(n_g-1)") で分散されていると考えられます。
正準判別分析
正準判別分析はp値において、グループ間の分散がグループ内の分散と違いが最大になる割合の線形結合を探すのに使われます。作成された正準変量はグループ間の分離をするのに使われます。
トレーニングデータから総平均を引いたものをXとし、そのランクをkとします。直交行列QはQR分解(すべての列ランクを行う)またはXからのSVDで計算されます。そして /math-309bc456893777dfc6d546a826f80222.png "Q_X\") はｋ列の最初のQになります。/math-7f2bb5755738fc0dc2bce690ffbd02fb.png "Q_g\") を /math-748eca54826a6c8ac743ce0e8d9de0da.png "n_g-1") によってnとし、直交行列によりグループを定義します。ｋｘ /math-bf9c1d397597752f55b120ef5a7d1a98.png "n_g-1\") の行列Vを 次のようにします。
/math-f1ba4d194a3d3c69ac32a0e4f172c72f.png "V=Q_X^TQ_g")
The SVD of V is:
/math-156accfd15fe0adbcb2d4dd525fb0797.png "V=U_X \triangle U_g^T")
ゼロでない行列 /math-8722439958298b3dfbd60cded8990f4e.png "\triangle") の対角要素はiの正準係数で、正準変量とは /math-b1e513a36c4f9d05731ce051a387fc4e.png "\delta_i\") i=1,2,...,と /math-8c794226d0a6faf42a11d59a9f5c7945.png "l=\mathrm{min}(k, n_g)\") で関係しています。
グループ内の二乗和行列における固有値は：
/math-5b9da385f048bd163e59dc5f058db912.png "\lambda_i=\frac{\delta_i^2}{1-\delta_i^2}")
-
- i より有意に大きくなる次元数を検定するには
/math-4e02cd8f2b56fe09ee5b91cfca052400.png "\Lambda_i=\prod_{j=i+1}^{l} 1/(1+\lambda_j)")
-
- 自由度
/math-556713eae767a14fe21eb016694c8238.png "(k-i)(n_g-1-i)\") を持った 統計が使用されて次のようになります。
/math-ec58bf195d61ef078e1bc352e76becac.png "(n-1-n_g-\frac{1}{2}(k-n_g))\sum_{j=i+1}^{l} \mathrm{log}(1+\lambda_j)\ i=0,1,...,l-1")
-
- 正準変量のため行列Bをロードするには
/math-d3f8540a0d3fb0b603051f1afdac58e8.png "U_X\") から計算されます。正準変量はグループ内分散からプールされた単位を持つようにようにスケールされています。つまり：
-
/math-ea1d8b4a1a9949c04ce897debcb0f3cb.png "B^TSB=I\")
- 固有ベクトルの符号はSVDの結果に特有ではありません。つまり、各列のBはー1によってかけることができます。Originは
/math-dd2ed2b7aef1cc29f38fcf893e186990.png "RB\") 内の各和がプラスになるように正規化しています。RはSのコレスキー因数分解結果です。
- 定数は次のように計算されます。
-
/math-842807b1d59b68029f9e47502204de56.png "C_0=-X_mB\")
-
/math-ad5fcaadac0238a2873254ce2bf39933.png "X_m\") は変数の平均を表す行ベクトルです。
-
-
-
/math-f83a67f933bd9989a9c82a155597d0f4.png "D=S_aB\")
-
/math-720199a1df04989fca46b40db9dc3ffa.png "S_a\") は対角行列で対角要素はグループ内の共分散行列Sの対角要素の平方根を取ったものになります。
-
-
-
/math-3d43ba2f89c50ad3d6e41f32eaad76dc.png "C=S_a^{-1}SB\")
-
-
-
/math-ae15be586008ad972cdb2e54d0a8c6c4.png "M_j=C_0+\bar{x}_jB\")
/math-163ee00f28eda2300d069cc196469322.png "M_j\") と /math-9a256199940d00c2536a5816e8562a68.png "\bar{x}_j\") は順に正準グループ平均とｊ番目のグループ平均の行ベクトルとなっています。
-
-
-
/math-e7722c1abe5d5c070d236c68e07ee5c8.png "A_i=C_0+X_iB\")
-
/math-5c89cbde69c9b5306689d7612cc3a31c.png "A_i\") はi番目の観測、 /math-364c437f9831f6eeb78ffd3643f5b2f0.png "X_i\") 、の正準スコアです。
- i 番目の観測はトレーニングデータでもテストデータでも構いません。
マハラノビス距離
マハラノビス距離は観測点とグループ間の距離を測ります。これには2種類の形があります。j 番目のグループと観測点 /math-22cf6d90476d6b87655fb83d74ee08f5.png "x_i\") の距離は次のようになります。
-
/math-4e610adaa8869efc930751889369fb0d.png "D_{ij}^2=(x_i-\bar{x}_j)S_j^{-1}(x_i-\bar{x}_j)^T")
-
/math-d92d0334e383634096e4117f51874db1.png "D_{ij}^2=(x_i-\bar{x}_j)S^{-1}(x_i-\bar{x}_j)^T")
分類
事前確率
事前確率は他のグループから観測データが集められた可能性のユーザ見解を反映します。Originは2種類の先験的確率をサポートしています。
-
-
-
/math-e1de4455b1c3e1e0a1ab05eda657f3ad.png "\pi_j=1/n_g\")
-
-
-
/math-97ca54b52364dcc3af70b85c99b12763.png "\pi_j=n_j/n\")
- はトレーニングデータのｊ番目のグループ内にある観測データの数です。
事後確率
観測データがj 番目のグループから出ている場合のp値は平均 /math-d02f177eaa99e609210550a8e168aa2f.png "\mu_j\") と共分散行列 /math-46ad204a4582a9cd6ba59c2724056d0b.png "\Sigma_j\") をもつ多変量正規分布に当てはまることを前提としています。観測点 がグループｊのから計測された確立は /math-381ca9d538e0de13d017ea1309e83ddb.png "p(x_i|\mu_j,\Sigma_j)\") となるとき、グループｊに所属する事後確率は次のようになります。
-
/math-55fe41c69c4f2dbb7fe75337bfa93c42.png "q_j=p(j|x_i,\mu_j,\Sigma_j)\propto p(x_i|\mu_j,\Sigma_j)\pi_j")
と のパラメータはテストデータ から推測されています。すると観測点は最も事後確率が高いグループに入れられます。Originは事後確率を計算するのに2通りの計算方法を準備しています。
-
- グループ内の共分散行列は等しいと推測されています。
-
-
-
/math-78b58f6de3a01695068779e99f80e92c.png "\mathrm{log}(q_j)=-\frac{1}{2}D_{ij}^2+\mathrm{log}(\pi_j)+c_0")
-
-
/math-73850f915f6390176ce13cfdd422f1a3.png "D_{ij}^2") はi 番目の観測点とプールされたグループ内共分散行列のマハラノビス距離で、 /math-ecd67c43f1d5cc8aa1cc153e7b57fd71.png "c_0\") は定数です。
-
- グループ内の共分散行列は等しくないと推測されています。
-
-
-
/math-97bb999f1c63118e1c15d0cd99f3ecd9.png "\mathrm{log}(q_j)=-\frac{1}{2}D_{ij}^2+\mathrm{log}(\pi_j)-\frac{1}{2}\mathrm{log}|S_j|+c_0")
-
- はi 番目の観測点とプールされたグループ内共分散行列のマハラノビス距離で、 は定数です。
/math-f680dc664fc7f822ae305c4d1586956d.png "q_j\") は次のように正規化され、 は正規化したものから算出されます。
-
/math-766af91331121006a73b3d401c5e661e.png "\sum_{j=1}^{n_g} q_j=1")
異型性インデックス
異型性インデックス/math-cc421199a07c0dc7f86b04779cf2a388.png "I_j(x_i)\") は観測点がi番目の観測データよりもグループjにより典型的に当てはまる確立を示しています。すべてのグループで１に近い場合、それはトレーニングデータに入っているグループから来ていない可能性が高いことを意味しています。異型性インデックスは次のように計算されます。
-
/math-732acf44e40ad9cac61088dcd141df05.png "I_j(x_i)=P(B\le z:\frac{1}{2}p,\frac{1}{2}(n_j-d))")
/math-8c5084405b04323ffd7f2875ca6a833b.png "P(B\le \beta:\ a, b)") はグループ内共分散行列が等しいもののベータ分布の下限破棄の可能性を示しています。
-
/math-cd6fb8799495aa07b0c4c32cb5b97836.png "z=D_{ij}^2/(D_{ij}^2+(n-n_g)(n_j-1)/n_j)")
グループ内共分散行列が等しくないものに関しては次のようになります。
-
/math-45a4cc9d19289a7db844343cc7755561.png "z=D_{ij}^2/(D_{ij}^2+(n_j^2-1)/n_j)")
線形判別関数係数
線形判別関数(またはフィッシャーの線形判別関数)は次のように計算されます。
-
-
-
/math-bba90644f0ee4526a124e9a060efad5f.png "b_j=S^{-1}\bar{x}_j^T")
/math-e10b85a99fd3e0319c1a12b37038bee1.png "b_j\") はpの大きさの列ベクトルです。
-
-
-
/math-b49fb8d88d19012d29f58a5ec2d3145e.png "a_j=\bar{x}_jb_j")
トレーニングデータを分類する
トレーニングデータ内の観測データは事後確率によって分類されることができます。つまり事後確率が一番高いグループに入れられることになります。各グループからのマハラノビス距離の二乗と各グループからの異型性インデックスもそれぞれ計算されます。
トレーニングデータの分類結果は与えられたグループメンバーシップと予想されるグループの結果を比較した状態でまとめられています。分類失敗された観測データに事前確率で重み付けされた分類失敗エラー率はパーセンテージ(％)で計算されます。
-
/math-23a8991fe54e13492cec65429649bf96.png "E=\sum_{j=1}^{n_g} e_j\pi_j")
/math-c6c1265e2ebb2dd779b29fea40842909.png "e_j\") はｊグループに分類失敗の観測点があるパーセントです。
トレーニングデータをクロス検証する
トレーニングデータを分類すると同じ手順に従ってください。ただし、トレーニングデータ内で観測点の所属を明らかにするとき、グループ内共分散行列またはプールされたグループ内共分散行列を計算するときに除外されたデータは使用できません。
テストデータを分類する
グループ内共分散行列とプールされたグループ内共分散行列はトレーニングデータから算出されます。テストデータ内の観測データは事後確率によって分類されることができます。つまり事後確率が一番高いグループに入れられることになります。
|