3

我已经使用 glmmTMB 运行了一组 23 个模型。（我已经将我的模型设置为list下面的示例代码）

cand.models<-list()
cand.models[[1]]<-glmmTMB(count~depth + slope + SST + (1|individual), family=list(family="truncated_nbinom1", link="log"), data=df)
cand.models[[2]]<-glmmTMB(count~depth + slope + (1|individual), family=list(family="truncated_nbinom1", link="log"), data=df)

我想创建一个汇总表，为包含在cand.models. 我尝试使用broom::glance()，它应该创建一个“单行”摘要，其中包括偏差，以及 AIC 和 BIC 等其他内容。

summ.table<-do.call(rbind, lapply(cand.models, broom::glance))

但是，输出不包括模型偏差！（它只包括 sigma、logLik、AIC、BIC 和 df.residual）。有谁知道它为什么不提供偏差（也许是 glmmTMB 特有的问题？）。或者，有没有人有替代解决方案来提取偏差？

4

1 回答 1

3

跟进评论：“主”版本broom不包括模型的glance方法glmmTMB，这意味着glance对模型不起作用glmmTMB。

broom您正在使用的版本是 Ben Bolker (bbolker) 的分叉版本，他在其中添加了一种适用于模型对象的新方法。但是，此方法不包含在输出中。不过，由于，您可以执行以下操作：glanceglmmTMBglancedeviancedeviance = -2*log(likelihood)

summ.table$deviance = -2 * summ.table$logLik
于 2017-10-22T23:26:15.390 回答