[R-sig-ME] mcmcsamp can't handle random effects syntax

Nicholas Lewin-Koh lewin-koh.nicholas at gene.com
Wed Sep 17 18:05:44 CEST 2008


Hi,
This may be well known, but it was mucking me up.

> fm1 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy)
> ttt<-mcmcsamp(fm1,n=1000)
> fm1 <- lmer(Reaction ~ Days + (1+Days|Subject), sleepstudy)
> ttt<-mcmcsamp(fm1,n=1000)
Error in .local(object, n, verbose, ...) : 
  Code for non-trivial theta_T not yet written
> fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
> ttt<-mcmcsamp(fm1,n=1000)
Error in .local(object, n, verbose, ...) : 
  Code for non-trivial theta_T not yet written

These are all the same model, right?

Thanks
Nicholas



sessionInfo() 
R version 2.7.2 Patched (2008-08-26 r46446) 
i386-pc-mingw32 

locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

attached base packages:
[1] splines   stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
[1] XML_1.96-0         lme4_0.999375-26   Matrix_0.999375-13 lattice_0.17-14


loaded via a namespace (and not attached):
[1] grid_2.7.2


 "Seek the company of those who seek the truth, and run away from those who
have found it." - Vaclav Havel




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