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.
>