#  Examples of Poisson Distributions
#  Edps 587
#  C.J.Anderson
#
#  
#

# create data
y1 <- seq(from=0, to=20)
y2 <- seq(from=0, to=20)
y10 <- seq(from=0, to=30)
y25 <- seq(from=0, to=50)


# find density values
d1 <- dpois(y1,1)           # mean=1
d2 <- dpois(y2,2)           #      2
d10<- dpois(y15,10)         #     10
d25 <- dpois(y25,25)        #     25   

# graph
par(mfrow=c(2,2))

plot(y1,d1,type="h",main="Poisson:  lambda=1")
plot(y2,d2,type="h",main="Poisson:  lambda=2")
plot(y10,d10,type="h",main="Poisson:  lambda=10")
plot(y25,d25,type="h",main="Poisson:  lambda=25")