EMCluster icon indicating copy to clipboard operation
EMCluster copied to clipboard

simple.init() runs forever with small matrices

Open torockel opened this issue 3 years ago • 0 comments

Hi everyone,

simple.init() seem to run forever with some small matrices. For example:

set.seed(123)
ds <- matrix(rnorm(50), ncol = 5)
EMCluster::simple.init(ds, nclass = 2)

never finished on my machine. But,

set.seed(123)
ds <- matrix(rnorm(500), ncol = 5)
EMCluster::simple.init(ds, nclass = 2)
#> Method: 
#>  n = , p = 0.82147030.1785297, nclass = , flag = , logL = .
#>  conv.iter = , conv.eps = ,
#> nc: 
#> NULL
#> pi: 
#> [1] 0.8215 0.1785

Created on 2022-03-22 by the reprex package (v2.0.1)

runs in under a second.

torockel avatar Mar 22 '22 13:03 torockel