行列変量ベータ分布

行列変量ベータ分布
表記Bp1つのb{\displaystyle {\rm {B}}_{p}(a,b)}
パラメータ1つのb{\displaystyle a,b}
サポートp×p{\displaystyle p\times p}正定値と負定値の両方を持つ行列あなた{\displaystyle U}pあなた{\displaystyle I_{p}-U}
PDF{βp1つのb}1詳細あなた1つのp+1/2詳細pあなたbp+1/2{\displaystyle \left\{\beta _{p}\left(a,b\right)\right\}^{-1}\det \left(U\right)^{a-(p+1)/2}\det \left(I_{p}-U\right)^{b-(p+1)/2}.}
CDF1F11つの;1つの+b;Z{\displaystyle {}_{1}F_{1}\left(a;a+b;iZ\right)}

統計学において、行列変量ベータ分布はベータ分布の一般化である。MANOVAアンサンブルヤコビアンサンブルとも呼ばれる。

が行列変量ベータ分布に従う 正定値行列であり、が実パラメータである場合、 ( と表記されることもある)の確率密度関数は以下の通りである。あなた{\displaystyle U}p×p{\displaystyle p\times p}1つのb>p1/2{\displaystyle a,b>(p-1)/2}あなたBp1つのb{\displaystyle U\sim B_{p}\left(a,b\right)}Bp1つのb{\displaystyle B_{p}^{I}\left(a,b\right)}あなた{\displaystyle U}{βp1つのb}1詳細あなた1つのp+1/2詳細pあなたbp+1/2{\displaystyle \left\{\beta _{p}\left(a,b\right)\right\}^{-1}\det \left(U\right)^{a-(p+1)/2}\det \left(I_{p}-U\right)^{b-(p+1)/2}.}

多変数ベータ関数は次のとおりです。 βp1つのb{\displaystyle \beta _{p}\left(a,b\right)}

βp1つのbΓp1つのΓpbΓp1つの+b{\displaystyle \beta_{p}\left(a,b\right)={\frac {\Gamma_{p}\left(a\right)\Gamma_{p}\left(b\right)}{\Gamma_{p}\left(a+b\right)}}}

ここで、多変数ガンマ関数は次のように与えられる。 Γp1つの{\displaystyle \Gamma _{p}\left(a\right)}

Γp1つのπpp1/41pΓ1つの1/2{\displaystyle \Gamma_{p}\left(a\right)=\pi^{p(p-1)/4}\prod_{i=1}^{p}\Gamma\left(a-(i-1)/2\right).}

定理

逆行列の分布

の密度は次のように与えられる。 あなたBp1つのb{\displaystyle U\sim B_{p}(a,b)}Xあなた1{\displaystyle X=U^{-1}}

1βp1つのb詳細X1つの+b詳細Xpbp+1/2{\displaystyle {\frac {1}{\beta _{p}\left(a,b\right)}}\det(X)^{-(a+b)}\det \left(X-I_{p}\right)^{b-(p+1)/2}}

ただし、 およびとする。 X>p{\displaystyle X>I_{p}}1つのb>p1/2{\displaystyle a,b>(p-1)/2}

直交変換

が定数直交行列である場合、あなたBp1つのb{\displaystyle U\sim B_{p}(a,b)}H{\displaystyle H}p×p{\displaystyle p\times p}HあなたHTB1つのb{\displaystyle HUH^{T}\sim B(a,b).}

また、が に依存しないランダム直交行列である場合、 はに依存しない分布になります。 H{\displaystyle H}p×p{\displaystyle p\times p}あなた{\displaystyle U}HあなたHTBp1つのb{\displaystyle HUH^{T}\sim B_{p}(a,b)}H{\displaystyle H}

が任意の定数で階数の行列である場合、 は一般化行列変量ベータ分布(具体的には )に従います。 {\displaystyle A}q×p{\displaystyle q\times p}qp{\displaystyle q\leq p}q{\displaystyle q}あなたT{\displaystyle AUA^{T}}あなたTGBq1つのb;T0{\displaystyle AUA^{T}\sim GB_{q}\left(a,b;AA^{T},0\right)}

分割された行列の結果

とすると、次のように 分割する。あなたBp1つのb{\displaystyle U\sim B_{p}\left(a,b\right)}あなた{\displaystyle U}

あなた[あなた11あなた12あなた21あなた22]{\displaystyle U={\begin{bmatrix}U_{11}&U_{12}\\U_{21}&U_{22}\end{bmatrix}}}

