avallecam / powder Goto Github PK

# 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)

library(pwr)

power.prop.test(p1 = 0.198, 
                p2 = 0.087, 
                sig.level = 0.05, 
                power = .80,
                alternative = "two.sided")

library(powder)

expand_grid(p1=0.198,p2=seq(0.18,0.08,-0.01)) %>% 
  mutate(h=purrr::pmap_dbl(.l = select(.,p1,p2),.f = ES.h)) %>% 
  mutate(power = 0.8,
         sig.level = 0.05,
         alternative = "two.sided") %>% 
  pwr_tidy(test_function = pwr.2p.test)