R中属性检验的应用

sulc1iza  于 2023-04-03  发布在  其他
关注(0)|答案(1)|浏览(120)

我需要一些关于在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。
为什么会这样呢?

vuktfyat

vuktfyat1#

不要除以组中的数字。这会导致样本量大幅减少,从而严重影响p值。

prop.test(x=c(122,99), n=c(600,400), correct=F)

    2-sample test for equality of proportions without continuity
    correction

data:  c(122, 99) out of c(600, 400)
X-squared = 2.7194, df = 1, p-value = 0.09914
alternative hypothesis: two.sided
95 percent confidence interval:
 -0.097324375  0.008991042
sample estimates:
   prop 1    prop 2 
0.2033333 0.2475000

您应该已经注意到您的呼叫的估计比例的奇怪结果:

prop 1       prop 2 
0.0003388889 0.0006187500

相关问题