Я пытаюсь оценить вероятность того, что наблюдение будет классифицировано как одна из двух возможных многомерных нормалей. Используя пакет R MixSim, вот мои настройки и вывод.
> Mu
[,1] [,2] [,3] [,4]
[1,] 1.071975 17.81168 122934.80 186201.2
[2,] 1.127719 23.85656 97991.84 211144.2
> S
, , 1
[,1] [,2] [,3] [,4]
[1,] 0.01025366 0.134281 -4.292023e+02 4.292023e+02
[2,] 0.13428099 5.066837 -1.220579e+04 1.220579e+04
[3,] -429.20227492 -12205.788493 3.213470e+07 -3.213470e+07
[4,] 429.20227492 12205.788493 -3.213470e+07 3.213470e+07
, , 2
[,1] [,2] [,3] [,4]
[1,] 6.638739e-03 -7.368134e-02 1.881672e+02 -1.881672e+02
[2,] -7.368134e-02 2.913884e+01 -7.487336e+04 7.487336e+04
[3,] 1.881672e+02 -7.487336e+04 2.039656e+08 -2.039656e+08
[4,] -1.881672e+02 7.487336e+04 -2.039656e+08 2.039656e+08
>
> overlap(Pi=c(0.5, 0.5), Mu=Mu, S=S)
$OmegaMap
[,1] [,2]
[1,] 1.0 -1
[2,] 0.5 1
$BarOmega
[1] -0.5
$MaxOmega
[1] 0
$rcMax
[1] 1 1
Почему OmegaMap имеет отрицательное значение?