[R-sig-ME] observation level random effects; estimated variance of variance component estimates
j.hadfield at ed.ac.uk
Fri Feb 20 13:10:23 CET 2009
I'm not sure how to fit the model you want in lmer, but hopefully
someone will have a better idea. However, if its of interest or use
I've released the package MCMCglmm that fits the model you're after.
By default MCMCglmm fits rand as a residual component:
m1<-MCMCglmm(Seizure ~ Base + Trt*Base + Age + V4, random=~id,
which I believe is the model you want.
gives CI's on the variance components (the second of which is you rand
Hope this helps,
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
More information about the R-sig-mixed-models