Мне нужно нарисовать график биномиального распределения с 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")
и получаю такое:
как можно увидеть с графика с виду все нормально, но это не так :( Мне нужно чтобы у гистограммы было ооочень много столбиков гистограммы, а получается не очень много и я не могу понять как увеличить число попыток и в гистограмме тоже. Может я что-то не так делаю или у меня в корне неправильный подход?
