OPTIONS LS=75;

DATA a;

  INFILE 'coronary.dat';
  
  INPUT Smoke AgeL AgeU PY Deaths;
  
  Ln_PY = Log(PY); age = (ageL+ageU)/2 - 60; product = age*smoke;

PROC GENMOD; MODEL deaths=age smoke product/ link= log  dist=poisson
                                             offset = Ln_PY;
  
PROC SORT; BY DESCENDING AgeL;

PROC GENMOD ORDER=DATA;

  CLASS AgeL;
  
  MODEL Deaths = AgeL Smoke / link   = log
                              dist   = Poisson
                              offset = Ln_PY
                              noint;