Resources

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;

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;