Displaying 2 results from an estimated 2 matches for "ratdrink".
Did you mean:
randrik
2006 Oct 21
0
Constructing predictions from HPDinterval() after lmer()
...ng up on http://finzi.psych.upenn.edu/R/Rhelp02a/archive/
81159.html where Douglas Bates gives a helpful application of lmer()
to data(sleepstudy, package = 'lme4'), I need a bit more help in
order to plot the correct confidence intervals of a designed
experiment such as:
> data(ratdrink, package = 'faraway')
I follow the steps Douglas took in his example:
> summary( rd.er <- lmer(wt ~ weeks*treat + (1 | subject), data =
ratdrink) )
Linear mixed-effects model fit by REML
Formula: wt ~ weeks * treat + (1 | subject)
Data: ratdrink
AIC BIC logLik MLdevian...
2008 Jun 30
1
Coda not providing summary on mcmc object
The object is a mcmc sample from lmer. I am using R v2.7.1. Please let me
know what additional information I can provide, hopefully I am just making a
simple mistake. Thanks in advance!
> data(ratdrink, package = 'faraway')
> rd.er <- lmer(wt ~ weeks*treat + (1 | subject), data = ratdrink)
> rd.mc <- mcmcsamp(rd.er, 10000)
> library(coda)
Loading required package: lattice
> summary(rd.mc)
Length Class Mode
1 merMCMC S4
> HPDinterval(rd.mc)
Error in...