Кривая плотности Пуассона, гистограмма и заштрихованная область с ggplot2?

Question

Как сделать:

1. кривая плотности и гистограмма, показывающая распределение Пуассона с лямбда = 2,5; а также
2. кривая плотности с заштрихованной площадью, показывающая P(X >= 4 с лямбда = 2.5)

ось х должна быть от 0 до 10

Answer 1

Распределение Пуассона - это дискретное распределение вероятностей (функция определяется только при целочисленных значениях). Таким образом, вместо линии лучше представить точки с целочисленными значениями. Чтобы раскрасить определенный диапазон под функцию можно использовать geom = "area" а также xlim = c(min(range), max(range):

ggplot(data.frame(x = 0:10), aes(x)) +
  stat_function(geom = "point", n = 11, fun = dpois, args = list(lambda = 2.5)) +
  stat_function(geom = "area", aes(x), n = 7, fun = dpois, args = list(lambda = 2.5), xlim = c(4,10), fill = "lightblue", alpha = 0.5)+
  theme_bw()+
  scale_x_continuous(breaks = 0:10)

Если n аргумент в stat_function не соответствует количеству целочисленных значений в диапазоне, который график будет выглядеть в стиле фанк.