|
[multipart/alternative]
Dear SAS-l,
I am trying to use Bonferroni and Sidak adjusted p-Values using PROC MULTTEST. But something is wrong.
Thanks for any help
Carlos
[text/html]
/* BONFERRONI AND SIDAK ADJUSTED p-VALUES USING PROC MULTTEST*/
data ONE;
input TEST PVAL @@;
BON_ADJP = min(1,10*PVAL);
SID_ADJP = 1 - (1-PVAL)**10;
datalines;
1 0.0911 2 0.8912
3 0.0001 4 0.5718
5 0.0132 6 0.9011
7 0.2012 8 0.0289
9 0.0498 10 0.0058
;
data ONE; set ONE;
rename PVAL=RAW_P;
drop BON_ADJP SID_ADJP;
proc multtest data=ONE bon sid out=OUTP;
proc sort data=OUTP out=OUTP;
by RAW_P;
proc print data=OUTP;
run;
|
