R MixSim Отрицательное Вероятность перекрытия отрицательна - PullRequest
0 голосов
/ 04 февраля 2020

Я пытаюсь оценить вероятность того, что наблюдение будет классифицировано как одна из двух возможных многомерных нормалей. Используя пакет 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 имеет отрицательное значение?

...