Displaying 1 result from an estimated 1 matches for "samp0".
Did you mean:
samp
2011 Feb 19
0
lmer, MCMCsamp and ranef samples?
I really hope sombody could help me with the following,
I'm having problems accessing the random effect samples following the
example on MCMCsamp:
(fm1 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy))
set.seed(101); samp0 <- mcmcsamp(fm1, n = 1000, saveb=TRUE)
str(samp0)
Formal class 'merMCMC' [package "lme4"] with 9 slots
..@ Gp : int [1:3] 0 18 36
..@ ST : num [1:2, 1:1000] 0.98 0.234 1.097 0.258 0.915 ...
..@ call : language lmer(formula = Reaction ~ Days + (1 |
Subj...