Documentation Example 3 for PROC STDRATE
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: SRATEX3 */
/* TITLE: Documentation Example 3 for PROC STDRATE */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: standardized rate */
/* PROCS: STDRATE */
/* DATA: */
/* */
/* SUPPORT: Yang Yuan UPDATE: April, 21, 2011 */
/* REF: PROC STDRATE, EXAMPLE 3 */
/* MISC: */
/****************************************************************/
data Factory;
input Age $ Event_E Count_E Event_NE Count_NE;
datalines;
20-29 31 352 143 2626
30-39 57 486 392 4124
40-49 62 538 459 4662
50-59 50 455 337 3622
60-69 38 322 199 2155
70+ 9 68 35 414
;
ods graphics on;
proc stdrate data=Factory
refdata=Factory
method=indirect(af)
stat=risk
plots(stratum=horizontal)
;
population event=Event_E total=Count_E;
reference event=Event_NE total=Count_NE;
strata Age / stats;
run;
ods graphics off;
data Factory1;
input Exposure $ Age $ Event Count;
datalines;
Yes 20-29 31 352
Yes 30-39 57 486
Yes 40-49 62 538
Yes 50-59 50 455
Yes 60-69 38 322
Yes 70+ 9 68
No 20-29 143 2626
No 30-39 392 4124
No 40-49 459 4662
No 50-59 337 3622
No 60-69 199 2155
No 70+ 35 414
;
proc stdrate data=Factory1
method=mh(af)
stat=risk
effect
;
population group(order=data exposed='Yes')=Exposure
event=Event total=Count;
strata Age;
run;