avallecam / powder

A Tidy Extension for Power Analysis
https://avallecam.github.io/powder/
Other
0 stars 0 forks source link

add power.prop.test #2

Open avallecam opened 3 years ago

avallecam commented 3 years ago
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)
avallecam commented 3 years ago

alternativa: habilitar argumentos para que usen la función ES.h desde el inicio

library(pwr)
pwr.p.test(h = ES.h(p1 = 0.28, p2 = 0.50), # effect size
           sig.level = 0.05,
           power = 0.80)

library(tidyverse)
library(powder)
expand_grid(p1=0.28,p2=seq(0.55,0.05,-0.1)) %>% 
  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.p.test)