r - 两个样本(不等方差)welch t 检验的统计功效和样本量确定？

Question

两个样本 t 检验(等方差合并标准偏差)的命令是

power.t.test(n=, delta=, sd=, type="two.sample")

给定不等方差和样本数的两个样本，我如何计算统计功效？

Answer 1

您可以使用 MESS 包中的 power_t_test() 函数。

例如，要计算个体比例为 1:2 时所需的样本量，sd 为 1 和 3，效应量为 1.2(对于 80% 的功效)

> MESS::power_t_test(n=NULL, sd=1, power=.8, ratio=2, sd.ratio=3, delta=1.2)

     Two-sample t test power calculation with unequal sizes 

              n = 14.35824, 28.71648
          delta = 1.2
             sd = 3
      sig.level = 0.05
          power = 0.8
    alternative = two.sided

NOTE: n is vector of number in each group

要计算功率，您可以设置 power=NULL 并指定 n 而不是常规的 power.t.test()函数。

> MESS::power_t_test(n=20, sd=1, power=NULL, ratio=2, sd.ratio=3, delta=1.2)

     Two-sample t test power calculation with unequal sizes 

              n = 20, 40
          delta = 1.2
             sd = 3
      sig.level = 0.05
          power = 0.9149541
    alternative = two.sided

NOTE: n is vector of number in each group