Cat
Cat

Reputation: 1

P-values in a "pairwise.prop.test" from the stats package on R using a Bonferroni correction

First time asking a question. I have done my best to conduct my own searches and have looked at the manual. I hope I am not posting a repeat or an off-topic.

I have an output from using the pairwise.prop.test with a bonferroni correction using:

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

It came out just fine and ran smoothly.

My question is more on the interpretation of the p-values printed. I know that the Bonferroni correction for my data would be that anything less than 0.01 is significant. These are my results:

I assume that I would interpret them with my known Bonferonni and determine their significance from there. I'm paranoid that I am assuming wrong.

Upvotes: 0

Views: 1223

Answers (1)

Matt
Matt

Reputation: 994

Actually, the pairwise comparison functions in R already correct the p values for you. The p.adjust.method argument you used told R to perform a multiple comparisons correction by the Bonferroni method. Since you performed 3 separate comparisons, the p values you see have already been multiplied by 3: you don't have to adjust your significance threshold. To see for yourself, look what happens when you run

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

Upvotes: 1

Related Questions