两个样本z检验的功效分析

Question

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

所以我用以下代码进行了z分数测试：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包，请解释它是如何工作的。

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

Answer 1

我认为你在使用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

  

＆＃34;效果大小将是平均值的平均值   偏差。效果尺寸越大，a的功率越大   给定样本大小。或者，效果大小越大，样本越小   达到相同功率所需的尺寸。所以，很好的估计效果   尺寸是良好功率分析的关键。但它并不总是一个   轻松确定效果大小的任务。对效果大小的良好估计   来自现有文献或试点研究＆＃34;