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;

    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; 
xtrapoint=4; bmi=15.9375; vitd=39.2; output; 
run;

data total;
set irishwomen extra;

title 'no extra points';
proc glm data=total; where xtrapoint=0;
model vitd=bmi / solution clparm;
estimate 'level at bmi=25' intercept 1 bmi 25;
run;

title 'extra point: bmi=25, vitd=130';
proc glm data=total; where xtrapoint in (0,1);
model vitd=bmi / solution clparm;
estimate 'level at bmi=25' intercept 1 bmi 25;
run;

title 'extra point: bmi=50, vitd=80';
proc glm data=total; where xtrapoint in (0,2);
model vitd=bmi / solution clparm;
estimate 'level at bmi=25' intercept 1 bmi 25;
run;

title 'extra point: bmi=50, vitd=10';
proc glm data=total; where xtrapoint in (0,3);
model vitd=bmi / solution clparm;
estimate 'level at bmi=25' intercept 1 bmi 25;
run;

title 'with under weight woman: bmi=15.9, vitd=39.2';
proc glm data=total; where xtrapoint in (0,4);
model vitd=bmi / solution clparm;
estimate 'level at bmi=25' intercept 1 bmi 25;
run;