En estadística , la distribución beta de matriz variable es una generalización de la distribución beta . Si es una matriz definida positiva con una distribución beta de matriz variable y son parámetros reales, escribimos (a veces ). La función de densidad de probabilidad para es:
tú
{\estilo de visualización U}
pag
×
pag
{\displaystyle p\times p}
a
,
b
>
(
pag
−
1
)
/
2
{\displaystyle a,b>(p-1)/2}
tú
∼
B
pag
(
a
,
b
)
{\displaystyle U\sim B_{p}\left(a,b\right)}
B
pag
I
(
a
,
b
)
{\displaystyle B_{p}^{I}\left(a,b\right)}
tú
{\estilo de visualización U}
{
β
pag
(
a
,
b
)
}
−
1
det
(
tú
)
a
−
(
pag
+
1
)
/
2
det
(
I
pag
−
tú
)
b
−
(
pag
+
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}.}
Aquí está la función beta multivariante:
β
pag
(
a
,
b
)
{\displaystyle \beta _{p}\left(a,b\right)}
β
pag
(
a
,
b
)
=
Γ
pag
(
a
)
Γ
pag
(
b
)
Γ
pag
(
a
+
b
)
{\displaystyle \beta _{p}(a,b\right)={\frac {\Gamma _{p}(a\right)\Gamma _{p}(b\right)}{\Gamma _{p}(a+b\right)}}}
¿Dónde está la función gamma multivariante dada por?
Γ
pag
(
a
)
{\displaystyle \Gamma _{p}\izquierda(a\derecha)}
Γ
pag
(
a
)
=
π
pag
(
pag
−
1
)
/
4
∏
i
=
1
pag
Γ
(
a
−
(
i
−
1
)
/
2
)
.
{\displaystyle \Gamma _{p}(a\right)=\pi ^{p(p-1)/4}\prod _{i=1}^{p}\Gamma \left(a-(i-1)/2\right).}
Teoremas
Distribución de la matriz inversa
Si entonces la densidad de está dada por
tú
∼
B
pag
(
a
,
b
)
{\displaystyle U\sim B_{p}(a,b)}
incógnita
=
tú
−
1
{\displaystyle X=U^{-1}}
1
β
pag
(
a
,
b
)
det
(
incógnita
)
−
(
a
+
b
)
det
(
incógnita
−
I
pag
)
b
−
(
pag
+
1
)
/
2
{\displaystyle {\frac {1}{\beta _{p}(a,b)}}\det(X)^{-(a+b)}\det \left(X-I_{p})^{b-(p+1)/2}}
siempre que y .
incógnita
>
I
pag
Estilo de visualización X>I_{p}}
a
,
b
>
(
pag
−
1
)
/
2
{\displaystyle a,b>(p-1)/2}
Si y es una matriz ortogonal constante , entonces
tú
∼
B
pag
(
a
,
b
)
{\displaystyle U\sim B_{p}(a,b)}
yo
{\estilo de visualización H}
pag
×
pag
{\displaystyle p\times p}
yo
tú
yo
yo
∼
B
(
a
,
b
)
.
{\displaystyle HUH^{T}\sim B(a,b).}
Además, si es una matriz ortogonal aleatoria que es independiente de , entonces , se distribuye independientemente de .
yo
{\estilo de visualización H}
pag
×
pag
{\displaystyle p\times p}
tú
{\estilo de visualización U}
yo
tú
yo
yo
∼
B
pag
(
a
,
b
)
{\displaystyle HUH^{T}\sim B_{p}(a,b)}
yo
{\estilo de visualización H}
Si es cualquier constante , matriz de rango , entonces tiene una distribución beta variable de matriz generalizada, específicamente .
A
{\estilo de visualización A}
q
×
pag
{\displaystyle q\veces p}
q
≤
pag
{\displaystyle q\leq p}
q
{\estilo de visualización q}
A
tú
A
yo
Estilo de visualización AUA
A
tú
A
yo
∼
GRAMO
B
q
(
a
,
b
;
A
A
yo
,
0
)
{\displaystyle AUA^{T}\sim GB_{q}\left(a,b;AA^{T},0\right)}
Resultados de la matriz particionada
Si y particionado como
tú
∼
B
pag
(
a
,
b
)
{\displaystyle U\sim B_{p}\left(a,b\right)}
tú
{\estilo de visualización U}
tú
=
[
tú
11
tú
12
tú
21
tú
22
]
{\displaystyle U={\begin{bmatrix}U_{11}&U_{12}\\U_{21}&U_{22}\end{bmatrix}}}
donde es y es , entonces definir el complemento de Schur como da los siguientes resultados:
tú
11
Estilo de visualización U_{11}}
pag
1
×
pag
1
{\displaystyle p_{1}\times p_{1}}
tú
22
Estilo de visualización U_{22}
pag
2
×
pag
2
{\displaystyle p_{2}\times p_{2}}
tú
22
⋅
1
{\displaystyle U_{22\cdot 1}}
tú
22
−
tú
21
tú
11
−
1
tú
12
{\ Displaystyle U_ {22} -U_ {21} {U_ {11}} ^ {-1} U_ {12}}
tú
11
Estilo de visualización U_{11}}
es independiente de
tú
22
⋅
1
{\displaystyle U_{22\cdot 1}}
tú
11
∼
B
pag
1
(
a
,
b
)
{\displaystyle U_{11}\sim B_{p_{1}}\left(a,b\right)}
tú
22
⋅
1
∼
B
pag
2
(
a
−
pag
1
/
2
,
b
)
{\displaystyle U_{22\cdot 1}\sim B_{p_{2}}\left(a-p_{1}/2,b\right)}
tú
21
∣
tú
11
,
tú
22
⋅
1
{\displaystyle U_{21}\mid U_{11},U_{22\cdot 1}}
tiene una distribución t variable de matriz invertida, específicamente
tú
21
∣
tú
11
,
tú
22
⋅
1
∼
I
yo
pag
2
,
pag
1
(
2
b
−
pag
+
1
,
0
,
I
pag
2
−
tú
22
⋅
1
,
tú
11
(
I
pag
1
−
tú
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).}
Resultados de Wishart
Mitra demuestra el siguiente teorema que ilustra una propiedad útil de la distribución beta variable de la matriz. Supongamos que son matrices Wishart independientes . Supongamos que es definida positiva y que . Si
S
1
,
S
2
Estilo de visualización S_{1},S_{2}
pag
×
pag
{\displaystyle p\times p}
S
1
∼
Yo
pag
(
norte
1
,
Σ
)
,
S
2
∼
Yo
pag
(
norte
2
,
Σ
)
{\displaystyle S_{1}\sim W_{p}(n_{1},\Sigma ),S_{2}\sim W_{p}(n_{2},\Sigma )}
Σ
{\estilo de visualización \Sigma}
norte
1
+
norte
2
≥
pag
{\displaystyle n_{1}+n_{2}\geq p}
tú
=
S
−
1
/
2
S
1
(
S
−
1
/
2
)
yo
,
{\displaystyle U=S^{-1/2}S_{1}\left(S^{-1/2}\right)^{T},}
donde , entonces tiene una distribución beta variable de matriz . En particular, es independiente de .
S
=
S
1
+
S
2
Estilo de visualización S=S_{1}+S_{2}}
tú
{\estilo de visualización U}
B
pag
(
norte
1
/
2
,
norte
2
/
2
)
Estilo de visualización Bp(n1/2,n2/2)
tú
{\estilo de visualización U}
Σ
{\estilo de visualización \Sigma}
Véase también
Referencias
Gupta, AK; Nagar, DK (1999). Distribuciones de variables matriciales . Chapman y Hall. ISBN 1-58488-046-5 .
Mitra, SK (1970). "Un enfoque sin densidad para la distribución beta de la matriz variable". The Indian Journal of Statistics . Serie A (1961–2002). 32 (1): 81–88. JSTOR 25049638.