我不确定如何生成使用GLMM时提示的置信区间 .

通常,如果它不是二项式,我会做这样的事情:

library(boot)
library(lme4)
library(dplyr)
dat <- data.frame(x = runif(100, -2,2),ind = gl(n = 10, k = 10))
dat$y <- 1 + 2 * dat$x + rnorm(10, 0, 1.2)[dat$ind] + rnorm(100, 0, 0.5)
m <- lmer(y ~ x + (1|ind), dat)

b_par <- bootMer(x = m, FUN = fixef, nsim = 200)
boot.ci(b_par, type = "perc", index = 1)

但是当 y 是二项式时,我不知道如何生成上面提到的 dat$y .

谢谢