FILENAME goptions URL "http://192.38.117.59/~linearpredictors/datafiles/goptions.sas";
%include goptions;

FILENAME url URL "http://192.38.117.59/~linearpredictors/datafiles/readFibrosis.sas";
%include url;

ods output Estimates=or23vs01;

data fibrosis;
set fibrosis;    

log2ykl40=log2(ykl40);
log2p3np=log2(p3np);
log2ha=log2(ha);

fibrosis3=(stage>2);
fibrosis23=(stage>1);
fibrosis123=(stage>0);
RUN;

ods output Estimates=or23vs01;

title 'Estimates from table on the two last pages';
proc genmod data=fibrosis descending;
  model fibrosis23=log2ha log2p3np log2ykl40
              / dist=binomial link=logit;
estimate 'log2ha OR' log2ha 1 / exp;
estimate 'log2p3np OR' log2p3np 1 / exp;
estimate 'log2ykl40 OR' log2ykl40 1 / exp;
run;

DATA or23vs01;
SET or23vs01;
    
type=substr(Label,1,3);
RUN;

title 'Estimates with confidence intervals';
PROC PRINT NOOBS DATA=or23vs01; where type ne 'Exp';
    VAR Label LBetaEstimate LBetaLowerCL LBetaUpperCL;
RUN;

title 'ORs with confidence intervals';
PROC PRINT NOOBS DATA=or23vs01; where type='Exp';
    VAR Label LBetaEstimate LBetaLowerCL LBetaUpperCL;
RUN;