[ADMB Users] ICC calculation in glmmADMB R package

Sen Li sen.li at uclouvain.be
Tue Sep 13 01:33:07 PDT 2011


Dear Ben,

Thanks for the reply. I tried VarCorr() before, but it turned out to be 
inapplicable with gimm.ADMB objects.

Here follows the outputs of my model, of VarCorr(), and of sessionInfo():

 > print(test2,sd_S_print=T)

GLMM's in R powered by AD Model Builder:

   Family: nbinom
   alpha = 1.0808

Fixed effects:
   Log-likelihood: -690.741
   AIC: 1395.482
   Formula: y ~ x1 + x2 + x3
  (Intercept)           x1                x2                    x3
    11.048000    -0.025435   -45.852000    -0.061126

Random effects:
   Grouping factor: Scode
   Formula: ~1
Structure: Diagonal matrix
(Intercept)
    1.275421

Covariance matrix of random effects vector (left) and corresponding 
standard deviations (right):

             (Intercept)     (Intercept)
(Intercept)      1.6267         0.77551

Note: The diagonal elements of the above left matrix are variances, NOT 
std's.

Number of Observations: 125
Number of Groups: 51
 > VarCorr(test2)
Erreur dans UseMethod("VarCorr") :
   pas de méthode pour 'VarCorr' applicable pour un objet de classe 
"glmm.admb" ## the method is not applicable to an object of class 
"glmm.admb"##
 > sessionInfo()
R version 2.13.0 (2011-04-13)
Platform: x86_64-pc-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=French_Belgium.1252  LC_CTYPE=French_Belgium.1252
[3] LC_MONETARY=French_Belgium.1252 LC_NUMERIC=C
[5] LC_TIME=French_Belgium.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] nlme_3.1-101   glmmADMB_0.5-2 MASS_7.3-12

loaded via a namespace (and not attached):
[1] grid_2.13.0     lattice_0.19-23 tools_2.13.0
 >

I'd appreciate any further help.

Best,
Sen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20110913/ec0b256f/attachment.html>


More information about the Users mailing list