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 <f...@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@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. > ______________________________________________ R-help@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.
