#' 1/2, matrix(c(1,-2,3,1),ncol=2), c(0,0), "logit")
#' μ <- computeMu(io$X, io$Y, list(K=2))
#' o <- optimParams(io$X, io$Y, 2, "logit")
#' 1/2, matrix(c(1,-2,3,1),ncol=2), c(0,0), "logit")
#' μ <- computeMu(io$X, io$Y, list(K=2))
#' o <- optimParams(io$X, io$Y, 2, "logit")