[R] lmer, mcmcsamp, coda, HPDinterval

[Gabor Grothendieck]

Detach the coda library.  Its HPDinterval is conflicting with the
HPDinterval in lme4.  You should have gotten a warning message about
the conflict when you tried to load coda.

On Sat, Jan 30, 2010 at 9:34 PM, Doug Adams <fog0 at gmx.com> wrote:
>
> Hi,
>
> I've got a linear mixed model created using lmer:
>
> A6mlm <- lmer(Score ~ division + (1|school), data=Age6m)
>
> (To those of you to whom this model looks familiar, thanks for your patience
> with this & my other questions.)  Anyway, I was trying this to look at the
> significance of my fixed effects:
>
> A6post <- mcmcsamp(A6mlm, 50000)
> library(coda)
> HPDinterval(A6post)
>
> ..but I got this message:
> "no applicable method for 'HPDinterval' applied to an object of class
> "merMCMC""
>
> Should I be coercing A6post to another type, or am I missing other steps
> altogether?
> Thanks   :)
>
> Doug Adams
>
> -----
> Doug Adams
> MStat Student
> University of Utah
> --
> View this message in context: http://n4.nabble.com/lmer-mcmcsamp-coda-HPDinterval-tp1457803p1457803.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>