我想知道是否有人知道从 R 中 glmmadmb 的输出中提取色散参数估计值的方法。我使用的是负二项式模型,并且想将此代码用于几个不同的物种而不必去并为代码的其余部分手动提取此值。
这是我的输出示例:
> summary(mod1)
Call:
glmmadmb(formula = species ~ (1 | year) + (1 | site), data = cs,
family = "nbinom2", link = "log")
AIC: 8131.7
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 4.05 0.19 21.3 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Number of observations: total=798, year=53, site=15
Random effect variance(s):
Group=year
Variance StdDev
(Intercept) 0.1925 0.4388
Group=site
Variance StdDev
(Intercept) 0.4748 0.689
Negative binomial dispersion parameter: 1.7211 (std. err.: 0.088936)
Log-likelihood: -4061.86
我还没有找到提取这个值的函数。
先感谢您!