* do-file for Supplementary Exercise 10.17
version 17 /* works also with versions 13-16 */

import delimited "R:\Chapter 10\ex10_017.csv"
scatter vo2 hr
twoway (scatter vo2 hr) (lfit vo2 hr)
regress vo2 hr
predict fit, xb
predict pe_fit, stdf
scalar tstar=invttail(17,.025)
generate pe_low=fit-tstar*pe_fit
generate pe_upp=fit+tstar*pe_fit
list, separator(5)