如何在一个图中绘制多个泊松分布

Question

我想绘制多个泊松（具有不同的Lambda（1:10））

我发现以下函数来绘制图

plot_pois = function(lambda = 5)
{
  plot(0:20, dpois( x=0:20, lambda=lambda ), xlim=c(-2,20))
  normden <- function(x){dnorm(x, mean= lambda, sd=sqrt(lambda))}
  curve(normden, from=-4, to=20, add=TRUE, col=lambda)
}
plot.new()
plot_pois(2)

但是我无法在它上面再画一个泊松。 我试图改变plot ，以points或lines ，但它完全改变剧情。 我还想为不同的lambda值添加包含不同颜色的图例。

如果我可以使用ggplot对其进行ggplot ，那将是一个更好的选择。

Answer 1

在ggplot2中，您可以使用lapply遍历不同的lambda：

library(ggplot2)
lambdas <- c(5, 2)
ggplot(data = data.frame(x = 0:20)) +
  lapply(lambdas, function(l) geom_point(aes(x = x, y = dpois(x, lambda = l), col = factor(l)))) +
  lapply(lambdas, function(l) stat_function(fun = dnorm, args = list(mean = l, sd = sqrt(l)), 
                                aes(x = x, col = factor(l))))

然后可以像往常一样在ggplot2中自定义轴标题和限制，图例标题等。

Answer 2

另一个可能的tidyverse解决方案：

library(tidyverse)

# Build Poisson distributions

p_dat <- map_df(1:10, ~ tibble(
  l = paste(.),
  x = 0:20,
  y = dpois(0:20, .)
))

# Build Normal distributions

n_dat <- map_df(1:10, ~ tibble(
  l = paste(.),
  x = seq(0, 20, by = 0.001),
  y = dnorm(seq(0, 20, by = 0.001), ., sqrt(.))
))

# Use ggplot2 to plot

ggplot(n_dat, aes(x, y, color = factor(l, levels = 1:10))) +
  geom_line() +
  geom_point(data = p_dat, aes(x, y, color = factor(l, levels = 1:10))) +
  labs(color = "Lambda:") +
  theme_minimal()

^{由reprex软件包 （v0.2.1）创建于2019-05-06}