[R] Plot multiple columns

[Noah Silverman]

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