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 bmi_score='21.75';
if bmi>=25 and bmi<30 then bmi_score='27.5';
if bmi>=30 and bmi<100 then bmi_score='32.5';
RUN;

PROC MEANS DATA=irishwomen;
    CLASS bmi_score;
    VAR vitd;
    OUTPUT OUT=averages MEAN=mvitd;
RUN;

PROC GPLOT DATA = averages; 
    PLOT mvitd*bmi_score /
        HAXIS=AXIS1 VAXIS=AXIS2;
    AXIS1
        MINOR = NONE
        OFFSET = (5 pct)
        LABEL = ('BMI score');
    AXIS2
        MINOR = NONE
        LABEL = (A=90 R=0 'Vitamin D');
    SYMBOL1 VALUE = CIRCLE I=NONE COLOR = BLACK;
RUN;
QUIT;