在R中使用prop.test的

Question

我需要一些澄清有關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) 

最後一行上面的例子返回一個0.09914的p值。

然而，當我們直接輸入上述數值，我們得到了一個完全不同的p值：

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

上面的代碼行返回0.8382的p值。

爲什麼會發生這種情況？

Answer 1

不要除以組中的數字。這將產生嚴重影響的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