R 如何编程阻塞日志模型

我在R中编程了一个随机变量离散时间危险模型，如下所示

(logit.full. <-
   glmer(event ~ 
                  + a1
                  + a2
                  + a3
                  + obsnum1 + obsnum2 + obsnum3
                  + (1 + obsnum1 + obsnum2 + obsnum3 | country_cluster),
         family=binomial("logit"), data=data.final))

（logit.full.见二项式：
“二项式”系列链接“logit”、“probit”、“cauchit”（分别对应于logistic、normal和Cauchy CDF）“log”和“cloglog”（补充日志）”

因此，您可以使用：

glmer(..., family=binomial("cloglog"), ...)