Hello,
I looked at page 9
https://support.sas.com/resources/papers/proceedings09/240-2009.pdf
SAS Global Forum 2009 -- Paper 240-2009 Beyond SAS/GENETICS™ Amina Barhdadi, Yassamin Feroz zada, Marie-Pierre Dubé Montreal Heart Institute and Université de Montréal, Montreal, Canada
The code is clear.
... ...
If cumulative>0.995 then Meff=number; /* cumulative and number are two provided variables in eigenvalues data set*/
Run;
/*Step4-Apply the Bonferroni correction formula to adjust point-wise significance level*/
Data Meff;
Set Meff;
alphaG=alphaE/Meff; /*alphaE is the experiment-wise error rate*/
Run;
Meff is an estimate of the effective number of independent tests. In your use case Meff = 7 , so you divide the test-wise error rate by 7 to get the experiment-wise error rate (Bonferroni correction). Like 0.05 / 7 = 0,007143 .
Important Note: The last comment line says : /*alphaE is the experiment-wise error rate*/
The experiment-wise error rate is on the left-hand side of the equality sign of course.
What you divide by 7 (Meff) is the test-wise error rate. So, the comment is wrong (should say alphaG) ... or alphaG and alphaE should switch places.
Koen
... View more