/******************************************************************
* Multiplicity: adjust one prespecified family of four raw p-values
* Holm controls FWER without an independence assumption.
******************************************************************/

proc import datafile="example.csv" out=pvalues dbms=csv replace;
  guessingrows=max;
run;

proc sort data=pvalues;
  by ordern;
run;

proc multtest inpvalues(raw_p)=pvalues holm bonferroni out=adjusted;
run;

proc print data=adjusted noobs;
  var endpoint family ordern raw_p stpbon_p bon_p;
run;
