[R-sig-ME] Gamma-shape parameter extraction
andreu blanco
andreu.blanco at gmail.com
Tue Oct 24 11:30:14 CEST 2017
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]]
More information about the R-sig-mixed-models
mailing list