[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