[R] lmer, mcmcsamp, coda, HPDinterval

[Doug Adams]

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.