Небольшая ошибка размера выборки при вычислении мощности для двух пропорций

Question

Я пытаюсь использовать анализ мощности для выбора размера выборки с помощью библиотеки pwr в R.

      library(pwr)
pwr.2p2n.test(h = 0.1, n1 = 78, power = 0.8, sig.level = 0.0125)

Я получаю следующую ошибку:

      Error in uniroot(function(n2) eval(p.body) - power, c(2 + 1e-10, 1e+09)) : f() values at end points not of opposite sign

Если я изменю размер выборки (больше) или размер эффекта на средний (0,4), он будет работать. Любые решения были бы замечательными, спасибо.

Answer 1

Я не думаю, что можно достичь мощности 0,8 с любым размером выборки во второй группе (независимо от того, насколько большой), с другими вашими ограничениями, например

      library(pwr)
pfun <- function(n2) {
    pwr.2p2n.test(h = 0.1, n1 = 78, n2 = n2, sig.level = 0.0125)$power
}
pfun2 <- Vectorize(pfun)

png("ppow.png")
curve(pfun2(x), from = 100, to = 1e9, log="x")
dev.off()

Я пробовал "только" размер выборки до 10^9, но кажется, что независимо от того, насколько большой вы n2, вы не можете получить мощность> 0,054 или около того с этой настройкой. Так что это не вычислительная проблема, это структурная / статистическая проблема.