An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110119/4900f07d/attachment.pl>
MCMC object indexing
2 messages · Robinson, David G, Ben Bolker
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