/*******************/
/* British doctors */
/*******************/

DATA Britdoc;
  SET Britdoc;
  Off = log(PersonYears/1000);
RUN;

PROC GENMOD DATA=Britdoc;
  CLASS Age Smoking;
  MODEL Deaths = Age Smoking / DIST=POI OFFSET=Off;
  ESTIMATE '35-44 Smoker' Intercept 1 Age 1 0 0 0 0 Smoking 0 1;
  ESTIMATE '55-64 Nonsmoker' Intercept 1 Age 0 0 1 0 0 Smoking 1 0;
RUN;

/*********************************/
/* AIDS Data with overdispersion */
/*********************************/
 
** Poisson;
PROC GENMOD DATA=aids;
  MODEL Cases = Year / DIST=POI;
  EFFECTPLOT FIT;
RUN;

** Quasi-poisson;
PROC GENMOD DATA=aids;
  MODEL Cases = Year / DIST=POI PSCALE;
  EFFECTPLOT FIT;
RUN;

** Negative binomial;
PROC GENMOD DATA=aids;
  MODEL Cases = Year / DIST=NB;
  EFFECTPLOT FIT;
RUN;