1

我需要计算对数据执行的统计检验的功效。我有 2 个 csv 文件，每个文件的样本大小为 50。均值的差异在 0.05 水平上具有统计显着性。样本来自方差未知的正态分布。

所以我使用以下代码执行了 az score 测试：X 和 Y 是两个样本，每个样本大小为 50。

zTest <- function(x, y) {
Difference <- (mean(x) - mean(y)) # difference between the two sample means
seDifference <- sqrt(((sd(x)^2)/length(x)) + ((sd(y)^2)/length(y))) #standard error for difference
zScore <- Difference/seDifference # z score
return(zScore) # return z score
}

我得到的 Z 分值为 -15.78006

现在我需要计算上面执行的统计测试的功效。我的问题是如何从这里找到力量。公式是什么？以及如何在 R 中应用它。如果您建议在 R 中使用 pwr 包，请解释它是如何工作的。

如果我含糊不清，请提前致谢并道歉。我是功率分析的新手。

4

1 回答 1

2

我认为您在使用 pwr 包时需要效果大小。

pwr.t.test(n = 50, d = NULL, sig.level = 0.05, power = NULL,
           type = c("two.sample"),
           alternative = c("two.sided"))

d = NULL你错过了什么

关于效果大小的一些解释：http: //www.ats.ucla.edu/stat/r/dae/t_test_power2.htm

“效应量将是均值与合并标准差的差异。效应量越大，给定样本量的功效就越大。或者，效应量越大，达到相同功效所需的样本量就越小。因此，对效果大小的良好估计是良好功效分析的关键。但确定效果大小并不总是一件容易的事。对效果大小的良好估计来自现有文献或试点研究"

于 2016-04-15T19:58:29.093 回答