[R-sig-ME] Bug in mcmcsamp

[Jarrod Hadfield]

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