Getting Started Example for PROC MULTTEST
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: MULTGS */
/* TITLE: Getting Started Example for PROC MULTTEST */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: multiple test, */
/* multiple comparisons */
/* PROCS: MULTTEST */
/* DATA: */
/* */
/* SUPPORT: Bob Derr */
/* REF: SAS/STAT User's Guide, PROC MULTTEST chapter */
/* MISC: */
/* */
/****************************************************************/
/*----------------------------------------------------------------
Made-up data to test effectiveness of a drug on 15 subjects.
3 balanced groups each receive 0, 1, and 2 mg of the drug, and
record presence or absence of 10 side effects.
Analyze each side effect individually, and adjust for multiplicity.
----------------------------------------------------------------*/
title 'Drug Example';
data Drug;
input Dose$ SideEff1-SideEff10;
datalines;
0MG 0 0 1 0 0 1 0 0 0 0
0MG 0 0 0 0 0 0 0 0 0 1
0MG 0 0 0 0 0 0 0 0 1 0
0MG 0 0 0 0 0 0 0 0 0 0
0MG 0 1 0 0 0 0 0 0 0 0
1MG 1 0 0 1 0 1 0 0 1 0
1MG 0 0 0 1 1 0 0 1 0 1
1MG 0 1 0 0 0 0 1 0 0 0
1MG 0 0 1 0 0 0 0 0 0 1
1MG 1 0 1 0 0 0 0 1 0 0
2MG 0 1 1 1 0 1 1 1 0 1
2MG 1 1 1 1 1 1 0 1 1 0
2MG 1 0 0 1 0 1 1 0 1 0
2MG 0 1 1 1 1 0 1 1 1 1
2MG 1 0 1 0 1 1 1 0 0 1
;
ods graphics on;
proc multtest bootstrap nsample=20000 seed=41287 notables
plots=PByTest(vref=0.05 0.1);
class Dose;
test ca(SideEff1-SideEff10);
contrast 'Trend' 0 1 2;
run;
ods graphics off;
