nishiru3の日記

備忘録です。ネットのゴミ。

偏相関係数の導出

概要

x,y,zの3変数を考える。目的変数をy、説明変数をxとし、背後にzがあると想定する。ここではzの影響を取り除いたxy相関係数、すなわち偏相関係数r_{xy,z}を導出することを目的とする。

導出

まず、x,zでの回帰を考える。x=az+bへ回帰するものとすると、係数a,bは次のとおりである。

f:id:nishiru3:20190429135539p:plain

ここに、\bar{x},\bar{z}は平均値である。この回帰式から求められるxの推定値を\hat{x}とすると残差x^{\prime}は、

f:id:nishiru3:20190429135852p:plain

である。同様にy,zでの回帰を考え、y=az+bへ回帰するものとすると、残差y^{\prime}は、

f:id:nishiru3:20190429140107p:plain

である。zの影響を取り除いたxy相関係数r_{xy,z}は上記の残差x^\primey^\prime相関係数を取ればよい。
定義式より変形して整理すると、

f:id:nishiru3:20190429140318p:plain


残差X^\primeY^\primeの平均は0であるので、

f:id:nishiru3:20190429140505p:plain

となり、E[X^\prime Y^\prime]E[X^{\prime 2}]およびE[Y^{\prime 2}]を求めればよい。

E[X^\prime Y^\prime]の計算

f:id:nishiru3:20190429140927p:plain

E[X^{\prime 2}]およびE[Y^{\prime 2}]の計算

f:id:nishiru3:20190429141159p:plain

相関係数の導出

f:id:nishiru3:20190429141313p:plain

f:id:nishiru3:20190429141350p:plain