《判別関数》

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

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

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

　D_h(\mbox{\boldmath $x$})=

 (\mbox{\boldmath $x$}-\mbox{\boldmath $m$}_h)^{\top}
 C_h(\mbox{\boldmath $x$}-\mbox{\boldmath $m$}_h)

グループが正規母集団とみなされ, 分散共分散行列がすべて等しいとき, 上の式で 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle {\boldmath $x$}\,=\,$\mbox{\boldmath $m$}_k\, } $とおいて得られる距離を, $$G_{k}\,$$と$$G_{h}\,$$の間のマハラノビス汎距離という. 平均や分散共分散行列は, 各グループに属していることがわかっている対象についての測定値より計算される. $構文解析に失敗 (構文エラー): {\displaystyle D_h(\mbox{\boldmath $x$}) (h=1,\ 2,\ \cdots,\ r)\, } $ の中で,$構文解析に失敗 (構文エラー): {\displaystyle D_k(\mbox{\boldmath $x$})\, } $ が最小であれば, この対象は,$$G_{k}\,$$に属していると判別すればよい. また, どれにも属さないという判別が許される場合は, あらかじめ上限を設定しておいて, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle D_k(\mbox{\boldmath $x$})\, } $ がそれを越えたときは, どれにも属さないと判別すればよい.

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

$${\mit\Delta}D_{12}(\mbox{\boldmath $x$}) =2(\mbox{\boldmath $m$}_2-\mbox{\boldmath $m$}_1)^{\top} C\mbox{\boldmath $x$}-(\mbox{\boldmath $m$}_2-\mbox{\boldmath $m$}_1)^{\top} C(\mbox{\boldmath $m$}_1+\mbox{\boldmath $m$}_2)$$

と変形できるので, $構文解析に失敗 (不明な関数「\mit」): {\displaystyle {\mit\Delta}D_{12}(\mbox{\boldmath $x$})\, } $は, $$x_{i}(i=1,\ 2,\ \cdots ,\ p)\,$$の線形式になる. したがって, これを($$G_{1}\,$$と$$G_{2}\,$$を判別する)線形判別関数という. $$r\,$$が3以上のときは, 線形判別関数は, $${}_{r}{\rm {C}}_{2}\,$$ 個できる. なお, 分散共分散行列が等しくないときは, $構文解析に失敗 (不明な関数「\mit」): {\displaystyle {\mit\Delta}D_{12}(\mbox{\boldmath $x$})\, } $ は, $$x_{i}(i=1,\ 2,\ \cdots ,\ p)\,$$ の2次式になる.

---

参考文献

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