*
* Example 1.3.6 from pp 13-14
*
open data deaths.dat
calendar(m) 1973
data(format=free,org=columns) 1973:1 1978:12 deaths
source bdtrends.src
*
* The procedure HarmonicFit performs the harmonic regression using the listed set
* of periodicities, producing as output the fitted values.
*
@harmonicfit(periods=||12,6||) deaths / hdeaths
*
* The graph in the text uses two different styles to show the actual data
* (squares) and fitted components (lines). The default behavior of GRAPH is to
* use distinct line patterns or colors instead. To mimic the appearance in the
* text, you can use the pair of option OVERLAY=DOTS and OVSAME. Overlay graphs
* do a pair of graphs in the same box. OVSAME requests that they use the same
* vertical scale, which is what we want here.
*
graph(footer="Figure 1-11 Estimated Harmonic Component of Accidental Deaths",overlay=dots,ovsame) 2
# hdeaths
# deaths