2015-05-16 54 views
0

尝试计算值对的比例。为什么使用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 

我做错了什么?

由于

+1

我不是比例测试专家,但你的问题困扰我。你从这个角度看到什么样的观点表明,添加数据不会改变数据的属性和后续的关键统计数据的计算? – EngBIRD

+0

注意输出的最后一行:'P值调整方法:holm'现在去做一些研究;) – Dason

回答

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")