*
* Example 8/3 from pp 418-421
*
open data ukdeaths.dat
calendar(m) 1975
data(format=free,org=columns) 1975:1 1984:12 deaths seatbelt
*
spgraph
graph(footer="Figure 8-14 Monthly total deaths and serious injuries",grid=(t==1983:2))
# deaths
grtext(entry=1983:4,y=2200,align=left) "Seat belt\\legislation\\introduced"
spgraph(done)
*
boxjenk(regressors,ar=1,maxl) deaths
# constant seatbelt
*
@regcorrs
*
boxjenk(sdiffs=1,regressors,apply,maxl) deaths
# seatbelt
@regcorrs
*
boxjenk(sdiffs=1,ar=1,ma=1,sma=1,regressors,apply,maxl) deaths
# seatbelt