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

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

/*
Women from Ireland 
*/ 

DATA irishwomen;
	SET vitamind (WHERE = (country=4 and category=2));
    if bmi > 0 and bmi < 25 then bmigroup = 0;
    if bmi >= 25 and bmi < 30 then bmigroup = 1;
    if bmi >= 30 then bmigroup = 2;

logvitd=log10(vitd);    
RUN;

ods trace on;
ods output estimates=comparisons;

PROC GLM DATA=irishwomen;
    CLASS bmigroup;
    MODEL logvitd=bmigroup / solution clparm alpha=0.025;
    ESTIMATE 'over weight vs. normal weight' bmigroup -1 1 0;
    ESTIMATE 'obese vs. over weight' bmigroup 0 -1 1;
RUN;

ods trace off;

data comparisons;
set comparisons;

ratio_estimate=10**Estimate;
ci_lower=10**LowerCL;
ci_upper=10**UpperCL;
RUN;

proc print data=comparisons;
    var Parameter ratio_estimate ci_lower ci_upper;
run;