我刚刚为R＆＃39 {s} power.t.test()遇到了一些奇怪的行为，如果我尝试运行它，它会返回负delta（＆＃34;真正意义上的差异＆＃34;）值对于power的低值。这是一个MWE：
power.t.test(n = 30,sd = 30,type = "paired",power = 0.01)$delta
[1] -2.073823
当我开始逐渐增加power时，结果delta会单调增加 - 如果所有值都是正数，这就是我所期望的，但是如果你考虑负增量值，这是非常奇怪的 - 对于type = "paired"和alternative = "two.sided"（标准）我认为delta应该作为绝对值返回。下面还有一些例子：
sapply(X = (1:10)/100,
FUN = function(x,...){power.t.test(power=x,...)$delta},
n = 30, sd = 30, type = "paired",
simplify = TRUE)
[1] -2.0738235 -0.5309280 0.4482013 1.1848257 1.7840618
[6] 2.2941601 2.7414432 3.1419327 3.5061993 3.8415337
power.t.test(n = 30,sd = 30,type = "paired",power = 0.02500001)$delta
[1] -4.235165e-22
我的猜测是power.t.test()使用的近似公式因power的低值而失败。对此问题的任何想法将不胜感激。
干杯， 菲利普