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));
        id=_N_;
RUN;

PROC GLM DATA=irishwomen;
    MODEL vitd=bmi;
    OUTPUT OUT=fittedline P=predicted;
RUN;

DATA fittedline;
SET fittedline;

type=0; OUTPUT;
type=id; OUTPUT;
vitd=predicted; type=id; OUTPUT;
RUN;

PROC GPLOT DATA = fittedline; 
    PLOT vitd*bmi=type /
        NOLEGEND HAXIS=AXIS1 VAXIS=AXIS2;
    AXIS1
        MINOR = NONE
        OFFSET = (5 pct)
        LABEL = ('BMI');
    AXIS2
        MINOR = NONE
        LABEL = (A=90 R=0 'Vitamin D');
    SYMBOL1 VALUE = CIRCLE I=RL COLOR = BLACK;
    SYMBOL2 VALUE = NONE I=JOIN L=33 COLOR = BLUE R=41;
RUN;
QUIT;