有没有办法在R（Phyr和MuMin软件包）中挖掘PGLMM_比较模型？

我正在做一个比较分析，我的反应变量是0或1，因此我需要用二项误差分布做一个系统发育校正分析。我使用phyr包（）中的PGLMM_compare函数创建了一个包含所有变量的完整模型，但MuMin不支持将此输出作为“全局模型”，因此我无法对其进行挖掘。我正在寻找一种方法，以找到最好的模型，并可能从这些执行模型平均，但似乎这些软件包是不兼容的。手工创建所有模型是很困难的，因为我有8个解释变量。有没有办法挖掘出具有二项式误差结构的系统发育模型？提前感谢。

您需要至少为

挖泥船

和

模型.avg实施以下方法才能与
pglmm\u compare
一起使用：

nobs.pglmm_比较（对象，…）
logLik.pglmm_比较（对象，…）
系数pglmm_比较（对象，…）
可系数。pglmm_比较（模型，…）