[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