SSM Procedure

Filtering Pass

The filtering pass sequentially computes the quantities shown in Table 5 for t equals 1 comma 2 comma ellipsis comma n and i equals 1 comma 2 comma ellipsis comma q asterisk p Subscript t Baseline.

Table 5: KFS: Filtering Phase

Quantity Description
ModifyingAbove y With caret Subscript t comma i Baseline equals normal upper E left-parenthesis y Subscript t comma i Baseline vertical-bar y Subscript t comma i minus 1 Baseline comma ellipsis y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis One-step-ahead prediction of the response values
nu Subscript t comma i Baseline equals y Subscript t comma i Baseline minus ModifyingAbove y With caret Subscript t comma i One-step-ahead prediction residuals
upper F Subscript t comma i Baseline equals normal upper V normal a normal r left-parenthesis y Subscript t comma i Baseline vertical-bar y Subscript t comma i minus 1 Baseline comma ellipsis y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis Variance of the one-step-ahead prediction
ModifyingAbove alpha alpha With caret Subscript t comma i Baseline equals normal upper E left-parenthesis alpha alpha Subscript t Baseline vertical-bar y Subscript t comma i minus 1 Baseline comma ellipsis y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis One-step-ahead prediction of the state vector
bold upper P Subscript t comma i Baseline equals normal upper C normal o normal v left-parenthesis alpha alpha Subscript t Baseline vertical-bar y Subscript t comma i minus 1 Baseline comma ellipsis y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis Covariance of ModifyingAbove alpha alpha With caret Subscript t comma i
bold b Subscript t comma i left-parenthesis d plus k plus g right-parenthesis-dimensional vector
bold upper S Subscript t comma i left-parenthesis d plus k plus g right-parenthesis-dimensional symmetric matrix
left-parenthesis ModifyingAbove delta delta With caret ModifyingAbove beta beta With caret ModifyingAbove gamma gamma With caret right-parenthesis Subscript t comma i Superscript prime Baseline equals bold upper S Subscript t comma i Superscript negative 1 Baseline bold b Subscript t comma i Estimates of delta delta, beta beta, and gamma gamma by using the data up to left-parenthesis t comma i right-parenthesis
bold upper S Subscript t comma i Superscript negative 1 Covariance of left-parenthesis ModifyingAbove delta delta With caret ModifyingAbove beta beta With caret ModifyingAbove gamma gamma With caret right-parenthesis Subscript t comma i Superscript prime
bold upper S Subscript t comma i Superscript asterisk left-parenthesis d plus k plus g right-parenthesis-dimensional symmetric matrix
needed in the marginal likelihood computation


Here the notation normal upper E left-parenthesis y Subscript t comma i Baseline vertical-bar y Subscript t comma i minus 1 Baseline comma ellipsis comma y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis denotes the conditional expectation of y Subscript t comma i given the history up to the index left-parenthesis t comma i minus 1 right-parenthesis: left-parenthesis y Subscript t comma i minus 1 Baseline comma ellipsis comma y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis. Similarly normal upper V normal a normal r left-parenthesis y Subscript t comma i Baseline vertical-bar y Subscript t comma i minus 1 Baseline comma ellipsis comma y Subscript t comma 1 Baseline comma bold upper Y Subscript t minus 1 Baseline comma ellipsis comma bold upper Y 1 right-parenthesis denotes the corresponding conditional variance. The quantity nu Subscript t comma i Baseline equals y Subscript t comma i Baseline minus ModifyingAbove y With caret Subscript t comma i is set to missing whenever y Subscript t comma i is missing. Note that ModifyingAbove y With caret Subscript t comma i are one-step-ahead forecasts only when the model has only one response variable and the data are a time series; in all other cases it is more appropriate to call them one-measurement-ahead forecasts (since the next measurement might be at the same time point). Despite this, ModifyingAbove y With caret Subscript t comma i are called one-step-ahead predictions (and nu Subscript t comma i are called one-step-ahead residuals) throughout this document. In the diffuse case, the conditional expectations must be appropriately interpreted. The vector bold b Subscript t comma i and the matrix bold upper S Subscript t comma i contain some accumulated quantities that are needed for the estimation of delta delta, beta beta, and gamma gamma. Of course, when left-parenthesis d plus k plus g right-parenthesis equals 0 (the nondiffuse case), these quantities are not needed. In the diffuse case, because the matrix bold upper S Subscript t comma i is sequentially accumulated (starting at t equals 1 comma i equals 1), it might not be invertible until some t equals t Subscript asterisk Baseline comma i equals i Subscript asterisk Baseline. The filtering process is called initialized after t equals t Subscript asterisk Baseline comma i equals i Subscript asterisk Baseline. In some situations, this initialization might not happen even after the entire sample is processed—that is, the filtering process remains uninitialized. This can happen if the regression variables are collinear or if the data are not sufficient to estimate the initial condition delta delta for some other reason. In the diffuse case if the marginal likelihood is to be computed (see the section Likelihood Computation and Model-Fitting Phase), an additional matrix (bold upper S Subscript t comma i Superscript asterisk) is computed at each step by sequential accumulation.

The filtering process is used for a variety of purposes. One important use of filtering is to compute the likelihood of the data. In the model-fitting phase, the unknown model parameters theta theta are estimated by maximum likelihood. This requires repeated evaluation of the likelihood at different trial values of theta theta. After theta theta is estimated, it is treated as a known vector. The filtering process is used again with the fitted model in the forecasting phase, when the one-step-ahead forecasts and residuals based on the fitted model are provided. In addition, this filtering output is needed by the smoothing phase to produce the full-sample component estimates and for the structural break analysis.

Last updated: June 19, 2025