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
この値を抽出する関数は見つかりませんでした。
前もって感謝します!