ODS RTF FILE = 'C:\Documents and Settings\rjchar2\My Documents\CPH931F09\WA3SAS931F09.rtf';

PROC IMPORT DATAFILE = 'C:\Documents and Settings\rjchar2\My Documents\CPH931F09\PostOpExcel.xls'
     OUT = PostOp  DBMS = EXCEL REPLACE;
     SHEET = Sheet1;
     GETNAMES = YES;
RUN;

PROC SORT DATA=PostOp OUT=PostOpSorted;
BY DESCENDING STATUS;
RUN;

PROC CATMOD ORDER=DATA DATA=PostOpSorted;
DIRECT   HIGHTEMPI HIGHTEMPS HIGHBP COMFORT;
MODEL    STATUS = HIGHTEMPI HIGHTEMPS HIGHBP COMFORT;
RUN; QUIT;


* Alternatively, one may use the following code, which I discovered
* only recently.  The results are the same except that now  -2 LOG L
* is included in the output.
*
* PROC LOGISTIC DATA=PostOpSorted;
* MODEL    STATUS(REFERENCE=FIRST) = HIGHTEMPI HIGHTEMPS HIGHBP COMFORT / LINK=GLOGIT;
* RUN; QUIT;


PROC LOGISTIC DATA=PostOpSorted DESCENDING;
MODEL STATUS = HIGHTEMPI HIGHTEMPS HIGHBP COMFORT;
RUN;


PROC IMPORT DATAFILE = 'C:\Documents and Settings\rjchar2\My Documents\CPH931F09\SARS.xls'
     OUT = GenLinMod  DBMS = EXCEL REPLACE;
     SHEET = SARS;
     GETNAMES = YES;
RUN;

PROC GENMOD DATA=GenLinMod;
WHERE COUNTRY='SINGAPORE';
MODEL DAILYINF = TIME  / LINK=LOG DIST=POISSON;
RUN;

PROC GPLOT DATA=GenLinMod; 
WHERE COUNTRY='SINGAPORE';
PLOT DAILYINF*TIME; 
RUN; 


PROC GENMOD DATA=GenLinMod;
WHERE COUNTRY='SINGAPORE';
MODEL DAILYINF = TIME TIME2 / LINK=LOG DIST=POISSON;
RUN;



PROC GENMOD DATA=GenLinMod;
WHERE COUNTRY='SINGAPORE';
MODEL DAILYINF = TIME TIME2 / LINK=LOG DIST=POISSON SCALE=DEVIANCE;
RUN;


PROC GENMOD DATA=GenLinMod;
WHERE COUNTRY='SINGAPORE';
MODEL DAILYINF = TIME TIME2 / LINK=LOG DIST=NEGBINOM;
RUN;

ODS RTF CLOSE;
RUN;