Hi, I used the term "run", as each iteration of the Gibbs sampler produces 22 variables (coefficients for Beta in a regression model)
The example wont work On 6/1/10 5:54 AM, Ben Bolker wrote: > Noah Silverman <noah <at> smartmediacorp.com> writes: > > >> I'm running a long MCMC chain that is generating samples for 22 variables. >> >> I have each run of the chain as a row in a matrix. >> > >> So: Chain[,1] is the column with all the samples for variable one. >> Chain[,2] is the column with all the samples for variable 2, etc. >> > The previous 2 paragraphs seemed contradictory until I realized > that in the first paragraph you are using "run" to mean what I would > usually call a "sample" ... > > >> I'd like to fit all 22 on a single page to print a nice summary. It is >> OK if the graphs are small, I just need to show the overall shape and >> convergence. >> > How about for example > > >> x <- matrix(runif(22000),ncol=22) >> library(coda) >> m <- as.mcmc(x) >> xyplot(m) >> xyplot(m,layout=c(4,6)) >> > ______________________________________________ > 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. > ______________________________________________ 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.