两个样本（不等方差）韦尔奇 t 检验的统计功效和样本量确定？答案

【问题标题】：Statistical power and sample size determination for two sample (unequal variance) welch t-test?两个样本（不等方差）韦尔奇 t 检验的统计功效和样本量确定？
【发布时间】：2018-07-17 16:06:07
【问题描述】：

两个样本t检验（等方差合并标准差）的命令是

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

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

【问题讨论】：

您是在尝试确定达到特定功效水平所需的样本数量，还是说您知道样本量并想要计算检验的功效？

标签： r

【解决方案1】：

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

例如，当您的个体比例为 1:2 时，计算所需的样本量，标准差为 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

【讨论】：

是的，这解决了我的问题。令人惊讶的是，Rstuio 帮助中没有提到“ratio”和“sd.ratio”
它们不是基础 R 附带的 power.t.test 函数的一部分。您需要一个类似于 power_t_test() 的扩展函数
我将对两个样本（不等方差）使用这些命令 pwr.t2n.test(n1=NULL, n2= NULL, d = NULL, sig.level = 0.05, power = NULL,alternative= c("two.lateral", "less","greater")) 结合 ES.t.two(m1 = NULL, m2 = NULL, sd1 = NULL, sd2 = NULL, n1 = NULL, n2 = NULL, t = NULL, se1 = NULL, se2 = NULL, df = NULL, alternative = c("two.lateral", "one.lateral")) 估计cohen d。