Another thought I had was to use Sweave.
I could, in theory, use Sweave to create a plot only 7 inches high, and
then a summary table 4 inches high below it, all on the same page of a
Latex document. Could be messy, but might work.
On 6/1/10 11:18 AM, Dennis Murphy wrote:
>
>
> On Tue, Jun 1, 2010 at 10:51 AM, Noah Silverman
> > wrote:
>
> You are correct,
>
> I initially missed the "as.mcmc" step. Without it, R doesn't want
> to squeeze so many plots onto a page.
>
>
> I've had that problem before with lattice...which makes me wonder
> whether you'll be able to put the plots
> and table on a single page. Here are some haphazard thoughts:
>
> * The gplots package has functions textplot and sinkplot which IIRC
> will allow you to add a table to a plot.
> I think you need to use layout() in base graphics to set it up
> properly, but you'd need to get 22 graphs
> and a table on one graphics page...
> * Given the number of graphs you want per page, perhaps sparkline
> plots might be less wasteful
> of space (see the YaleToolkit package for one approach to sparklines)
> * The problem is the time it would take to coordinate all of this...
> * In a similar vein, there's always grid... :)
>
> Hopefully Deepayan or Felix Andrews can share some of their Lattice
> wisdom and show us how to
> do what you want in their system. That would be the most convenient
> solution.
>
> Dennis
>
>
> Thanks!
>
>
> On 6/1/10 10:48 AM, Dennis Murphy wrote:
>> Hi:
>>
>> Ben's example worked for me as follows:
>>
>> x <- matrix(runif(22000),ncol=22)
>> library(lattice)
>> library(coda)
>> m <- as.mcmc(x)
>> xyplot(m, layout = c(2, 11))
>>
>> layout = uses columns first, then rows...I always forget that
>> :) It also has an optional
>> third argument for pages, which might come in handy here. If you
>> want the order of the
>> panels to be reversed, use the additional option as.table = TRUE.
>>
>> HTH,
>> Dennis
>>
>> On Tue, Jun 1, 2010 at 9:37 AM, Noah Silverman
>> > wrote:
>>
>> 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 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.
>>
>>
>
[[alternative HTML version deleted]]