Distribuciones de Wishart e invertidas de Wishart
La distribución de Wishart es la distribución de matriz de covarianza de muestras de vectores aleatorios multinormales independientes. Es una generalización de la distribución en múltiples dimensiones. La distribución aparece naturalmente en estadística multivariante, tal como regresión, covarianza, etc.
Genere una matriz definida positiva para usar como parámetro de la distribución de Wishart.
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_24.png)
\[CapitalSigma] = DiagonalMatrix[RandomReal[10, 5]];
Las matrices de la distribución de Wishart son simétricas y definidas positivas. »
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_25.png)
dist = WishartMatrixDistribution[30, \[CapitalSigma]];
mat = RandomVariate[dist];
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_26.png)
SymmetricMatrixQ[mat] && PositiveDefiniteMatrixQ[mat]
![](assets.es/wishart-and-inverse-wishart-distributions/O_12.png)
La distribución inversa de Wishart es la distribución de matrices inversas de la distribución de Wishart. »
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_27.png)
invdist =
InverseWishartMatrixDistribution[30, Inverse[\[CapitalSigma]]];
invmat = RandomVariate[invdist];
Las matrices de la distribución inversa de Wishart son simétricas y definidas positivas.
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_28.png)
SymmetricMatrixQ[invmat] && PositiveDefiniteMatrixQ[invmat]
![](assets.es/wishart-and-inverse-wishart-distributions/O_13.png)
Compare la distribución de valores propios para matrices de distribuciones de Wishart e inversas de Wishart.
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_29.png)
eigs = Flatten[
RandomVariate[
MatrixPropertyDistribution[Eigenvalues[x], x \[Distributed] dist],
10^4]];
inveigs =
Flatten[RandomVariate[
MatrixPropertyDistribution[Eigenvalues[x]^-1,
x \[Distributed] invdist], 10^4]];
![](assets.es/wishart-and-inverse-wishart-distributions/O_14.png)
Para un vector distinto de cero matriz de Wishart
con matriz de escala
, la estadística
es distribuida en
.
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_31.png)
y = #/Sqrt[#.\[CapitalSigma].#] &[RandomReal[1, 5]];
data = RandomVariate[
MatrixPropertyDistribution[y.w.y,
w \[Distributed] WishartMatrixDistribution[30, \[CapitalSigma]]],
10^4];
![Click for copyable input](assets.es/wishart-and-inverse-wishart-distributions/In_32.png)
Show[Histogram[data, Automatic, PDF, PlotTheme -> "Detailed"],
Plot[PDF[ChiSquareDistribution[30], x], {x, 0, 80}],
ImageSize -> Medium]
![](assets.es/wishart-and-inverse-wishart-distributions/O_15.png)