7

したがって、サンプルサイズの計算（サンプルサイズが等しくない2つのサンプルに依存しない比率の場合）が、SASのprocpowerとrのいくつかのサンプルサイズ関数の間でどれほど近いかを確認しようとしています。私はUCLAのウェブサイトでここにあるデータを使用しています。

UCLAサイトは、次のようにパラメーターを提供します。

p1 = .3、p2 = .15、power = .8、null Difference = 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は、差が小さすぎる場合にもエラーを出しますか？nullの差もゼロであるため、この例には含まれていません。

使用すると上記のエラーも発生します

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でprocpowerの結果を複製できるかどうかを確認しようとしています。

最初の例では、サンプルサイズが等しく、両側の代替案を使用して、ほぼ同じ結果を得ることができました。 

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

Hmiscパッケージから。しかし、サンプルサイズが等しくない片側テストを行う場合、それらを再現することはできません。

等しくないグループサイズの片側サンプルサイズ計算のためにrのプロセスを複製する方法はありますか？

乾杯。

4

1 に答える 1

3

とそのpwr.t.test導関数dは、ヌルの差（ゼロと見なされる）ではなく、2つの母集団間の効果量/仮説の差です。母平均間の差がゼロの場合、存在しない差を検出できるサンプルサイズはありません。

母集団Aの比率が15％で、母集団Bの比率が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
于 2013-03-14T14:14:37.473 に答える