[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