[R] multiple fitted curves plot

Douglas Bates bates at stat.wisc.edu
Sat Sep 8 17:41:06 CEST 2001


Stewart Schultz <schultz at fig.cox.miami.edu> writes:

> How would I plot several (24) fitted curves (nonlinear least squares)
> within a single plot?  

Are you sure you want to do that?  Overlaying 24 curves on a single
plot tends to make for a busy plot.  You may want to consider a
trellis type of plot instead.  See xyplot from the lattice package.

> To get the parameters I use:
> 
> resff <- list()
> for (s in levels(PairID)) {
>       resff[[s]] <- nls(Photo ~ SSasymp(Ci, Asym, lrc, c0),
>       subset = (PairID == s))
> }

An alternative would be to use the nlsList function from the nlme
package.

If you really do want to plot all the fitted nls models on a single
plot, it is best to use predict on the fitted model with a closely
spaced set of Ci values that cover the range of the observed data.  It
could look like


plot(Ci, Photo, type = "n")  # sets up the plotting scales
crange <- range(Ci)
cseq <- seq(crange[1], crange[2], len = 100)
for (s in levels(PairID)) {
  lines(cseq, predict(resff[[s]], list(Ci = cseq)))
}

> Is it also possible to plot two different groups (male, female) of curves
> in the same plot, with different style lines?

Yes.

>From your questions it seems that you may want to consider some of the
techniques available in the nlme package.  The book

@Book{pinh:bate:2000,
  author =	 {Jos\'{e} C. Pinheiro and Douglas M. Bates},
  title = 	 {Mixed-Effects Models in \textsf{S} and \textsf{S-PLUS}},
  publisher = 	 {Springer},
  year = 	 2000,
  series =	 {Statistics and Computing}
}

has several examples of these techniques.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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