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;

DATA irlwomen;
	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;
RUN;

/*********************************
    Adding groupwise means to the dataset
    */

PROC MEANS mean;
    CLASS bmigroup;
    VAR vitd;
    OUTPUT OUT=means mean=mean;
RUN;

DATA irlwomen;
    SET irlwomen means;
RUN;

/*********************************
    Mean level of Vitamin D plotted against BMI-groups 
    */

PROC GPLOT DATA = irlwomen;
    PLOT vitd*bmigroup mean*bmigroup/ OVERLAY
        HAXIS=AXIS1 VAXIS=AXIS2;
    AXIS1
        MINOR = NONE
        OFFSET = (20 pct)
        ORDER = (0 to 2 by 1)
        LABEL = ('BMI group');
    AXIS2
        MINOR = NONE
        LABEL = (A=90 R=0 'Vitamin D');
    SYMBOL1 VALUE = CIRCLE COLOR = BLACK;
    SYMBOL2 VALUE = DOT COLOR = RED;
RUN;
QUIT;