比例测试

时间:2015-05-16 17:05:16

标签: r

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

答案 0 :(得分:1)

删除&#34; Holm&#34;来自pairwise.prop.test的方法将返回相同的结果,这会调整p值 - 类似于&#34; Bonferroni&#34;。注意,后者舍入到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")