Prueba multivariable de normalidad
BaringhausHenzeTest es una prueba multivariable de normalidad con la prueba estadística basada en la función de características empíricas. También se le conoce como la prueba de Baringhaus–Henze–Epps–Pulley (BHEP).
In[1]:=
![Click for copyable input](assets.es/multivariate-normality-test/In_92.png)
data = RandomVariate[NormalDistribution[], {10^3, 3}];
In[2]:=
![Click for copyable input](assets.es/multivariate-normality-test/In_93.png)
BaringhausHenzeTest[data]
Out[2]=
![](assets.es/multivariate-normality-test/O_71.png)
La prueba estadística es invariante bajo transformaciones afines de los datos.
In[3]:=
![Click for copyable input](assets.es/multivariate-normality-test/In_94.png)
data2 = AffineTransform[{RandomReal[1, {3, 3}], RandomReal[1, 3]}][
data];
{BaringhausHenzeTest[data2, "TestStatistic"],
BaringhausHenzeTest[data, "TestStatistic"]}
Out[3]=
![](assets.es/multivariate-normality-test/O_72.png)
La prueba estadística también es consistente bajo cada distribución alternativa; es decir, esta crece ilimitadamente con el tamaño de muestra a menos que provenga de una distribución gaussiana.
In[4]:=
![Click for copyable input](assets.es/multivariate-normality-test/In_95.png)
covm = {{2, 1, 0}, {1, 3, -1}, {0, -1, 2}};
ng\[ScriptCapitalD] = MultivariateTDistribution[covm, 12];
g\[ScriptCapitalD] = MultinormalDistribution[{0, 0, 0}, covm];
Extraiga muestras de una distribución t multivariable y de una distribución normal multivariable.
In[5]:=
![Click for copyable input](assets.es/multivariate-normality-test/In_96.png)
nongaussianData = RandomVariate[ng\[ScriptCapitalD], 2^14];
gaussianData = RandomVariate[g\[ScriptCapitalD], 2^14];
Calcule y visualice la prueba estadística de BHEP como una función de tamaño de muestra.
In[6]:=
![Click for copyable input](assets.es/multivariate-normality-test/In_97.png)
sizes = 2^Range[3, 14];
{nonGaussianTestStatistic, GaussianTestStatistic} =
ParallelTable[{size,
BaringhausHenzeTest[Take[d, size],
"TestStatistic"]}, {d, {nongaussianData, gaussianData}}, {size,
sizes}];
muestre la entrada completa de Wolfram Language
Out[7]=
![](assets.es/multivariate-normality-test/O_73.png)