R help - May 2010 - Need help on having multiple distributions in one graph

R-listers:

I have searched the help files and everything I have related to R 
graphics. I cannot find how to graph y against
several distributions on a single graph. Here is code for creating 4 
Poisson distributions with different mean values, although I would 
prefer having it in a loop: The top of the y axis for the first 
distribution, with count of 0, is .6, which is the highest point for 
any  of the distributions.

obs <- 1:20 
y <- obs-1
m0 <- (exp(-.5) * .5^y)/factorial(y)
m1 <- (exp(-1) * 1^y)/factorial(y)
m3 <- (exp(-3) * 3^y)/factorial(y)
m4 <- (exp(-5) * 5^y)/factorial(y)

How do I plot the graph of each distribution on y, all on a single 
graph? I have spent so many hours on this,
which is really quite simple in applications such as Stata. Thanks very 
much for the assistance:

Joseph Hilbe
hilbe at asu.edu  or jhilbe at aol.com

Hi Joseph,

How about this?

matplot(cbind(m0, m1, m3, m4), type = 'l', lty = 1)
legend('topright', paste('m', c(0, 1, 3, 4), sep =
""), lty = 1, col = 1:4)

See ?matplot and ?legend for details.

HTH,
Jorge


On Mon, May 3, 2010 at 6:42 PM, <> wrote:
> R-listers:
>
> I have searched the help files and everything I have related to R graphics.
> I cannot find how to graph y against
> several distributions on a single graph. Here is code for creating 4
> Poisson distributions with different mean values, although I would prefer
> having it in a loop: The top of the y axis for the first distribution, with
> count of 0, is .6, which is the highest point for any  of the
distributions.
>
> obs <- 1:20 y <- obs-1
> m0 <- (exp(-.5) * .5^y)/factorial(y)
> m1 <- (exp(-1) * 1^y)/factorial(y)
> m3 <- (exp(-3) * 3^y)/factorial(y)
> m4 <- (exp(-5) * 5^y)/factorial(y)
>
> How do I plot the graph of each distribution on y, all on a single graph? I
> have spent so many hours on this,
> which is really quite simple in applications such as Stata. Thanks very
> much for the assistance:
>
> Joseph Hilbe
> hilbe@asu.edu  or jhilbe@aol.com
>
> ______________________________________________
> R-help@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.
>
	[[alternative HTML version deleted]]