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 that sort of thing doesn't work for you, please
> post an example of what doesn't work ...
> 
>   Ben Bolker

Ben, thanks for the suggestion and it got me going the right direction.
FWIW-I did get an error from your suggestion, but it led me to find that I
had an mcmc.list object and not an mcmc object.  Once I got past that, your
suggestion worked perfectly. Thanks!

______________________________________________
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.

Reply via email to