- op <- optimParams(K,link,fargs$optargs)
- x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))
- do.call(rbind, op$run(x_init))
- },
- # flexmix
- function(fargs) {
- library(flexmix)
- source("../patch_Bettina/FLXMRglm.R")
- K <- fargs$optargs$K
- dat <- as.data.frame( cbind(fargs$Y,fargs$X) )
- 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] )
- } )
- rbind(p, beta_b[2:nrow(beta_b),], beta_b[1,])
- } ),
+ res2 <- NULL
+ tryCatch({
+ op <- optimParams(K,link,fargs$optargs)
+ x_init <- list(p=rep(1/K,K-1), beta=mu, b=rep(0,K))
+ res2 <- do.call(rbind, op$run(x_init))
+ }, error = function(e) {
+ res2 <- NA
+ })
+ res2
+ }
+# ,
+# # flexmix
+# function(fargs) {
+# library(flexmix)
+# source("../patch_Bettina/FLXMRglm.R")
+# K <- fargs$optargs$K
+# 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
+# }
+ ),