Question

这个问题是部分方法论，部分编码

我的人口中，由于没有人开始疼痛，我发现13％的人出现疼痛（95％CI：9,18）。换句话说，预期比例= 0.13。

我想开发一个可以告诉我的表格或图表：

如果我进行了大小为N（N = 10,15,20等）的研究，为了发现与预期比例的显着差异，有多少比例的人需要发展疼痛（0.13）？ Alpha = 0.05，功率= 0.8。我不需要代码来制作表格/图表，但是如何找到这些数字！

显然，在较小的研究中（例如n = 10），人们需要更大的显着性差异。我很抱歉缺少代码，我试图在G * Power中解决这个问题，但我更愿意在R中这样做。谢谢!!

Answer 1

我会发布可能需要调整的旧代码（小n的问题，...）：

library(pwr)

estimatedP <- 0.13
estES <- pwr.p.test(n=150, sig.level=0.05, power=0.8, alternative="greater")$h

estES包含两个比例差异的效果大小。您的比例为estimatedP。现在有两种选择。您的estimatedP低于当时或高于第二组的比例。您正在寻找EShigh（比例高于13％）：

ESlow <- function(lowProp, estimP, estimES){
 # Estimates lower proportion
 2*asin(sqrt(estimP)) -2*asin(sqrt(lowProp))  - estimES
}

EShigh <- function(highProp, estimP, estimES){
 # Estimates higher proportion
 2*asin(sqrt(highProp)) -2*asin(sqrt(estimP))  - estimES
}

在下一步中，计算与效果大小相对应的比例：

highP <- uniroot(EShigh, c(0,1), tol=0.0001, estimES=estES, estimP=estimatedP)$root
lowP  <- uniroot(ESlow , c(0,1), tol=0.0001, estimES=estES, estimP=estimatedP)$root

看看比例：

highP # Proportion > 0.205 to be detectable
lowP  # Proportion < 0.07 to be detectable

您可以检查效果大小是否具有几乎相同的幅度但方向不同：

ES.h(lowP, estimatedP)
ES.h(estimatedP, highP)

R：样本大小和比例差异

1 个答案: