-# ,
-# # flexmix
-# function(fargs) {
-# library(flexmix)
-# source("../patch_Bettina/FLXMRglm.R")
-# K <- ncol(fargs$beta)
-# dat <- as.data.frame( cbind(fargs$Y,fargs$X) )
-# res2 <- NULL
-# tryCatch({
-# fm <- flexmix( cbind(V1, 1-V1) ~ .-V1, data=dat, k=K,
-# model = FLXMRglm(family = binomial(link = link)) )
-# p <- mean(fm@posterior[["scaled"]][,1])
-# out <- refit(fm)
-# beta_b <- sapply( seq_len(K), function(i) {
-# as.double( out@components[[1]][[i]][,1] )
-# } )
-# res2 <- rbind(p, beta_b[2:nrow(beta_b),], beta_b[1,])
-# }, error = function(e) {
-# res2 <- NA
-# })
-# res2
-# }