相関比を最大にすることによる判別係数の求め方
- 群の数を k,各群のケース数を n1,n2, ... ,nk とする。
- p 個の変数を X1,X2, ... ,Xp として,任意の重み係数 a1, a2, ... , ap を用いて作られる合成変量を Z とする。
- 第 j 群,第 i ケースの合成変量を Zij( j = 1,2, ... ,k; i = 1,2, ... ,nj )とする。
- 全体の平均値を
,第 j 群における平均値を
とすれば,Z の平方和 St は 群内平方和 Sw と 群間平方和 Sb に分解できる。
- Z により各群がよく判別できるということは,相関比 η2 = Sb / St が大きいということに対応するので( あるいは,Sb / Sw の比が大きいことに対応すると考えてもよい ),
相関比 が最大になるように重み係数 a1, a2, ... , ap を決定すればよい。これを,判別係数 と呼ぶ。
- l 群の m 番目のケースの i 番目の変数の測定値を
Ximl,l 群の i 番目の変数の平均値を
,i 番目の変数の全体の平均値を
としたとき,観察値の群内および群間平方和・積和行列の要素を,
とすると,
…… ( 1 )
- ( 1 )式を,ai で偏微分して 0 とおき行列を用いて表すと,( 2 )式のようになる。
…… ( 2 )
- ( 2 )式が自明の解以外を持つためには係数行列式が( 3 )式の条件を満たさなければならない。
…… ( 3 )
- W の逆行列を W - 1 とすると,( 4 ) 式となり,θ は行列 W - 1 B の固有値であることがわかる。
…… ( 4 )
- ( 4 )式を満たす固有値は複数個( m = min( p,k - 1 )個 )存在するが,( 1 )式を最大にするのはその内の最大の固有値である。また,係数ベクトル a はその固有値に対応する固有ベクトルである。
- 2 番目以降に大きい固有値とその固有ベクトルも判別関数を構成することができる。
- それぞれの判別関数の寄与率は( 5 )式のようになり,m より少ない
判別関数で十分な判別を行うことができる。すなわち,次元の減少を伴う判別 であるといわれる(正準判別分析)。
…… ( 5 )
注:二群の判別の場合には,次ページ後半に示すように計算は簡単になる。
- 相関比を最大にすることによる判別係数の求め方のページへのリンク