SSM Procedure

Trend Models for Irregular Data

A good reference for these models is De Jong and Mazzi (2001). Throughout this section h Subscript t Baseline equals left-parenthesis tau Subscript t plus 1 Baseline minus tau Subscript t Baseline right-parenthesis denotes the difference between the successive time points. The system matrices bold upper T Subscript t and bold upper Q Subscript t that govern these models depend on h Subscript t. However, whenever the notation is unambiguous, the subscript t is omitted.

Polynomial Spline Trend

This model is a general-purpose tool for extracting a smooth trend from the noisy data. The order of the spline governs the order of the local polynomial that defines the spline. The order-1 spline corresponds to Brownian motion (continuous-time random walk), the order-2 spline corresponds to integrated Brownian motion (continuous-time integrated random walk), and the order-3 spline provides a locally quadratic trend; the default order is 1. The dimension of the state underlying this component is the same as the order of the spline. The system matrices for the orders up to 3 are described as follows (in all the cases the initial condition is fully diffuse):

  • order-1 spline: bold upper Z equals left-parenthesis 1 right-parenthesis, bold upper T equals left-parenthesis 1 right-parenthesis, and bold upper Q equals sigma squared left-parenthesis h right-parenthesis

  • order-2 spline: bold upper Z equals left-parenthesis 1 0 right-parenthesis, bold upper T equals left-parenthesis 1 h comma 0 1 right-parenthesis, and bold upper Q equals sigma squared left-parenthesis StartFraction h cubed Over 3 EndFraction StartFraction h squared Over 2 EndFraction comma StartFraction h squared Over 2 EndFraction h right-parenthesis

  • order-3 spline: bold upper Z equals left-parenthesis 1 0 0 right-parenthesis, bold upper T equals left-parenthesis 1 h StartFraction h squared Over 2 EndFraction comma 0 1 h comma 0 0 1 right-parenthesis, and

    bold upper Q left-bracket i comma j right-bracket equals sigma squared asterisk StartFraction h Superscript 6 minus i minus j plus 1 Baseline Over left-parenthesis 6 minus i minus j plus 1 right-parenthesis left-parenthesis 3 minus i right-parenthesis factorial left-parenthesis 3 minus j right-parenthesis factorial EndFraction 1 less-than-or-equal-to i comma j less-than-or-equal-to 3

The system matrices for higher orders are similarly defined (for more information, see De Jong and Mazzi (2001)).

Note that, in addition to providing an estimate of the trend, this methodology can provide estimates of the higher-order derivatives of the trend. If alpha alpha denotes the k-dimensional subsection that is associated with a polynomial spline of order k, then its normal jth element (1 less-than-or-equal-to j less-than-or-equal-to k), alpha alpha left-bracket j right-bracket, corresponds to the derivative of order left-parenthesis j minus 1 right-parenthesis of this polynomial spline. For an example of the estimation of the first derivative of a trend component, see Example 33.12. For additional information about using these types of trend patterns in data analysis, see Eubank, Huang, and Wang (2003); Kohn, Ansley, and Tharm (1991); Selukar (2015).

Decay and Growth Trends

There are two choices for the decay trend: DECAY and DECAY(OU). Similarly, there are two choices for the growth trend: GROWTH and GROWTH(OU). The "OU" stands for the Ornstein-Uhlenbeck form of these models. The decay trend is a sum of two correlated components: one component is a random walk, and the other component is a stationary autoregression. In its Ornstein-Uhlenbeck form, the random walk component is replaced by a constant. The growth trend (and its Ornstein-Uhlenbeck variant) has the same form as the decay trend except that the autoregression is nonstationary (in fact, it is explosive). For growth trend models, floating-point errors can result for even moderately long forecast horizons because of the explosive growth in the trend values.

The system matrices for the decay and the growth types in their respective cases are identical, except for the sign of the rate parameter phi: phi less-than 0.0 for the decay type, and phi greater-than 0.0 for the growth type. In addition, the initial conditions for the growth models are fully diffuse; they are only partially diffuse for the decay models. The underlying state vector for all these models is two-dimensional.

The system matrices for the DECAY type are

StartLayout 1st Row 1st Column bold upper Z 2nd Column equals 3rd Column left-parenthesis 1 1 right-parenthesis 2nd Row 1st Column bold upper T 2nd Column equals 3rd Column normal upper D normal i normal a normal g left-parenthesis 1 comma exp left-parenthesis h phi right-parenthesis right-parenthesis 3rd Row 1st Column bold upper Q 2nd Column equals 3rd Column StartFraction sigma squared Over phi cubed EndFraction left-parenthesis h phi Baseline 1 minus exp left-parenthesis h phi right-parenthesis comma 1 minus exp left-parenthesis h phi right-parenthesis left-parenthesis exp left-parenthesis 2 h phi right-parenthesis minus 1 right-parenthesis slash 2 right-parenthesis EndLayout

The initial condition is partially diffuse with bold upper Q bold 1 equals normal upper D normal i normal a normal g left-parenthesis 0 comma StartFraction minus sigma squared Over 2 phi cubed EndFraction right-parenthesis. The system matrices for the GROWTH type are the same (with phi greater-than 0.0), except that the initial condition is fully diffuse; so bold upper Q bold 1 equals 0.

For the DECAY(OU) type, bold upper Z and bold upper T are the same as DECAY, whereas

bold upper Q equals normal upper D normal i normal a normal g left-parenthesis 0 comma sigma squared StartFraction left-parenthesis exp left-parenthesis 2 h phi right-parenthesis minus 1 right-parenthesis Over 2 phi EndFraction right-parenthesis normal a normal n normal d bold upper Q bold 1 equals normal upper D normal i normal a normal g left-parenthesis 0 comma StartFraction minus sigma squared Over 2 phi EndFraction right-parenthesis

The system matrices for the GROWTH(OU) type are the same (with phi greater-than 0.0), except that the initial condition is fully diffuse; so bold upper Q bold 1 equals 0.

Last updated: June 19, 2025