ここで、はであり、シュアー補数を と定義すると、次の結果が得られます。 あなた11{\displaystyle U_{11}}p1×p1{\displaystyle p_{1}\times p_{1}}あなた22{\displaystyle U_{22}}p2×p2{\displaystyle p_{2}\times p_{2}}あなた221{\displaystyle U_{22\cdot 1}}あなた22あなた21あなた111あなた12{\displaystyle U_{22}-U_{21}{U_{11}}^{-1}U_{12}}

  • あなた11{\displaystyle U_{11}}独立しているあなた221{\displaystyle U_{22\cdot 1}}
  • あなた11Bp11つのb{\displaystyle U_{11}\sim B_{p_{1}}\left(a,b\right)}
  • あなた221Bp21つのp1/2b{\displaystyle U_{22\cdot 1}\sim B_{p_{2}}\left(a-p_{1}/2,b\right)}
  • あなた21あなた11あなた221{\displaystyle U_{21}\mid U_{11},U_{22\cdot 1}}逆行列変量t分布を持ち、具体的にはあなた21あなた11あなた221Tp2p12bp+10p2あなた221あなた11p1あなた11{\displaystyle U_{21}\mid U_{11},U_{22\cdot 1}\sim IT_{p_{2},p_{1}}\left(2b-p+1,0,I_{p_{2}}-U_{22\cdot 1},U_{11}(I_{p_{1}}-U_{11})\right).}

ウィシャートの結果

ミトラは、行列変量ベータ分布の有用な性質を示す次の定理を証明している。独立なウィシャート行列とする。 が正定値行列であり、 であるとする。もし S1S2{\displaystyle S_{1},S_{2}}p×p{\displaystyle p\times p}S1Wpn1ΣS2Wpn2Σ{\displaystyle S_{1}\sim W_{p}(n_{1},\Sigma ),S_{2}\sim W_{p}(n_{2},\Sigma )}Σ{\displaystyle \Sigma }n1+n2p{\displaystyle n_{1}+n_{2}\geq p}

あなたS1/2S1S1/2T{\displaystyle U=S^{-1/2}S_{1}\left(S^{-1/2}\right)^{T},}

ここでは行列変量ベータ分布 に従う。特に、は とは独立である。 SS1+S2{\displaystyle S=S_{1}+S_{2}}あなた{\displaystyle U}Bpn1/2n2/2{\displaystyle B_{p}(n_{1}/2,n_{2}/2)}あなた{\displaystyle U}Σ{\displaystyle \Sigma }

スペクトル密度

スペクトル密度はヤコビ多項式で表される。[ 1 ]

極値分布

最大固有値の分布は、トレイシー・ウィダム分布の変換によってよく近似されます。[ 2 ]

参照

参考文献

  • ポッターズ、マーク;ブショー、ジャン=フィリップ( 2020年11月30日)「7. ヤコビアンサンブル」ランダム行列理論入門:物理学者、エンジニア、データサイエンティスト向け。ケンブリッジ大学出版局。doi 10.1017 /9781108768900。ISBN 978-1-108-76890-0
  • フォレスター、ピーター (2010). 「3. ラゲール・アンサンブルとヤコビ・アンサンブル」.対数気体とランダム行列. ロンドン数学会モノグラフ. プリンストン: プリンストン大学出版局. ISBN 978-0-691-12829-0
  • アンダーソン, GW; ギオネ, A.; ゼイトゥニ, O. (2010). 「4. いくつかの一般論」.ランダム行列入門. ケンブリッジ: ケンブリッジ大学出版局. ISBN 978-0-521-19452-5
  • Mehta, ML (2004). 「19. 行列アンサンブルと古典直交多項式」.ランダム行列. アムステルダム: Elsevier/Academic Press. ISBN 0-12-088409-7
  • Gupta, AK; Nagar, DK (1999).行列変量分布. Chapman and Hall. ISBN 1-58488-046-5
  • カトリ, CG (1992). 「行列ベータ分布とその線形モデル、検定、歪度、尖度への応用」. ヴェヌゴパル, N. (編). 『確率論への貢献』 . John Wiley & Sons. pp.  26– 34. ISBN 0-470-22050-3
  • ミトラ, SK (1970). 「行列変量ベータ分布への密度フリーアプローチ」.インド統計ジャーナル. シリーズA (1961–2002). 32 (1): 81– 88. JSTOR  25049638 .