* PS11.2, using DATA10-5, for evaluating Model B - Table 11.2 
time 1960 1
sample 1960. 1994.
read(data10-5) year calwage uswage
genr t=time(0)
print t
genr tsq = t*t
sample 1960. 1989. 
dim yhat 50
auto calwage  t tsq / max drop
fc / max predict=yhat iblup beg=1960. end=1994.  
smpl 1990. 1994.
ols calwage yhat 
sample 1961. 1994.
genr yf = int(100*yhat+0.5)/100 
genr ape = int(100*((100*abs(calwage-yhat)/calwage)+0.5))/100
print year calwage yf ape
delete / all
stop