我需要一些关于在R中使用prop.test命令的说明。
请看下面的例子:
pill <- matrix(c(122,478,99,301), nrow=2, byrow=TRUE)
dimnames(pill) <- list(c("Pill", "Placebo"), c("Positive", "Negative"))
pill
Positive Negative
Pill 122 478
Placebo 99 301
prop.test(pill, correct=F)
上例中的最后一行代码返回p值0.09914。
然而,当我们直接输入上述值时,我们会得到一个完全不同的p值:
prop.test(x=c(122/600,99/400), n=c(600,400), correct=F)
上面的代码行返回的p值为0.8382。
为什么会这样呢?
1条答案
按热度按时间vuktfyat1#
不要除以组中的数字。这会导致样本量大幅减少,从而严重影响p值。
您应该已经注意到您的呼叫的估计比例的奇怪结果: