- op <- optimParams(fargs$X, fargs$Y, K, fargs$link, fargs$M)
- x_init <- list(p=rep(1/K,K-1), beta=fargs$mu, b=rep(0,K))
+ mu <- computeMu(fargs$X, fargs$Y, list(K=K, M=fargs$M))
+ op <- optimParams(fargs$X, fargs$Y, K, fargs$link, fargs$M, 1)
+ x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))