Glucose1 <- glmer(GlucoseYN ~ 1 + (1|Subject), data=LongFormat, family = "binomial", glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 10000000)))

getME(Glucose1, c("theta","beta"))
# theta = 3.276187, beta = -0.3515557

Glucose1b <- glmer(GlucoseYN ~ 1 + (1|Subject), data=LongFormat, family = "binomial", start = 3.276187, glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 10000000)))