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

PROC GLM DATA=irishwomen;
    MODEL vitd=bmigroup;
    OUTPUT OUT=modelcontrol R=residual;
RUN;

PROC SORT DATA=modelcontrol; BY bmigroup;
RUN;

%begin_multilayout(2,2);

PROC UNIVARIATE DATA=modelcontrol GOUT=plots; BY bmigroup;
  VAR residual;  
  qqplot /  normal(mu=EST sigma=EST l=33);
RUN;

PROC UNIVARIATE DATA=modelcontrol GOUT=plots; 
  VAR residual;  
  qqplot /  normal(mu=EST sigma=EST l=33);
  HISTOGRAM / ENDPOINTS=-60 TO 60 BY 10 NORMAL;
RUN;
QUIT;

%end_multilayout(2,2);