*
* Example SERIESF.DAT from pp 187-189
*
open data seriesf.dat
calendar(m) 1995
data(format=free,org=columns) 1995:1 1998:12 demand
*
* This time, we use the ClassicalDecomp procedure
*
@ClassicalDecomp(trend=linear,seasonal=mult,tc=tc,factors=factors) demand
set fitted = tc*factors
set error  = demand-fitted
@uforeerrors demand fitted
*
* %%ferrmse is produced by @uforeerrors. To match the "z" scores in the text, we
* need to multiply the standard error by sqrt(48.0/46.0) so we're using the same
* divisor.
*
set z = error/%%ferrmse*sqrt(46.0/48.0)
*
print(picture="*.###") / demand factors tc error z