Probabilidade condicional do nascimento

probabilidade_nascimento

prob_diaria <- function(dias_passados){
 
  dia_quantile <- seq(from = -2, 
                      to = 2, 
                      length.out = 30)
  prob <- pnorm(dia_quantile)
  prob <- c(prob[1], prob[2:29] - prob[1:28], prob[1])
  d <- data.frame(dia = (40 * 7 - 15):(40 * 7 + 14), 
                  prob)
  
  d <- d[(dias_passados + 1):nrow(d),]
  d$prob <- d$prob / sum(d$prob)
  
  return(d)
} 

prob_diaria(7)