gam plots

I would like postscript gam plots for the two smoothed terms to be 
produced on separate pages, but am having trouble getting these.  I 
can get the two plots output to a single page using

 postscript("gams.ps")
 plot(mymodel.gam, pages=1, se=T)
 dev.off()

I don't want this, but when I try replacing pages=1 with pages=2, I 
still get only one page in the postscript file with the second plot 
overlaid on the first one:

 plot(mymodel.gam, pages=2, se=T)

Press return for next page....

I also noticed that the plot gives tick marks for the values of the 
variable along the horizontal axis, but I would like a rug of *all* 
the values (Splus5 produces a rug using jittering to show the density 
of the observations, but the version of R that I have doesn't do 
this, as far as I can tell).

Simon Wood  snw at st-and.ac.uk  http://www.ruwpa.st-and.ac.uk/simon.html
The Mathematical Institute, North Haugh, St. Andrews, Fife KY16 9SS UK
Direct telephone: (0)1334 463799          Indirect fax: (0)1334 463748 

Dear R users,

Using the library(mgcv) and running R under MacOSX, I have fitted a
generalised additive model with binomial errors in order to check the
linearity of two continuous variables ap2mm and diffdaysm in a glm:

 mymodel.gam <- gam(diedhos~ s(ap2mm) + Dweekm + s(diffdaysm) +

  Dweekm:diffdaysm + ap2mm:Dweekm, binomial)

I would like postscript gam plots for the two smoothed terms to be
produced on separate pages, but am having trouble getting these.  I
can get the two plots output to a single page using

 postscript("gams.ps")
 plot(mymodel.gam, pages=1, se=T)
 dev.off()

I don't want this, but when I try replacing pages=1 with pages=2, I
still get only one page in the postscript file with the second plot
overlaid on the first one:

 plot(mymodel.gam, pages=2, se=T)

Press return for next page....


	I wondered if anyone knows of a way to get each plot on a
separate page, i.e. two postscript files?

I also noticed that the plot gives tick marks for the values of the
variable along the horizontal axis, but I would like a rug of *all*
the values (Splus5 produces a rug using jittering to show the density
of the observations, but the version of R that I have doesn't do
this, as far as I can tell).