Re: [R] two plots on the same page

From: Steven McKinney <smckinney_at_bccrc.ca>
Date: Fri, 28 Dec 2007 17:09:45 -0800

> -----Original Message-----
> From: r-help-bounces_at_r-project.org on behalf of Maura E Monville
> Sent: Fri 12/28/2007 2:18 PM
> To: r-help_at_stat.math.ethz.ch
> Subject: [R] two plots on the same page
>
> I'd like to know why I cannot get a plot and the QQnorm in the same sheet.
> The commands are simple but:
>
> library(nlme)
> glmod1 <- gls(upfmla,correlation=corAR1(),method="ML")
> summary(glmod1)
> par(mfrow = c(2,1))
> plot(glmod1, main="GLS Residuals vs. GLS Fitted")
> qqnorm(glmod1)
>
> No matter what (I tried different permutations of the plotting commands) the
> second drawing is overlapped to the first one instead of being placed
> underneath it on the same page.
> How come?
>

Because glmod1 is of class "gls"

> class(glmod1)

[1] "gls"

so when you plot it you are invoking
plot.gls()

See the help for function plot.gls() in the nlme library.

> ?plot.gls

reports that the value of a call is

   Value

   a diagnostic Trellis plot.

so you are dealing with Trellis plots handled by the lattice package. Trellis plots do not use the par() settings used by regular plots.

Here's one way to do it. (I substituted an example from nlme since I don't have your data and "upfmla" object.

library(nlme)

glmod1 <- gls(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), Ovary,

           correlation = corAR1(form = ~ 1 | Mare)) summary(glmod1)

plot1 <- plot(glmod1, main="GLS Residuals vs. GLS Fitted") plot2 <- qqnorm(glmod1)
print(plot1, split = c(1, 1, 1, 2), more = TRUE) print(plot2, split = c(1, 2, 1, 2))

Try something like this. The 'split' mechanism works for me and I get both plots on one page.

Remember also that you can save the output of trellis plot commands and print them later. Trellis plots do not always print right away - in many situations you must explicitly wrap your trellis plot command with print() before it will render on a graphics device.

You can read more about putting multiple trellis plots on a page in the help page for
print.trellis
in the lattice package.

Hope this helps.

After you master the lattice package
you can go to that New Years party
and celebrate!

Steven McKinney

Statistician
Molecular Oncology and Breast Cancer Program British Columbia Cancer Research Centre

email: smckinney +at+ bccrc +dot+ ca

tel: 604-675-8000 x7561

BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C.
V5Z 1L3
Canada

> Thank you very much.
> Happy New Year,
> Maura
>
>
>
>
>
> --
> Maura E.M
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.
>
>



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. Received on Sat 29 Dec 2007 - 01:15:27 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 29 Dec 2007 - 04:30:21 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.