*
* Example 1.5.4 from page 32
*
open data deaths.dat
calendar(m) 1973
data(format=free,org=columns) 1973:1 1978:12 deaths
*
* @ClassicalDecomp has options SPAN for the seasonal span (here not necessary,
* since the data have already been declared to be monthly),
* TREND=NONE/LINEAR/QUADRATIC to choose the type of trend, and FACTORS to get as
* a return the seasonal component.
*
@ClassicalDecomp(trend=quadratic,factors=s) deaths
set deseas = deaths-s
graph(footer="Figure 1-24 Deseasonalized Accidental Deaths")
# deseas
graph(footer="Figure 1-25 Seasonal Component")
# s