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 then bmigroup = 0;
    if bmi >= 25 then bmigroup = 1;

logvitd=log10(vitd);    
RUN;

PROC MEANS N MEAN MEDIAN DATA=irishwomen NWAY;
    CLASS bmigroup;
    VAR logvitd;
    OUTPUT OUT=summary N=count MEAN=mlogvitd MEDIAN=medlogvitd;
RUN;

DATA summary;
SET summary;

geom_average=10**mlogvitd;
back_median=10**medlogvitd;
RUN;

PROC PRINT DATA=summary;
RUN;

title 'for comparison';
PROC MEANS N MEAN MEDIAN DATA=irishwomen NWAY;
    CLASS bmigroup;
    VAR vitd;
RUN;