p_valor_vacina.R

chance_de_covid = 0.5

doentes_no_placebo = 1
doentes_na_vacina = 0
total_de_pessoas_por_grupo = 2


doentes_totais = doentes_na_vacina + doentes_no_placebo

## Isso só vale se o número de participantes for igual nos dois grupos!
eficácia_no_estudo = 1 - doentes_na_vacina/doentes_no_placebo

ver_se_pegou_covid <- function(chance_de_covid){
    return(runif(1)<=chance_de_covid)
}

resultado_total_das_simulações = c()
for (i in 1:100){
  
  pessoa_1 = ver_se_pegou_covid(chance_de_covid)
  pessoa_2 = ver_se_pegou_covid(chance_de_covid)
  pessoa_3 = ver_se_pegou_covid(chance_de_covid)
  pessoa_4 =  ver_se_pegou_covid(chance_de_covid)
    
  print(pessoa_4)
  casos_grupo_1 = pessoa_1 + pessoa_2
  casos_grupo_2 = pessoa_3 + pessoa_4
  
  # Evitando dividir por 0
  if (casos_grupo_2 == 0){
    resultado_dessa_simulacao = "não foi tão bem"
  } else if(casos_grupo_1/casos_grupo_2 >= eficácia_no_estudo){
    resultado_dessa_simulacao = "foi tão alto ou melhor"
  }

  resultado_total_das_simulações = c(resultado_total_das_simulações, resultado_dessa_simulacao)
}


print(table(resultado_total_das_simulações))