[R-sig-ME] Comparison of crossed ranom effects: lmer vs. MCMCglmm

Linus Holtermann holtermann at hwwi.org
Mon Jan 19 18:39:52 CET 2015


Hello,

I read that lmer can handle independent (often labelled as crossed) random effets in mixed models. It seems to be possible with MCMCglmm as long as groups for the random effects are uniquely labelled. I use the "Penicllin" data in the lme4-package to compare both approaches:

library(lme4)
library(MCMCglmm)

str(Penicillin)
attach(Penicillin)

ml <- lmer(diameter~ 1 + (1|plate)+ (1|sample))
summary(ml)

mcmc <- MCMCglmm(diameter~ 1, random=~ plate + sample,verbose=F, nitt=110000,burn=10000,thin=10,data=Penicillin)
summary(mcmc)

Why are the result for the plate-variance differ by a large amount? Is it because MCMCglmm applies Gibbs sampling? Or is MCMCglmm doing something else here, instead of fitting independent random effects?


Best regards,


Linus Holtermann
Hamburgisches WeltWirtschaftsInstitut gemeinnützige GmbH (HWWI)
Heimhuder Straße 71
20148 Hamburg
Tel +49-(0)40-340576-336
Fax+49-(0)40-340576-776
Internet: www.hwwi.org
Email: holtermann at hwwi.org
 
Amtsgericht Hamburg HRB 94303
Geschäftsführer: PD Dr. Christian Growitsch | Prof. Dr. Henning Vöpel
Prokura: Dipl. Kauffrau Alexis Malchin
Umsatzsteuer-ID: DE 241849425


More information about the R-sig-mixed-models mailing list