log using "fig1.log", replace
*
*  fig1.do
*
*      N.B.  The gprefs statements given below will 
*            overide your Custom 1 color settings!
*
*  This is a Stata program that produces figure 1 from
*  the manuscript.
*
use "http://www.mc.vanderbilt.edu/prevmed/wddtext/data/2.20.Framingham.dta", clear
*
*  Drop ten extreme outliers.
*
drop if bmi>55
gprefs set window scheme custom1
gprefs set custom1 background_color 255 255 255
gprefs set custom1 pen1_color   0   0   0
gprefs set custom1 pen2_color   0   0 255
*
*  Plot conventional scatter plot with overstrike.
*
graph dbp bmi,symbol(o) xlabel(20 25 to 50) ylabel(50 70 to 150) gap(3)
log off