《判別関数》

【はんべつかんすう (discriminant function) 】

　いくつかの変数(特性)についての測定値が得られている対象に対して, それが属している可能性があるグループが複数考えられるときに, それらの変数の関数を用いて対象の属するグループを判別することにする. このときに用いる関数を判別関数という.

　いくつかの特性の値からグループを判別するから,特性が説明変数であり, グループが(質的)目的変数である.説明変数を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_i(i=1,\ 2,\ \cdots,\ p)\, } , 目的変数を$y\,$で表す.また, $y\,$のとりうる値(グループ名)を$G_{h}(h=1,\ 2,\ \cdots ,\ r)\,$とする.すなわち, $r\,$ 個のグループが考えられているとする. グループの判別には, ${\boldsymbol {x}}(x_{i}(i=1,\ 2,\ \cdots ,\ p)\,$を並べたベクトル)と$G_{h}(h=1,\ 2,\ \cdots ,\ r)\,$の中心(平均)の間の距離$D_{h}({\boldsymbol {x}})\,$を用いる. $G_{h}\,$における平均ベクトル $(x_{i}(i=1,\ 2,\ \cdots ,\ p)\,$ の平均を並べたベクトル) を${\boldsymbol {m}}_{h}\,$, 分散共分散行列の逆行列を $C_{h}\,$ とする. このとき, $D_{h}({\boldsymbol {x}})\,$ は, 次式で計算される.

$D_{h}({\boldsymbol {x}})=({\boldsymbol {x}}{\boldsymbol {m}}_{h})^{\top }C_{h}({\boldsymbol {x}}-{\boldsymbol {m}}_{h})$

グループが正規母集団とみなされ, 分散共分散行列がすべて等しいとき, 上の式で ${\boldsymbol {x}}\,=\,{\boldsymbol {m}}_{k}\,$とおいて得られる距離を, $G_{k}\,$と$G_{h}\,$の間のマハラノビス汎距離という. 平均や分散共分散行列は, 各グループに属していることがわかっている対象についての測定値より計算される. $D_{h}({\boldsymbol {x}})(h=1,\ 2,\ \cdots ,\ r)\,$ の中で,$D_{k}({\boldsymbol {x}})\,$ が最小であれば, この対象は,$G_{k}\,$に属していると判別すればよい. また, どれにも属さないという判別が許される場合は, あらかじめ上限を設定しておいて, $D_{k}({\boldsymbol {x}})\,$ がそれを越えたときは, どれにも属さないと判別すればよい.

　構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle r=2\, } のときは, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathit{\Delta}D_{12}(\boldsymbol{x})=D_1(\boldsymbol{x})-D_2(\boldsymbol{x})\, } を計算して, ${\mathit {\Delta }}D_{12}({\boldsymbol {x}})>0\,$であれば 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle G_2\, } に属し, ${\mathit {\Delta }}D_{12}({\boldsymbol {x}})<0\,$であれば 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle G_1\, } に属すると判別すればよい. 分散共分散行列が等しいとき, すなわち, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle C_1=C_2=C\, } であるとき,

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathit{\Delta}D_{12}(\boldsymbol{x}) =2(\boldsymbol{m}_2-\boldsymbol{m}_1)^{\top} C\boldsymbol{x}-(\boldsymbol{m}_2-\boldsymbol{m}_1)^{\top} C(\boldsymbol{m}_1+\boldsymbol{m}_2) }

と変形できるので, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathit{\Delta}D_{12}(\boldsymbol{x})\, } は, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x_i(i=1,\ 2,\ \cdots,\ p)\, } の線形式になる. したがって, これを(構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle G_1\, } と構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle G_2\, } を判別する)線形判別関数という. 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle r\, } が3以上のときは, 線形判別関数は, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {}_r{\rm C}_2\, } 個できる. なお, 分散共分散行列が等しくないときは, 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathit{\Delta}D_{12}(\boldsymbol{x})\, } は, $x_{i}(i=1,\ 2,\ \cdots ,\ p)\,$ の2次式になる.

---

参考文献

[1] 奥野忠一, 久米均, 芳賀敏郎, 吉澤正, 『多変量解析法(改訂版)』, 日科技連出版, 1981.