r - r 中的样本大小和功率计算可以作为 SAS 中 proc power 的可行替代方案吗？

Question

因此，我试图了解 SAS 中的 proc power 与 r 中的某些样本大小函数之间的样本大小计算(对于样本大小不等的两个样本独立比例)有多接近。我正在使用此处找到的数据 at a UCLA website.

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 也会给出错误吗？在示例中并非如此，因为它们的零差值也为零。

我在使用时也遇到了上面的错误

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 中复制 proc power 结果。

对于第一个示例，我能够在相同的样本大小和使用双面替代方案的情况下获得几乎相同的结果

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

来自Hmisc包。但当他们用不等的样本量进行单方面测试时，我无法复制这些。

有没有办法在 r 中复制该过程，以计算不相等的组大小的单边样本大小？

干杯。

Answer 1

在pwr.t.test及其衍生产品中，d不是零差异(假设为零)，而是效果大小/假设差异两个人口。如果总体均值之间的差异为零，则任何样本量都无法让您检测到不存在的差异。

如果人群 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