[R] gam plots
Patty Solomon
psolomon at maths.adelaide.edu.au
Sat Dec 22 06:17:57 CET 2001
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).
I would be very grateful if anyone has any suggestions that would
help sort this out.
Many thanks,
Patty
--
--------------------------------------------------------------------------------
Assoc Prof Patty Solomon phone: (08) 8303 3033
Department of Applied Mathematics and fax: (08) 8303 3696
Centre for the Molecular Genetics of Development
Adelaide University
Adelaide SA 5005
AUSTRALIA.
email: patty.solomon at adelaide.edu.au
http://www.maths.adelaide.edu.au/people/psolomon
--------------------------------------------------------------------------------
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list