如何在“GLMMadaptive”的“mixed_model”函数的“random”参数中指定两个随机截距?

How to specify two random intercepts in "random" argument of "mixed_model" function of "GLMMadaptive"?

提问人:Farzin Shamloo 提问时间:10/29/2023 更新时间:10/29/2023 访问量:22

问:

我想获得具有两个随机截距的二项式模型的边际估计及其置信区间,但我不知道如何指定“GLMMadaptive”的“mixed_model”函数的随机参数。 这是我将它用于具有一个随机截距的模型的方式:

m3 <- mixed_model(fixed = cbind(resp,N-resp)~1,random=~1|Subject, data = df,
                  family = binomial(link = "probit"))
M_adp <- marginal_coefs(m3,std_errors = TRUE, cores = 5)
pnorm(confint(M_adp))

如果我想添加另一个随机截距项,我应该如何编码?

我尝试了不同的符号(例如,),但它不起作用。random = list(subjid = ~ 1, condition = ~1)

置信区间 边际效应

评论


答: 暂无答案