# Function to run multilevel logistic regression ----
glmmixed <- function(df.in, dependent, explanatory, random_effect){
lme4::glmer(paste0(dependent, "~", paste(explanatory, collapse="+"), " + (1|", random_effect, ")"),
data=df.in, family="binomial", control=lme4::glmerControl(optimizer="bobyqa",
optCtrl=list(maxfun=200000)))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.