使用Bonferroni校正从R上的stats包中的“pairwise.prop.test”中的P值

时间:2017-06-08 19:10:25

标签: r

第一次问一个问题。我已经尽力进行自己的搜索并查看了手册。我希望我不会发布重复或偏离主题。

我使用pairwise.prop.test和bonferroni校正输出:

pairwise.prop.test(fishingsept13,p.adjust.method = 'bonferroni')

它出来很好,运行顺利。

我的问题更多的是对印刷的p值的解释。我知道我的数据的Bonferroni校正将是任何小于0.01的重要事项。这些是我的结果:

我认为我会用我已知的Bonferonni解释它们并从那里确定它们的意义。我是偏执,我假设错了。

1 个答案:

答案 0 :(得分:1)

实际上,R中的成对比较函数已经为你纠正了p值。您使用的p.adjust.method参数告诉R通过Bonferroni方法执行多重比较校正。由于您执行了3次单独的比较,因此您看到的p值已经乘以3:您不必调整显着性阈值。要亲眼看看,看看运行时会发生什么

pairwise.prop.test(fishingsept13,p.adjust.method = 'none')