# ejemplo reproducible ----------------------------------------------------
#' seguimos el ejemplo disponible en
#' https://cran.r-project.org/web/packages/pwr/vignettes/pwr-vignette.html
#' para la función "pwr.2p.test" para
#' test de igualdad de proporciones para dos muestras iguales
#' muestreo aleatorio de
#' hombres y mujeres de
#' dentro de una universidad
#' preguntamos: ¿consume alcohol al menos una vez a la semana?
#'
#' hipótesis nula:
#' no hay diferencia en la proporción de personas que responden "sí".
#' hipótesis alternativa:
#' sí hay diferencia en la proporción de personas que responden "sí".
#'
#' esta es una hipótesis alternativa de "dos-colas"
#' un sexo (asignado al nacer) tiene una proporción más alta pero no sabemos cuál.
#'
#' si quisieramos detectar una _diferencia tan pequeña como 5%_.
#' ¿cuántos estudiantes necesitamos muestrean en _cada grupo_
#' si queremos un _poder de 80%_ y
#' un _nivel de significancia de de 0.05_?
#'
#' si pensamos que en _un grupo la proporción es de 55% y en el otro de 50%_:
pwr.2p.test(h=ES.h(p1 = 0.55, p2 = 0.50),
power = 0.8,
sig.level = 0.05,
alternative = "two.sided")
#' OJO: el tamaño de muestra es por grupo
#' Necesitamos muestrear 1,565 hombres y 1,565 mujeres
#' para detectar 5% de diferencia con 80% poder.
#'
#' si pensamos que en _un grupo la proporción es de 10% y en el otro de 5%_:
pwr.2p.test(h=ES.h(p1 = 0.10, p2 = 0.05),
power = 0.8,
sig.level = 0.05,
alternative = "two.sided")
#' Necesitamos muestrear 424 hombres y 424 mujeres
#' para detectar 5% de diferencia con 80% poder.
# alternativa directa -----------------------------------------------------
power.prop.test(p1 = 0.55,
p2 = 0.50,
sig.level = 0.05,
power = .80)