【徹底解説】多変量正規分布の独立性

本記事は「これなら分かる!はじめての数理統計学」シリーズに含まれます。

不適切な内容があれば,記事下のコメント欄またはお問い合わせフォームよりご連絡下さい。

多変量正規分布の独立性

XN(μ,Σ)とする。X1RD1X2RD2が独立となる必要十分条件は,

(1)Σ12=OD1,D2

である。ただし,以下のようにX,μ,t,ΣD1次元の部分とD2=DD1次元の部分に区分けした。

(2)X=[X1X2]μ=[μ1μ2]t=[t1t2]Σ=[Σ11Σ12Σ21Σ22]

多変量正規分布の独立性に関する定理です。D1=1を逐次適用していくことで,全ての要素が独立である必要十分条件は分散共分散行列が対角行列となることが分かります。

証明

多変量正規分布に従う確率変数が独立となる必要十分条件を考えますので,必要性と十分性に分けて考える必要があります。一方で,X1X2が独立となる必要十分条件は「Xのモーメント母関数がX1のモーメント母関数とX2のモーメント母関数の積で表されること」ですので,この必要十分条件を同値変形していきます。

多変量正規分布のモーメント母関数の定義より,Xのモーメント母関数は以下のように表されます。

(3)MX(t)=exp(μTt+12tTΣt)

同様に,X1X2のモーメント母関数は以下のようになります。

(4)MX1(t1)=exp(μ1Tt1+12t1TΣ11t1)(5)MX2(t2)=exp(μ2Tt2+12t2TΣ22t2)

冒頭にある通り,以下の恒等式を考えます。

(6)MX(t)=MX1(t1)MX2(t2)

すなわち,X1X2が独立となる条件は以下です。

(7)exp(μTt+12tTΣt)=exp(μ1Tt1+12t1TΣ11t1)exp(μ2Tt2+12t2TΣ22t2)

指数部分のみを考えます。

(8)μTt+12tTΣt=μ1Tt1+μ2Tt2+12t1TΣ11t1+12t2TΣ22t2

左辺を行列の分解形式で計算します。

[μ1μ2]T[t1t2]+12[t1t2]T[Σ11Σ12Σ21Σ22][t1t2](9)=μ1t1+μ2t2+12t1TΣ11t1+12t2TΣ22t2+12t1TΣ12t1+12t2TΣ21t2

式(8)に式(9)を代入すると,以下が得られます。

(10)12t1TΣ12t1+12t2TΣ21t2=0

t1t2はモーメント母関数の変数ですので,ゼロベクトル以外の値も取り得ます。したがって,以下が得られます。

(11)Σ12=OD1,D2(12)Σ21=OD2,D1

以上より,X1X2が独立となる必要十分条件は「Σ12Σ21がゼロ行列であること」を証明できました。

参考文献

本稿の執筆にあたり参考にした文献は,以下でリストアップしております。

シェアはこちらからお願いします!

コメント

コメントする

※ Please enter your comments in Japanese to distinguish from spam.