Programming the statistical procedures from SAS

Using Bonferroni to adjust Spearman Rank correlation p-values

Accepted Solution Solved
Reply
New Contributor
Posts: 2
Accepted Solution

Using Bonferroni to adjust Spearman Rank correlation p-values

I have run a correlation using the Spearman rank procedure, with the code below. I would like to use Bonferroni to correct the resulting p-values but am not sure how to integrate this into the code. Everything I have read applies Bonferroni to ANOVA.


proc corr data=three spearman;
var unanfl invas soprfl unprfl soprgr unprgr;
by trt;
run;

Accepted Solutions
Solution
‎07-26-2017 08:44 AM
Valued Guide
Valued Guide
Posts: 684

Re: Using Bonferroni to adjust Spearman Rank correlation p-values

You may want to seriously think of alternatives to Bonferroni if you have many correlations. This is because the Bonferroni method will result in very low power to detect nonzero correlations if you have many correlations. This approach actually sets the bound for the type I error rate, and is based on k independent tests. But the correlations are likely (highly) correlated with each other.

With a collection of k p-values (including from t tests of correlations), there are multiplicity adjustments that have a much higher power. These are not automatically done, but these are easily done using a data step or using PROC MULT with output from the CORR procedure. I suggest you read chapter 2 in the book Westfall, Tobias, Rom, Wolfinger, and Hochberg. 1999. Multiple Comparisons and Multiple Tests. (SAS Publishing). This shows how to do the analysis. You first need to decide on the adjustment and then follow the directions.

View solution in original post


All Replies
Trusted Advisor
Posts: 2,114

Re: Using Bonferroni to adjust Spearman Rank correlation p-values

Kristine,

The Bonferonni adjustment for multiple comparisons is easily done by the user (rather than the program). Just divide the alpha level by the number of tests under consideration and use that as the new alpha.

Since CORR does not put any *'s to indicate significance, that is really all there is to do.

Doc Muhlbaier
Duke
Solution
‎07-26-2017 08:44 AM
Valued Guide
Valued Guide
Posts: 684

Re: Using Bonferroni to adjust Spearman Rank correlation p-values

You may want to seriously think of alternatives to Bonferroni if you have many correlations. This is because the Bonferroni method will result in very low power to detect nonzero correlations if you have many correlations. This approach actually sets the bound for the type I error rate, and is based on k independent tests. But the correlations are likely (highly) correlated with each other.

With a collection of k p-values (including from t tests of correlations), there are multiplicity adjustments that have a much higher power. These are not automatically done, but these are easily done using a data step or using PROC MULT with output from the CORR procedure. I suggest you read chapter 2 in the book Westfall, Tobias, Rom, Wolfinger, and Hochberg. 1999. Multiple Comparisons and Multiple Tests. (SAS Publishing). This shows how to do the analysis. You first need to decide on the adjustment and then follow the directions.
🔒 This topic is solved and locked.

Need further help from the community? Please ask a new question.

Discussion stats
  • 2 replies
  • 1703 views
  • 0 likes
  • 3 in conversation