生物統計学 学習ノート概念・導出から実践的な解釈まで

CHAPTER 07 · TOPIC 03

Fisher Z変換と相関係数の比較

相関係数を正規化し、母相関の信頼区間と独立二群の相関差を推測します。

このページの内容
  1. Fisher Z変換
  2. ρが指定値に等しいかを検定する
  3. 母相関の信頼区間
  4. 独立二群の相関を比較する
  5. この二群式を使えない場合

標本相関係数rの標本分布は、特に|ρ|が大きいと非対称で分散もρに依存します。Fisher Z変換はrを実数全体へ写し、分布をより正規かつ分散一定に近づけ、信頼区間や相関係数の比較を可能にします。

Fisher Z変換#

標本相関zr=12ln(1+r1r)=arctanh(r)z_r=\frac12\ln\left(\frac{1+r}{1-r}\right)=\operatorname{arctanh}(r)
母相関zρ=12ln(1+ρ1ρ)z_\rho=\frac12\ln\left(\frac{1+\rho}{1-\rho}\right)
逆変換r=tanh(zr)=e2zr1e2zr+1r=\tanh(z_r)=\frac{e^{2z_r}-1}{e^{2z_r}+1}
z(r)N(z(ρ),1n3),SE[z(r)]1n3z(r)\approx N\left(z(\rho),\frac1{n-3}\right),\qquad SE[z(r)]\approx\frac1{\sqrt{n-3}}

ρが指定値に等しいかを検定する#

仮説H0:ρ=ρ0,H1:ρρ0H_0:\rho=\rho_0,\qquad H_1:\rho\ne\rho_0
Z統計量Z=z(r)z(ρ0)1/n3Z=\frac{z(r)-z(\rho_0)}{1/\sqrt{n-3}}

ρ₀=0の検定には前頁のt検定が正確な標準法です。Fisher Zは非零の指定値、信頼区間、独立相関の比較で特に有用です。

母相関の信頼区間#

下限Lz=z(r)z1α/21n3L_z=z(r)-z_{1-\alpha/2}\frac1{\sqrt{n-3}}
上限Uz=z(r)+z1α/21n3U_z=z(r)+z_{1-\alpha/2}\frac1{\sqrt{n-3}}
r尺度へ戻すCIρ=[tanh(Lz),tanh(Uz)]CI_\rho=[\tanh(L_z),\tanh(U_z)]

独立二群の相関を比較する#

仮説H0:ρ1=ρ2,H1:ρ1ρ2H_0:\rho_1=\rho_2,\qquad H_1:\rho_1\ne\rho_2
各群を変換z1=z(r1),z2=z(r2)z_1=z(r_1),\qquad z_2=z(r_2)
差の標準誤差SE(z1z2)=1n13+1n23SE(z_1-z_2)=\sqrt{\frac1{n_1-3}+\frac1{n_2-3}}
比較ZZ=z(r1)z(r2)1/(n13)+1/(n23)Z=\frac{z(r_1)-z(r_2)}{\sqrt{1/(n_1-3)+1/(n_2-3)}}
(z1z2)±z1α/21n13+1n23(z_1-z_2)\pm z_{1-\alpha/2}\sqrt{\frac1{n_1-3}+\frac1{n_2-3}}

この二群式を使えない場合#

  • 同じ対象から得た二つの相関など、相関係数が互いに依存する
  • 二相関が一つの変数を共有する
  • クラスター、反復測定、重み付き標本がある
  • Pearson相関の仮定が著しく破れている
  • Spearman相関へ標準誤差式をそのまま流用する