sysuse lifeexp, clear 
keep if region==3
gen lngnp = ln(gnppc)
quietly regress le lngnp 
predict r, resid 
twoway dropline r gnp, ///
  yline(0, lstyle(foreground)) mlabel(country) mlabpos(9) ///
  ylab(-6(1)6) ///
  subtitle("Regression of life expectancy on ln(gnp)" "Residuals:" " ", pos(11)) ///
  note("Residuals in years; positive values indicate" "longer than predicted life expectancy")