Loi de Marchenko–Pastur
La loi de Marchenko-Pastur est la loi limite des valeurs propres des matrices de Wishart lorsque la dimension de la matrice et les degrés de liberté
tendent tous deux vers l'infini avec le rapport
. Pour
, la distribution n'a pas de masse ponctuelle et la fonction de densité de probabilité est bien définie.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_70.png)
PDF[MarchenkoPasturDistribution[1/2], x]
![](assets.fr/marchenko-pastur-distribution/O_31.png)
![](assets.fr/marchenko-pastur-distribution/O_32.png)
Échantillonnez à partir d'une loi de Wishart avec une matrice d'échelle identique et calculez les valeurs propres mises à l'échelle.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_72.png)
n = 10^4;
m = 10^3;
eigs = RandomVariate[
MatrixPropertyDistribution[Eigenvalues[x]/n,
x \[Distributed]
WishartMatrixDistribution[n, IdentityMatrix[m]]]];
Comparez le résultat échantillonné avec la fonction de densité de Marchenko–Pastur.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_73.png)
Show[Histogram[eigs, {0.05}, "PDF", ImageSize -> Medium,
PlotTheme -> "Detailed"],
Plot[PDF[MarchenkoPasturDistribution[m/n], x], {x, 0, 1.8},
PlotTheme -> "Detailed", PlotLegends -> None, Exclusions -> None]]
![](assets.fr/marchenko-pastur-distribution/O_33.png)
Pour , la matrice de Wishart est singulière. Avec la probabilité
, la distribution a une masse en
.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_74.png)
m = 500; n = 2 m;
CDF[MarchenkoPasturDistribution[n/m], 0]
![](assets.fr/marchenko-pastur-distribution/O_34.png)
Générez une matrice de Wishart singulière avec une covariance identique et calculez les valeurs propres échelonnées.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_75.png)
matrix = Transpose[#].# &[RandomVariate[NormalDistribution[], {m, n}]];
eigvs = Chop[Eigenvalues[matrix]/m];
Il y a un écart dans la densité des valeurs propres près de 0, et la case à 0 a une grande densité.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_76.png)
Histogram[eigvs, {0.05}, PDF, PlotRange -> 1, ChartStyle -> Orange,
ImageSize -> Medium]
![](assets.fr/marchenko-pastur-distribution/O_35.png)
Ajustez MarchenkoPasturDistribution aux valeurs propres.
![Click for copyable input](assets.fr/marchenko-pastur-distribution/In_77.png)
edist = EstimatedDistribution[eigvs,
MarchenkoPasturDistribution[\[Lambda], 1]]
![](assets.fr/marchenko-pastur-distribution/O_36.png)
La fonction de distribution cumulative de la loi ajustée montre une discontinuité de saut à l'origine.
![](assets.fr/marchenko-pastur-distribution/O_37.png)