我的prop.test函数中每个单独测试的p.value是多少?

问题描述 投票:0回答:1

prop.test函数中每个单独测试的p值是多少? (见下面的代码)。

当进行多次测试(在这种情况下k = 10000次测试)时,我想在prop.test中找到每个单独测试的alpha,因为它显然不是0.05。函数prop.test必须以某种方式调整p值。

在互联网和R中阅读(在控制台中键入prop.test)时,我没有找到p值是什么以及如何获得它的答案。我意识到在控制台中键入prop.test可以提供帮助,但我尝试了它并且无法理解prop.test-code,足以理解alpha(对于单个测试)是什么。因此,如果有人可以解释我如何获得alpha(个别测试)的数字,我将不胜感激。

非常感谢,/ Pedram

k <- 10000

pH0 <- 1:k
pH1 <- 1:k
nA <- 4000
nB <- 4000
p0 <- 0.01
p1 <- 2*p0

a = 0.05
pD = a

for(i in 1:k)
{
      x1 <- rbinom(1,nA,p0)
      x2 <- rbinom(1,nB,p0)
      y <- rbinom(1,nA,p1)

      pH0[i] <- prop.test(c(x1,x2),c(nA,nB))$p.value < pD
      pH1[i] <-prop.test(c(x1,y),c(nA,nA))$p.value < pD 
}
print(paste("FDR =",sum(pH0)/sum(pH0+pH1)))
print(paste("True positives =",sum(pH1)/k))
print(paste("FWER =",(1-((1-pD)^k)))
r bioinformatics p-value
1个回答
0
投票

除非您将conf.level指定为0.95以外的其他值,否则alpha级别默认为0.5:

prop.test(x, n, p = NULL,
          alternative = c("two.sided", "less", "greater"),
          conf.level = 0.95, correct = TRUE)

conf.level1-alpha相同。上面的代码块来自RStudio帮助菜单。

© www.soinside.com 2019 - 2024. All rights reserved.