Simulação simples de P valor

p_valor_vacina.R

doentes_no_placebo = 1
doentes_na_vacina = 0
total_de_pessoas_por_grupo = 2

chance_de_covid = doentes_no_placebo/total_de_pessoas_por_grupo



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
  
   if (casos_grupo_1 == 0 ){
    # É um pouco migué, mas vamos considerar que se ninguém 
    # pegou no grupo, é 100% de eficácia, se não, não!
    resultado_dessa_simulacao = "foi 100%"
  } else {
    resultado_dessa_simulacao = "foi menos que 100%"
   }

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


print(table(resultado_total_das_simulações))