# gerador_binomial
# Autor
# 13/03/2020
# Gerar valores com distribuição binomial com parâmetros n e p
rm(list = ls())

# Semente
set.seed(37)

# ---- valores fixos ----
# tamanho da amostra
m <- 100 
# número de ensaios
n <- 30
# P(sucesso)
p <- 0.25

# vetor de zeros
y <- c(rep(0,m))
# loop para gerar os valores
for(i in 1:m){
  k <- 0
  p0 <- (1-p)^n
  Fy <- p0
  u <- runif(1)
  while(u>Fy){
    py <- (((n-k)*p)/((k+1)*(1-p)))*p0
    p0 <- py
    Fy <- Fy+py
    k <- k+1
  }
  y[i] <- k
}

plot(table(y)/m)
lines(sort(y),dbinom(sort(y),n,p))