Question

我正在尝试按照Anderson等人的描述制作一张AIC表。（Suggestions for Presenting the Results of Data Analyses. J.Wildl Manage. 65(3): 2001）。为此，我一直在ICtab包中使用bbmle。

aic2 <- ICtab(du1500o,du1500a,du1500b,du1500c,du1500d,du1500e,du1500f,
              du1500g,du1500h,du1500i,du1500j,du1500k,du1500l,
              type="AICc", weights=T,delta=T,sort=T,base=T,nobs=87)

这给了我一个很好的小排名表。

        AICc  df dAICc weight
du1500a 402.4 6    0.0 0.4201
du1500k 403.4 6    1.0 0.2580
du1500j 404.4 6    2.0 0.1520
du1500f 405.6 6    3.2 0.0842
du1500g 406.8 6    4.4 0.0459
du1500e 408.7 6    6.3 0.0176
du1500d 410.2 6    7.8 0.0084
du1500l 410.7 6    8.3 0.0065
du1500i 412.3 6    9.9 0.0030
du1500o 412.4 4   10.0 0.0029
du1500c 415.2 6   12.8 <0.001
du1500h 416.6 6   14.2 <0.001
du1500b 416.6 6   14.2 <0.001

我想将每个模型的对数似然性添加到此表中。我可以使用logLik函数一次提取一个对数似然。但是，我无法找到更有效的方法从所有模型中提取LL。

head(logLik(du1500a))
[1] -194.6726

我应该将aic2变成数据框吗？怎么样？
如果是这样，我怎样才能有效地将每个模型的LL附加到该数据框？

非常感谢，纳瓦

将多个glmer模型中的logLik提取到数据框中

0 个答案: