cap log close
clear
log using StataLog, replace 

include "http://192.38.117.59/~linearpredictors/datafiles/readVitaminD.do"
use "`irlwomen'"
  
regress vitd bmi
predict residuals, residuals
predict yhat, xb

/*******************************
Scatter plot of Residuals vs BMI         
********************************/

scatter residuals bmi, name(residualsBmi)

/*******************************
Scatter plot of Residuals vs
predicted values of Vitamin D         
********************************/

scatter residuals yhat, name(residualsVitd)

/*******************************
QQplot of residuals
********************************/

qnorm residuals, name(QQplot)

/*******************************
Histogram of residuals with overlayed
normal density
********************************/

histogram residuals, /*
*/ start(-30) width(15) xlabel(-22.5 (15) 52.5) normal /*
*/ percent name(histogram)

log close