Marchenko–Pastur 분포
Marchenko–Pastur 분포는 행렬 차원 , 그리고 자유도
, 둘다
의 비율의 무한대 경향을 가진 위샤트 행렬의 고유값의 극한 분포입니다.
에 대해 분포는 질점을 가지지않으며 확률 밀도 함수는 명확하게 정의됩니다.
In[1]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_71.png)
PDF[MarchenkoPasturDistribution[1/2], x]
Out[1]=
![](assets.ko/marchenko-pastur-distribution/O_31.png)
전체 Wolfram 언어 입력 표시하기
Out[2]=
![](assets.ko/marchenko-pastur-distribution/O_32.png)
항등 척도 행렬이있는 위샤트 분포에서 샘플을 채취하고 스케일된 고유값을 계산합니다.
In[3]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_73.png)
n = 10^4;
m = 10^3;
eigs = RandomVariate[
MatrixPropertyDistribution[Eigenvalues[x]/n,
x \[Distributed]
WishartMatrixDistribution[n, IdentityMatrix[m]]]];
샘플된 결과를 Marchenko–Pastur 밀도 함수와 비교합니다.
In[4]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_74.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]]
Out[4]=
![](assets.ko/marchenko-pastur-distribution/O_33.png)
에대해 위샤트 행렬은 특이 행렬입니다. 확률
의 경우, 분포는
에서 질점을 가집니다.
In[5]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_75.png)
m = 500; n = 2 m;
CDF[MarchenkoPasturDistribution[n/m], 0]
Out[5]=
![](assets.ko/marchenko-pastur-distribution/O_34.png)
항등 공분산을 가지는 특이 위샤트 행렬을 생성하고 스케일 된 고유값을 계산합니다.
In[6]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_76.png)
matrix = Transpose[#].# &[RandomVariate[NormalDistribution[], {m, n}]];
eigvs = Chop[Eigenvalues[matrix]/m];
0 부근에서 고유값의 밀도에 차이가 있고, 0의 빈에서 밀도가 높아지고 있습니다.
In[7]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_77.png)
Histogram[eigvs, {0.05}, PDF, PlotRange -> 1, ChartStyle -> Orange,
ImageSize -> Medium]
Out[7]=
![](assets.ko/marchenko-pastur-distribution/O_35.png)
MarchenkoPasturDistribution을 고유값에 맞춥니다.
In[8]:=
![Click for copyable input](assets.ko/marchenko-pastur-distribution/In_78.png)
edist = EstimatedDistribution[eigvs,
MarchenkoPasturDistribution[\[Lambda], 1]]
Out[8]=
![](assets.ko/marchenko-pastur-distribution/O_36.png)
적합된 분포의 누적 분포 함수는 원점의 불연속 점프를 보입니다.
전체 Wolfram 언어 입력 표시하기
Out[9]=
![](assets.ko/marchenko-pastur-distribution/O_37.png)