[R] Plotting multiple curves with lattice graphs

Karen Chiswell karenc2204 at yahoo.com
Thu Apr 5 19:12:44 CEST 2007


Thanks for two nice & simple solutions.

I did not know about recordPlot() or replayPlot(). 
Unfortunately these do not give me quite what I want
because I need (for use in another function) an object
of class "lattice".  I should have mentioned this part
earlier!

However, your first solution using xyplot() works and
I think it is a bit better than what I'm currently
using.

Thanks!
Karen

--- Gabor Grothendieck <ggrothendieck at gmail.com>
wrote:

> Try this:
> 
> x <- 1:100
> xyplot(dnorm(x, 50, 10) + dnorm(x, 55, 12) +
> dnorm(x, 60, 15) ~ x, type = "l")
> 
> By the way, you can save classic graphics, at least
> on windows, like this:
> 
> ### your code
> parms <-
>
data.frame(ID=c(1,2,3),mu=c(50,55,60),sigma=c(10,12,15))
>
curve(dnorm(x,mean=parms$mu[1],sd=parms$sigma[1]),from=0,
> to=150, ylab="density", col="red")
>
curve(dnorm(x,mean=parms$mu[2],sd=parms$sigma[2]),from=0,
> to=150, ylab="density", col="blue", add=TRUE)
>
curve(dnorm(x,mean=parms$mu[3],sd=parms$sigma[3]),from=0,
> to=150, ylab="density", col="forestgreen", add=TRUE)
> ###
> 
> myplot <- recordPlot()
> title("Title A")
> 
> # replay the plot with Title B
> replayPlot(myplot)
> title("Title B")
> 
> 
> 
> On 4/5/07, Karen Chiswell <karenc2204 at yahoo.com>
> wrote:
> > Hi List,
> >
> > I would like to plot multiple curves (parametric
> > density curves) in one plot.
> >
> > For example:
> > # parameters for three normal density curves
> > parms =
> >
>
data.frame(ID=c(1,2,3),mu=c(50,55,60),sigma=c(10,12,15))
> >
> > # I can easily draw three normal density curves
> using
> > curve():
> >
>
curve(dnorm(x,mean=parms$mu[1],sd=parms$sigma[1]),from=0,
> > to=150, ylab="density", col="red")
> >
>
curve(dnorm(x,mean=parms$mu[2],sd=parms$sigma[2]),from=0,
> > to=150, ylab="density", col="blue", add=TRUE)
> >
>
curve(dnorm(x,mean=parms$mu[3],sd=parms$sigma[3]),from=0,
> > to=150, ylab="density", col="forestgreen",
> add=TRUE)
> >
> > but I would like to use lattice graphics (even
> though
> > I only want one panel) because I want to save the
> > resulting plot in an object.
> >
> > I could create a dataframe that contains the
> values I
> > want to plot, and then plot using xyplot().  For
> > example,
> >
> > x = seq(from=0, to=150, length=100)
> > plot.df = data.frame(ID=rep(parms$ID,each=100),
> >
>
x=rep(x,3),mu=rep(parms$mu,each=100),sigma=rep(parms$sigma,each=100))
> > plot.df$d =
> > dnorm(plot.df$x,mean=plot.df$mu,sd=plot.df$sigma)
> >
> > require(lattice)
> > p=xyplot(d~x, data=plot.df, groups=ID, type="l")
> >
> > However...  is there a more elegant/efficient way
> to
> > do this?  I have tried using the panel.curve panel
> > function (with xyplot() and histogram()), but I am
> not
> > very experienced with using panel functions (or R
> > actually), and I haven't had any success with this
> > yet.
> >
> > Many thanks for any suggestions,
> > Karen
> >
> > -------------------------------
> > Karen Chiswell
> > Department of Statistics
> > North Carolina State University
> > -------------------------------
> >
> >
> >
> >
> >
>
____________________________________________________________________________________
> > It's here! Your new message!
> > Get new email alerts with the free Yahoo! Toolbar.
> >
> > ______________________________________________
> > R-help at stat.math.ethz.ch 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.
> >
> 



 
____________________________________________________________________________________
The fish are biting.



More information about the R-help mailing list