測試比例

Question

嘗試計算值對的比例。爲什麼使用prop.test一對或幾個（pairwise.prop.test）我得到了不同的結果？

如果我只對第一對我得到運行：P值= 0.08181

prop.test(c(73,68),c(86,93),alternative ="two.sided") 
## 
## 2-sample test for equality of proportions with continuity correction 
## 
## data: c(73, 68) out of c(86, 93) 
## X-squared = 3.0286, df = 1, p-value = 0.08181 
## alternative hypothesis: two.sided 
## 95 percent confidence interval: 
## -0.01122346 0.24653229 
## sample estimates: 
## prop 1 prop 2 
## 0.8488372 0.7311828 

當我超過1對運行我得到差異。 P值

smokers <- c(73, 68, 98, 70) 
patients <- c(86, 93, 136, 182) 
pairwise.prop.test(smokers, patients) 
## 
## Pairwise comparisons using Pairwise comparison of proportions 
## 
## data: smokers out of patients 
## 
## 1  2  3 
## 2 0.16 -  - 
## 3 0.12 0.98 - 
## 4 1.8e-11 4.4e-07 2.9e-08 
## 
## P value adjustment method: holm 

我做錯了什麼？

由於

Answer 1

卸下從pairwise.prop.test的「Holm的」方法將返回相同的結果，其調整的p值 - 類似於「邦費羅尼」。請注意，後者輪到3位數字。

prop.test(c(73,68),c(86,93), alternative ="two.sided") 
smokers <- c(73, 68, 98, 70) 
patients <- c(86, 93, 136, 182) 
pairwise.prop.test(smokers, patients, p.adjust.method = "none")