Programming the statistical procedures from SAS

Using Bonferroni to adjust Spearman Rank correlation p-values

Reply
New Contributor
Posts: 2

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;
Valued Guide
Posts: 2,108

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
Valued Guide
Valued Guide
Posts: 679

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.
Ask a Question
Discussion stats
  • 2 replies
  • 1413 views
  • 0 likes
  • 3 in conversation