[ADMB Users] ICC calculation in glmmADMB R package
Rafael Mares
crm53 at cam.ac.uk
Tue Sep 13 03:20:04 PDT 2011
Is test2$S what you are looking for? I think it is similar to (or the
same as?) VarCorr() but for glmm.ADMB models.
--
Rafael Mares
Large Animal Research Group (LARG)
Department of Zoology
University of Cambridge
Downing Street
Cambridge
CB2 3EJ
On 13 September 2011 09:33, Sen Li <sen.li at uclouvain.be> wrote:
> 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
>
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users
>
>
More information about the Users
mailing list