Re: [R] MCMC object indexing

2011-01-20 Thread Robinson, David G
Stoker, Ben wrote: > 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

Re: [R] MCMC object indexing

2011-01-20 Thread Ben Bolker
Robinson, David G 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

[R] MCMC object indexing

2011-01-19 Thread Robinson, David G
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 f