[英]What is the p.value for each individual test in my prop.test function?
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)))
除非您將conf.level
指定為0.95
以外的conf.level
,否則alpha級別默認為0.5:
prop.test(x, n, p = NULL,
alternative = c("two.sided", "less", "greater"),
conf.level = 0.95, correct = TRUE)
conf.level
與1-alpha
相同。 上面的代碼塊來自RStudio幫助菜單。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.