Размер выборки и расчет мощности в r как жизнеспособная альтернатива мощности в SAS?

Question

Поэтому я пытаюсь увидеть, насколько близки вычисления размера выборки (для двух независимых пропорций выборки с неравными размерами выборки) между мощностью процессора в SAS и некоторыми функциями размера выборки в r. Я использую данные, найденные здесь на сайте UCLA.

Сайт UCLA предоставляет параметры следующим образом:

p1 =.3, p2 =.15, мощность =.8, нулевая разница =0, и для двусторонних тестов предполагается, что выборки одинакового размера;

для тестов с неравным размером выборки параметры одинаковы: вес группы 1 для группы 1 и 2 для группы 2, а тесты, которые они выполняют, односторонние.

Я использую функцию r

pwr.t.test(n=NULL,d=0,sig.level=0.05,type="two.sample",alternative="two.sided")

от pwr пакет.

Поэтому, если я ввожу выбор параметров, как сайт UCLA для их первого примера, я получаю следующую ошибку:

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

Похоже, это связано с тем, что разница не обнаруживается с помощью r. Я установил d=.5, и он побежал. Будет ли SAS также давать ошибку из-за слишком маленькой разницы? В этом примере это не так, поскольку их нулевая разница также равна нулю.

Я также получаю ошибку выше при использовании

pwr.2p.test(h = 0, n = , sig.level =.05, power = .8) 

а также

pwr.chisq.test(w =0, N = , df =1 , sig.level =.05, power =.8 ).

Возможно, я делаю что-то ужасно неправильное, но, похоже, не могу найти способ, если предполагаемая разница равна 0.

Я понимаю, что SAS и r используют разные методы для расчета мощности, поэтому я не должен ожидать, что получу одинаковый результат. Я действительно просто пытаюсь понять, смогу ли я воспроизвести результаты силы власти в r.

Мне удалось получить почти идентичные результаты для первого примера с одинаковыми размерами выборки и двусторонней альтернативой, используя

 bsamsize(p1=.30,p2=.15,fraction=.5, alpha=.05, power=.8)

от Hmisc пакет. Но когда они делают односторонние тесты с неравным размером выборки, я не могу их повторить.

Есть ли способ повторить процесс в r для односторонних вычислений размера выборки для неравных размеров группы?

Приветствия.

Answer 1

В pwr.t.test и его производные, d это не нулевая разница (предполагается, что она равна нулю), а величина эффекта / предполагаемая разница между двумя популяциями. Если разница между средними значениями равна нулю, размер выборки не позволит обнаружить несуществующую разницу.

Если доля населения А составляет 15%, а доля населения Б - 30%, используйте функцию pwr::ES.h рассчитать величину эффекта и сделать тест пропорций, таких как:

> pwr.2p.test(h=ES.h(0.30,0.15),power=0.80,sig.level=0.05)

     Difference of proportion power calculation for binomial distribution (arcsine transformation) 

              h = 0.3638807
              n = 118.5547
      sig.level = 0.05
          power = 0.8
    alternative = two.sided

NOTE: same sample sizes

> pwr.chisq.test(w=ES.w1(0.3,0.15),df=1,sig.level=0.05,power=0.80)

     Chi squared power calculation 

              w = 0.2738613
              N = 104.6515
             df = 1
      sig.level = 0.05
          power = 0.8

NOTE: N is the number of observations