r - 如何编写堵塞日志模型

Question

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

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

其中 obsnum 1-3 是基线风险函数，a1-a3 是随机效应。

现在我想用一个互补的log-log（clog-log）链接来计算模型，但到目前为止我还没有发现如何正确指定链接。现在有人吗，怎么做？

任何帮助将不胜感激！

Answer 1

参见 ?binomial：“'binomial' 系列链接 'logit'、'probit'、'cauchit'，（分别对应于逻辑、正常和 Cauchy CDF）'log' 和 'cloglog'（互补 log-log）”

所以你会使用：

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