Как нарисовать график биномиального распределения в R?

Question

[Андрей Горошко]

Мне нужно нарисовать график биномиального распределения с 100000 попыток. График должен состоять из самого графика линии и гистограммы. Я делаю так:

library(tidyverse)
summarize.binomial.observations <- function(trials, theta, observations) {
  sample.space <- c(1,0)
  results <- 1:observations %>%
  map_int(function(x)
    as.integer(
      sum(
        sample(sample.space, size = trials, replace = TRUE, prob = c(theta, 1 - theta))
      )
    )
  )
  return(results)
}
summary <- summarize.binomial.observations(80, 0.6, 5000)
resulting.df <- data.frame(flips <- summary) # ggplot only works with data frames
names(resulting.df) <- c("flips")
#binPDF <- data.frame(x<-1:80, y<-dbinom(1:80, size=80, prob = 0.6))
#names(binPDF) <- c("flips", "prob")
ggplot(resulting.df) + 
  geom_histogram(breaks=seq(1, 80, 2), aes(x=flips, y=..density..), position="identity",,colour = 'blue', fill = 'white') + 
  geom_density(aes(x=flips, y = ..density..), colour="red")

и получаю такое:



как можно увидеть с графика с виду все нормально, но это не так :( Мне нужно чтобы у гистограммы было ооочень много столбиков гистограммы, а получается не очень много и я не могу понять как увеличить число попыток и в гистограмме тоже. Может я что-то не так делаю или у меня в корне неправильный подход?

Answer 1

[zexer]

Чтобы увеличить количество столбцов в гистограмме распределения (для ggplot), используйте параметр bins внутри geom_histogram.