fargs$optargs$M <- M
mu <- computeMu(fargs$X, fargs$Y, fargs$optargs)
op <- optimParams(K,link,fargs$optargs)
- x_init <- c( rep(1/K,K-1), as.double(mu), rep(0,K) )
+ x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))
do.call(rbind, op$run(x_init))
},
# flexmix