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=ha p3np ykl40
              / dist=binomial link=logit;
estimate 'ha OR' ha 1 / exp;
estimate 'p3np OR' p3np 1 / exp;
estimate 'ykl40 OR' ykl40 1 / exp;
run;

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

title 'Estimates with SDs';
PROC PRINT NOOBS DATA=or23vs01; where type ne 'Exp';
    VAR Label LBetaEstimate StdErr;
RUN;

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