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