Message-ID: <loom.20110120T150247-276@post.gmane.org>
Date: 2011-01-20T14:10:08Z
From: Ben Bolker
Subject: MCMC object indexing
Robinson, David G <drobin <at> sandia.gov> writes:
>
> I have an mcmc object and I'm trying to plot the quantiles of the
> variables - and not as a function of the
> iterations as in cumuplot.
>
> I cannot seem to find the right combination of indexing to
> access the variables; after which I'm sure I can
> plot all the statistics I could hope for.
>
> Any hints for accessing the mcmc object would be appreciated.
>
> =Dave
A reproducible example would be very useful. In general
I find that I can just index an mcmc object as though it
were a matrix, e.g.:
z <- cbind(a=runif(20),b=runif(20),c=runif(20))
library(coda)
m <- mcmc(data=z)
t(apply(m,2,quantile,c(0.025,0.975)))
If that sort of thing doesn't work for you, please
post an example of what doesn't work ...
Ben Bolker