*
* Example 1.3.5 from pp 10-12
*
open data lake.dat
calendar 1875
data(format=free,org=columns) 1875:1 1972:1 lake
*
set trend = t
linreg lake
# constant trend
prj fitted
*
graph(footer="Figure 1-9 Level of Lake Huron Showing Least Squares Fit") 2
# lake
# fitted
*
* The series %RESIDS has the residuals from the most recent regression
*
set y = %resids
graph(footer="Figure 1-10 Residuals from Fitting a Line to Lake Data")
# y