0

我正在尝试按照 Anderson 等人的描述制作 AIC 表。（提出数据分析结果的建议。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
  1. 我应该aic2变成数据框吗？如何？
  2. 如果是这样，我怎样才能有效地将每个模型的 LL 附加到该数据帧？

非常感谢，纳瓦

4

0 回答 0