[R] AIC on GLMM pscl package
Caribu
hugoazina at gmail.com
Fri Mar 4 15:54:27 CET 2011
Hello,
I'm using GLMM on the pscl package and i'm not getting the AIC on the
summary.
The code i'm using is (example) :
mmall3 <-glmmPQL(allclues ~ cycloc + male, data=dados, family=poisson,
random=~1|animal/idfid)
and the results:
Linear mixed-effects model fit by maximum likelihood
Data: dados
AIC BIC logLik
NA NA NA
Random effects:
Formula: ~1 | animal
(Intercept)
StdDev: 0.4235518
Formula: ~1 | idfid %in% animal
(Intercept) Residual
StdDev: 0.947683 1.752526
Variance function:
Structure: fixed weights
Formula: ~invwt
Fixed effects: allclues ~ cycloc + male
Value Std.Error DF t-value p-value
(Intercept) 1.8050720 0.2653779 333 6.801892 0.0000
cycloc 0.0718826 0.0128099 181 5.611469 0.0000
male1 -0.6254748 0.3552453 5 -1.760684 0.1386
Correlation:
(Intr) cycloc
cycloc -0.060
male1 -0.744 -0.008
Standardized Within-Group Residuals:
Min Q1 Med Q3 Max
-2.8968532 -0.6758671 -0.2687695 0.3338190 4.5704000
Number of Observations: 522
Number of Groups:
Am i doing something wrong? Or theres a code to extract the AIC of the
model?
Thanks,
Best regards
Hugo Caribu
--
View this message in context: http://r.789695.n4.nabble.com/AIC-on-GLMM-pscl-package-tp3335371p3335371.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list