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));

    xtrapoint=0;
RUN;

data extra;    
xtrapoint=1; bmi=25; vitd=130; output;
xtrapoint=2; bmi=50; vitd=80; output; 
xtrapoint=3; bmi=50; vitd=10; output; 
run;

data total;
set irishwomen extra;

type='without extra point';
if xtrapoint>0 then type='with extra point';
RUN;

%begin_multilayout(2, 2);

PROC GPLOT DATA = total GOUT=plots; where xtrapoint=0;
    PLOT vitd*bmi /
        NOLEGEND HAXIS=AXIS1 VAXIS=AXIS2;
    AXIS1 ORDER=(15 TO 50 by 5) 
        MINOR = NONE        
        OFFSET = (5 pct)
        LABEL = ( 'BMI');
    AXIS2 ORDER=(0 to 140 BY 20)        
        MINOR = NONE
        LABEL = (A=90 R=0  'Vitamin D');
    SYMBOL1 VALUE = CIRCLE I=sm99s COLOR = black;
RUN;

PROC GPLOT DATA = total GOUT=plots; where xtrapoint in (0,1);
    PLOT vitd*bmi=type /
        NOLEGEND HAXIS=AXIS1 VAXIS=AXIS2;
    PLOT2 vitd*bmi /
        VAXIS=AXIS3;
    AXIS1 ORDER=(15 TO 50 by 5) 
        MINOR = NONE        
        OFFSET = (5 pct)
        LABEL = ( 'BMI');
    AXIS2 ORDER=(0 to 140 BY 20)
        MINOR = NONE
        LABEL = (A=90 R=0  'Vitamin D');
    AXIS3 ORDER=(0 to 140 BY 20)
        VALUE=NONE
        MINOR = NONE
        LABEL = (A=90 R=0  '');
    SYMBOL1 VALUE = STAR I=NONE COLOR = red;
    SYMBOL2 VALUE = CIRCLE I=sm99s COLOR = black;
    SYMBOL3 VALUE = NONE I=sm99s COLOR = red;
RUN;

PROC GPLOT DATA = total GOUT=plots; where xtrapoint in (0,3);
    PLOT vitd*bmi=type /
        NOLEGEND HAXIS=AXIS1 VAXIS=AXIS2;
    PLOT2 vitd*bmi /
        VAXIS=AXIS3;
    AXIS1 ORDER=(15 TO 50 by 5) 
        MINOR = NONE
        OFFSET = (5 pct)
        LABEL = ( 'BMI');
    AXIS2 ORDER=(0 to 140 BY 20)
        MINOR = NONE
        LABEL = (A=90 R=0  'Vitamin D');
    AXIS3 ORDER=(0 to 140 BY 20)
        VALUE=NONE
        MINOR = NONE
        LABEL = (A=90 R=0  '');
    SYMBOL1 VALUE = STAR I=NONE COLOR = red;
    SYMBOL2 VALUE = CIRCLE I=sm99s COLOR = black;
    SYMBOL3 VALUE = NONE I=sm99s COLOR = red;
RUN;

PROC GPLOT DATA = total GOUT=plots; where xtrapoint in (0,2);
    PLOT vitd*bmi=type /
        NOLEGEND HAXIS=AXIS1 VAXIS=AXIS2;
    PLOT2 vitd*bmi /
        VAXIS=AXIS3;
    AXIS1 ORDER=(15 TO 50 by 5) 
        MINOR = NONE    
        OFFSET = (5 pct)
        LABEL = ( 'BMI');
    AXIS2 ORDER=(0 to 140 BY 20)
        MINOR = NONE
        LABEL = (A=90 R=0  'Vitamin D');
    AXIS3 ORDER=(0 to 140 BY 20)
        VALUE=NONE
        MINOR = NONE
        LABEL = (A=90 R=0  '');
    SYMBOL1 VALUE = STAR I=NONE COLOR = red;
    SYMBOL2 VALUE = CIRCLE I=sm99s COLOR = black;
    SYMBOL3 VALUE = NONE I=sm99s COLOR = red;
RUN;

%end_multilayout(2,2);