Displaying 3 results from an estimated 3 matches for "sbjc".
Did you mean:
objc
2007 May 11
0
incorrect MCMC CIs in pvals.fnc (languageR) ?
...languageR)
fit = lmer(Reaction~Days + (1|Subject) + (0+Days|Subject),
data=sleepstudy)
pvals.fnc(fit)$random
# compare with...
samp = mcmcsamp(fit, n=10000, trans=FALSE)
HPDinterval(samp)
densityplot(samp, plot=F)
# 'pvals.fnc' reports sigma instead of sigma^2, but it looks like the
# Sbjc.(In) and Sbjc.Days are also sqrt compared with the corresponding
# results from HPDinterval and densityplot. Is that intended?
best wishes, Mike Allerhand
2006 Oct 21
0
Constructing predictions from HPDinterval() after lmer()
...46.420404 59.406398
weeks 25.070131 27.930363
treatthiouracil -4.420942 14.009291
treatthyroxine -10.758369 9.435761
weeks:treatthiouracil -11.404620 -7.337025
weeks:treatthyroxine -1.603858 2.842704
log(sigma^2) 3.683085 4.226153
log(sbjc.(In)) 3.633853 4.955867
deviance 965.750351 980.825613
attr(,"Probability")
[1] 0.95
*******************To make sure my request is clear*******************
What I need is the analog of what is produced by all.effects() in the
effects package:
> summary(rd...
2005 Aug 17
4
How to assess significance of random effect in lme4
Dear All,
With kind help from several friends on the list, I am getting close.
Now here are something interesting I just realized: for random
effects, lmer reports standard deviation instead of standard error! Is
there a hidden option that tells lmer to report standard error of
random effects, like most other multilevel or mixed modeling software,
so that we can say something like "randome