次の表 6 は、同品種のコーヒー豆の品質に関わる3項目(苦味、酸味、香り)と、その生産地の関係を調べたものである。

表 6 同品種のコーヒー豆の品質と生産地の関係
 
No.苦味($x_1$)酸味($x_2$)香り($x_3$)生産地($y$)
13.73.03.4A
23.43.12.6A
33.93.53.9A
43.42.63.8A
53.53.73.5A
64.43.63.2A
74.93.13.5A
84.63.42.9A
93.03.93.9A
104.43.83.9A
112.92.83.9B
123.31.73.8B
133.22.24.4B
143.91.43.3B
152.61.04.7B
163.72.74.6B
173.72.94.6B
183.51.24.2B
193.71.43.5B
203.21.03.7B
 

層別散布図により、3 つの変数中の任意の 2 つの変数の組み合わせによる生産地の判別がどの程度可能であるかを検討し、生産地を判別するための式を求めよ。

考え方と適用手法

ある観測対象が所属するグループを予測するための手法のひとつに判別分析がある。グループの予測とは、たとえば、良品グループと不良品グループ、購入者と非購入者といったように、2 つ以上に分類されるグループの中で、どのグループに属するかを予測することである。

また、変数が複数($\geqq 3$)ある場合には、層別散布図を用いて、判別の可能性を視覚的に検討することができる。

層別散布図

変数が 2 つ(以下)の場合には、散布図により観測データ間の関係を視覚的に把握することができる。しかし、変数が 3 つ以上になると、3 次元、4 次元のグラフとなり、もはや視覚的に把握することは困難となる。そこで、任意の 2 変数をとり出し、複数の散布図から把握することを考える。このような散布図を層別散布図という。

plots1.png
図 8a 苦味 $x_1$ と酸味 $x_2$ の層別散布図
plots2.png
図 8b 酸味 $x_2$ と香り $x_3$ の層別散布図
plots3.png
図 8c 香り $x_3$ と苦味 $x_1$ の層別散布図
3d_plots.png
図 8d 苦味と酸味と香りの 3D 散布図
 

判別分析

$k$ 個($k\geqq 2$)のグループが存在していて、どのグループに属するかわかっている観測対象について、何らかのデータが収集されているとする。このときに、どのグループに属するか不明な対象が得られた場合、すでに収集されているデータにもとづいて、その対象が属するグループを予測するのが判別分析である。ここでは、グループが 2 つの場合についてのみ考える。

判別分析にはいくつかの方法があるが、ここでは回帰分析の応用として、重回帰分析を用いた方法をとりあげる。すなわち、変数 $x_1 \sim x_3$ と $y$ の関係を超平面の式($y=0$ とおいたときの残りの変数が 2 つなら 2 次元(平面)を 2 つに分かつ直線の式、3 つなら 3 次元(空間)を 2 つに分かつ平面の式であり、ここでは 図8d の白丸と黒丸を分かつ平面の式となる)

$$ y=b_0+b_1x_1+b_2x_2+b_3x_3 $$

で表すことを考える。判別分析ではこの式を判別式とよぶ。

しかし、ここで問題になるのが従属変数 $y$ である。$y$ は生産地 A、B といった質的変量であるため、このままでは利用できない。そこで、生産地 A ならば $y=1$、B ならば $y=-1$ というように量的変量に変換して利用する。

Excel では [ツール]-[分析ツール]-[回帰分析] を利用すれば、上式の $b_0\sim b_3$ を簡単に求めることができる。また、判別式の係数や予測値は Excel 関数 : LINEST および TREND で求めることができる。

[授業計画に戻る]

添付ファイル: file3D_plots.3mf 55件 [詳細] file3d_plots.png 45件 [詳細] fileplots3.png 44件 [詳細] fileplots2.png 44件 [詳細] fileplots1.png 45件 [詳細]

トップ   編集 凍結解除 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2023-11-30 (木) 22:07:51