/*-----------------------------------------------------------------
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;