例題:
「表 1 に示すデータにおいて,変数 $Z$ の影響を取り除いたときの,変数 $X$ と変数 $Y$ の間の偏相関係数を求めなさい。」
$X$ | $Y$ | $Z$ |
---|---|---|
61.0 | 52.0 | 58.0 |
57.3 | 44.2 | 45.7 |
55.5 | 63.3 | 62.6 |
40.8 | 29.5 | 28.0 |
27.9 | 46.8 | 52.5 |
57.0 | 50.6 | 67.1 |
57.5 | 41.6 | 49.1 |
47.7 | 67.8 | 51.6 |
57.1 | 43.1 | 54.4 |
56.6 | 47.8 | 44.9 |
43.7 | 46.4 | 39.0 |
61.7 | 60.0 | 62.1 |
32.0 | 45.7 | 45.0 |
44.4 | 49.3 | 38.3 |
50.1 | 65.3 | 50.9 |
説明
2 変数の相関関係を調べるとき,それぞれの変数が第 3 の変数と相関を持っていることにより,観察された相関が真の相関関係を表さないことがある。
そこで,対象としている 2 変数以外の変数の影響を取り除いた相関係数を求めることが必要になることがある。
$p$ 個の変数があるとき,$q$ 個の変数を固定した(影響を取り除いた)$p-q$ 個の変数間の偏相関係数は「 $q$ 次のオーダーの偏相関係数 」と呼ばれる。
2 変数 $i$,$j$ 間の相関係数を $r_{ij}$ とし,固定される変数を $k_{1}, k_{2}, \dots , k_{q}$ の $q$ 個としたとき,2 変数 $i$,$j$ 間の $q$ 次のオーダーの偏相関係数は $(1)$ 式,$(2)$ 式の漸化式で計算できる。
\[
\begin{align*}
r(ii\cdot k_1) = \frac{r(ij) - r(ik_1)\ r(jk_1)} {\sqrt{1-r(ik_1)^2}\ \sqrt{1-r(jk_1)^2}} \tag{1}\\
r(ij){}_q = \frac{r(ij){}_{q-1} - r(ik_q){}_{q-1}\ r(jk_q){}_{q-1}}{\sqrt{1-r(ik_q){}_{q-1}^2}\ \sqrt{1-r(jk_q){}_{q-1}^2}} \tag{2}
\end{align*}
\]
なお,$(2)$ 式では,$r(ij){}_{q-1} =r(ij\cdot k_1\,k_2\dots k_{q-1})$ などと略記した。
2 変数 $i$,$j$ 以外の全ての変数を固定したときの偏相関係数は,相関係数行列 $\mathbf{r}$ の逆行列の要素を $r^{ij}$ としたとき,$(3)$ 式で求めることができる。 \[ r(ij\cdot k\ne i, j) = \frac{-r^{ij}}{\sqrt{r^{ii}\ r^{jj}}} \tag{3} \] 注
3 変数の場合に,1 変数を固定した 2 変数間の偏相関係数を求めるとき,$(1)$ 式は複雑であるので $(4)$ 式の方がわかりやすいかもしれない。 \[ r_{xy\cdot z} = \frac{r_{xy}-r_{xz}\ r_{yz}} {\sqrt{1-r_{xz}^2}\ \sqrt{1-r_{yz}^2}} \tag{4} \]
計算手順:
$X$ | $Y$ | $Z$ | |
---|---|---|---|
$X$ | 1.00000 | ||
$Y$ | 0.25058 | 1.00000 | |
$Z$ | 0.49484 | 0.59803 | 1.00000 |
$\displaystyle \frac{0.25058 - 0.49484 \times 0.59803} {\sqrt{ 1 - 0.49484^{2}} \ \sqrt{ 1 - 0.59803^{2}}} ≒ -0.06511$
(注:コンピュータを用いて計算したときは,小数点以下 5 桁までの答えは $-0.06512$ である)。
演習問題:
例題において,「変数 $Z$の影響を取り除いたときの」という表現がある。変数 $X$を従属変数,変数 $Z$を独立変数として直線回帰を行うとする。そのときの変数 $X$の残差(実測値から予測値を引いたもの)は,「変数 $Z$の影響を取り除いた(変数 $Z$によらない)値になる。同じように,変数 $Y$を従属変数,変数 $Z$を独立変数 として,変数 $Y$の残差を求める。このようにして求めた 2 つの残差の間の相関係数を求めると,それは,変数 $Z$の影響を取り除いたときの変数 $X$と変数 $Y$の偏相関係数になる。表 1.のデータを使ってこれを確かめよ。
応用問題: