【過去問解答】2021年統計検定1級<数理統計問5>

統計検定1級の過去問解答解説を行います。目次は以下をご覧ください。

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

問題

統計検定1級の過去問からの出題になります。統計検定の問題の著作権は日本統計学会に帰属していますので,本稿にて記載することはできません。「演習問題を俯瞰する」で詳しく紹介している公式の過去問題集をご購入いただきますようお願い致します。

解答

多変量正規分布の独立性に関する出題です。もはや線型代数の問題と言っても過言ではないです。以下のページが参考になるかもしれません。

(1)

(1)V[y]=LLT

xの各要素は独立に標準正規分布に従うことから,Onn×n次元ゼロ行列,Inn次元単位行列とすると,x多変量正規分布であるN(On,In)に従います。したがって,分散共分散行列V[x]Inとなることに注意すると,分散共分散行列の性質より,

(2)V[y]=V[Lx]=LV[x]LT=LInLT=LLT

が得られます。

(2)

LMT=Oが必要十分条件である

まず,多変量正規分布と線形変換の性質から,LxMxも多変量正規分布に従います。次に,多変量正規分布における無相関と独立の性質から,LxMxが無相関であることと独立であることは同値となります。無相関の定義より,行列表記の共分散を用いると,LxMxが独立となる必要十分条件は

(3)Cov[Lx,Mx]=On

となります。ここで,行列表記の共分散の性質Cov[x,x]xの分散共分散行列であることを用いると,

(4)Cov[Lx,Mx]=LCov[x,x]MT=LInMT=LMT

が得られます。式(3)と式(4)より,LMT=Oが必要十分条件であることが導かれます。

行列表記の共分散の性質は定義を同値変形していますので,必要十分条件として利用することができます。

(3)

(5)A=LMT(MMT)1

前問(2)と同様に,LxAMxMxが独立となる必要十分条件は,行列表記の共分散の性質より

(6)Cov[LxAMx,Mx]=LCov[x,x]MTAMCov[x,x]MT=LMTAMMT=On

となります。すなわち,

(7)AMMT=LMT

となります。ここで,フルランクな行列は正則であることから,Mは正則になります。また,ある行列が正則であることと行列式が0でないことは同値であること,積の行列式は行列式の積となること,および転置行列の行列式ともとの行列式は等しくなることから,

(8)det(MMT)=det(M)det(MT)=det(M)20

が得られます。したがって,MMTは正則となりますので,式(7)の右からMMTの逆行列を掛けることにより,

(9)A=LMT(MMT)1

が得られます。

(4)

任意の行列は特異値分解が可能であるため,Mの特異値分解により,任意のL,Mに対して条件を満たすAが存在する。

任意の行列は特異値分解が可能であることから,Mの左特異ベクトルを列ベクトルとする行列URm×r,右特異ベクトルを列ベクトルとする行列VRn×r,特異値を対角要素にもつ対角行列ΣRr×rを用いてM特異値分解すると,

(10)M=UΣVT

が得られます。これを式(7)に代入すると,

(11)AUΣVTVΣUT=LVΣUT

が得られます。いま,特異値分解の存在の証明によると,UVは直交行列となります。ゆえに,直交行列の転置行列と逆行列が等しくなることを利用すると,

(12)A=(LVΣUT)(UΣ1VTVΣ1UT)=LVΣ1UT

が得られます。ただし,Σは特異値を対角要素にもつ対角行列ですので,r個の特異値のうち,s個の非ゼロの特異値を大きい順にσ1,,σsとおくと,

(13)Σ1=[1/σ11/σs00]

となります。したがって,任意のL,Mに対して条件を満たすAが式(12)の形で存在することが示されました。

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

コメント

コメントする

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

目次
  1. 問題
  2. 解答
    1. (1)
    2. (2)
    3. (3)
    4. (4)