# vector de colores "colors"
colors <- c("Red", "Blue", "Gold", "Black", "Pink", "Green")

# declara una lista que tenga los valores de la distribucion de poisson que vamos a crear
poisson.dist <- list()

# Guarda el vector de la distribucion de cada valor
a <- c(1, 2, 3, 4, 5, 6) 
for (i in 1:6) {
  poisson.dist[[i]] <- c(dpois(0:20, i)) 
}

# Grafica cada vector de la lista usando los colores para representar cada valor
plot(unlist(poisson.dist[1]), type = "o", xlab="y", ylab = "P(y)",
     col = colors[i])
for (i in 1:6) {
  lines(unlist(poisson.dist[i]), type = "o", col = colors[i])
}
# Añade leyenda al grafico creado
legend("topright", legend = a, inset = 0.02, cex = 0.8, 
       horiz = TRUE, fill = colors, title = "Valores de lambda (la media)")

dpois(7,3)
dpois(0,3)
dpois(1,3)
dpois(2,3)
dpois(0,3) + dpois(1,3) + dpois(2,3)
ppois(2,3)