R和单尾t检验的功率

时间:2017-05-10 20:31:35

标签: r

我将进行单尾t检验,以确定一个均值是否明显低于另一个。问题是,当我使用 R pwr 包来确定 n = 30时我能期望的功率时,我得到了一个极端的即使对于大型效果也能降低功率所以,例如:

> pwr.t.test(d=0.8,sig.level=.05,n=30,alternative="less")

  Two-sample t test power calculation 

             n = 30
             d = 0.8
     sig.level = 0.05
         power = 1.251823e-06
   alternative = less

  NOTE: n is number in *each* group

更奇怪的是,当我增加 n 时,电源会下降。因此,例如,将 n 提升到300就可以了:

> pwr.t.test(d=0.8,sig.level=.05,n=300,alternative="less")

  Two-sample t test power calculation 

             n = 300
             d = 0.8
     sig.level = 0.05
         power = 0
   alternative = less

  NOTE: n is number in *each* group

我错过了什么?

1 个答案:

答案 0 :(得分:3)

我想这是因为d和替代='较少'是在不同的方向'。 试试这个,你会明白我的意思。

pwr.t.test(d= - 0.8,sig.level=.05,n=300,alternative="less")

 Two-sample t test power calculation 

          n = 300
          d = -0.8
  sig.level = 0.05
      power = 1
alternative = less

NOTE: n is number in *each* group