Hi, I am dealing with hurdle models for which I need to extract the dispersion parameter of my gamma model, which I did with glmmadmb. However I am not able to extract that value from the summary. Of course I can copy paste, but I rather want to know how to do it. Thanks! That's my summary:
summary(zag1)
Call:
glmmadmb(formula = Biomass ~ Protection + Exposure + Protection:Exposure +
(1 | Protection:Location), data = Aarmata.pos, family = "Gamma",
link = "log")
AIC: 644.9
Coefficients:
Estimate Std. Error z value
Pr(>|z|)
(Intercept) 4.184 0.452 9.25
<2e-16 ***
ProtectionProtected 0.201 0.647 0.31
0.756
ExposureSemiexposed -0.609 0.349 -1.75
0.081 .
ProtectionProtected:ExposureSemiexposed -0.674 0.575 -1.17
0.241
---
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
Number of observations: total=61, Protection:Location=8
Random effect variance(s):
Group=Protection:Location
Variance StdDev
(Intercept) 0.6592 0.8119
Gamma shape parameter: 1.4473 (std. err.: 0.25759)
Log-likelihood: -316.472
Andreu Blanco Cartagena Si no ?s imprescindible, no imprimeixis aquest e-mail. Estalviar paper ajuda a protegir el medi ambient. Si no es imprescindible, no imprimas este e-mail. Ahorrar papel ayuda a proteger el medio ambiente. [[alternative HTML version deleted]]