[R-sig-ME] Bug in mcmcsamp
Jarrod Hadfield
j.hadfield at ed.ac.uk
Sun Aug 12 14:03:08 CEST 2007
Hi,
I'm using lme4 0.99875-6 and I think there may be a bug in mcmcsamp.
When I fit 2 random effects, the posterior variances always seem to
be approximately the same irrespective of the true variances (or the
REML estimates of the true variances). For example:
A<-gl(100,10)
B<-gl(100,1,1000)
y<-rnorm(100,0,sqrt(1))[A]+rnorm(100,0,sqrt(3))[B]+rnorm
(1000,0,sqrt(3))
model1<-lmer(y~(1|A)+(1|B))
model1MCMC<-mcmcsamp(model1,n=5000,trans=FALSE)
plot(model1MCMC)
Sorry, if its my mistake.
Jarrod
More information about the R-sig-mixed-models
mailing list