error checking, try again

#' @export 

find_markers_glmm_single_gene <- function(dge_formula, design, y, main_effect, nsim) {

    ## Estimate model 

    tryCatch({

        glmer_res <- lme4::glmer(

            dge_formula, cbind(design, y), 'poisson',

            control = lme4::glmerControl(

        res$dge$sigma <- as.numeric(apply(as.data.frame(arm::sim(glmer_res, n.sims=nsim)@ranef[main_effect]), 2, sd))

        return(res)

    }, error = function(e) {

        return(NA)

    })

}