DYNAMICLINEAR Procedure

Simultaneous Graphical Dynamic Linear Models

The classical dynamic linear model (DLM) often falters when confronted with high-dimensional time series that are characterized by intricate correlations. Recognizing this limitation, Gruber and West (2016) and Gruber and West (2017) introduced an enhanced version known as a simultaneous graphical dynamic linear model (SGDLM). This advanced model is adept at efficiently managing large-scale, correlated time series.

Consider an m-dimensional time series that is denoted by bold y Subscript t Baseline equals left parenthesis y Subscript 1 t Baseline comma y Subscript 2 t Baseline comma ellipsis comma y Subscript m t Baseline right parenthesis prime. The SGDLM can be expressed as

In this representation, nu Subscript j t Baseline tilde upper N left parenthesis 0 comma lamda Subscript j t Superscript negative 1 Baseline right parenthesis and bold w Subscript j t Baseline tilde upper N left parenthesis bold 0 comma bold upper W Subscript j t Baseline right parenthesis. The expression bold y Subscript normal s normal p left parenthesis j right parenthesis comma t denotes a vector that represents values of other series indexed by normal s normal p left parenthesis j right parenthesis subset of or equal to 1 colon m minus j, where normal s normal p left parenthesis j right parenthesis denotes the simultaneous parental set of series j. The vector bold x Subscript j t Baseline equals left parenthesis x Subscript j comma 1 comma t Baseline comma x Subscript j comma 2 comma t Baseline comma ellipsis right parenthesis prime denotes the vector of external predictors for series j at time t; it is important to note that bold x Subscript j t might vary in dimensions across the m variables.

First, here are the definitions of certain matrices and vectors that are related to SGDLMs:

  • Let bold upper Theta Subscript t represent the collection of bold italic theta Subscript j t values for every j, formally written as

    bold upper Theta Subscript t Baseline equals StartSet bold italic theta Subscript 1 t Baseline comma bold italic theta Subscript 2 t Baseline comma ellipsis comma bold italic theta Subscript m t Baseline EndSet
  • Define bold upper Lamda Subscript t as a diagonal matrix that contains lamda Subscript j t values:

    bold upper Lamda Subscript t Baseline equals normal d normal i normal a normal g StartSet lamda Subscript 1 t Baseline comma lamda Subscript 2 t Baseline comma ellipsis comma lamda Subscript m t Baseline EndSet
  • bold italic mu Subscript t is the vector given by

    bold italic mu Subscript t Baseline equals left parenthesis mu Subscript 1 t Baseline comma mu Subscript 2 t Baseline comma ellipsis comma mu Subscript m t Baseline right parenthesis prime

    where the individual components are defined as bold italic mu Subscript j t Baseline equals bold x prime Subscript j t Baseline bold italic phi Subscript j t.

The matrix bold upper Gamma Subscript t captures the relationships between various time series. It is defined as

bold upper Gamma Subscript t Baseline equals Start 4 By 1 Matrix 1st Row  bold italic gamma Subscript 1 t Superscript prime Baseline 2nd Row  bold italic gamma Subscript 2 t Superscript prime Baseline 3rd Row  vertical ellipsis 4th Row  bold italic gamma Subscript m t Superscript prime Baseline EndMatrix equals Start 4 By 5 Matrix 1st Row 1st Column 0 2nd Column gamma Subscript 1 comma 2 comma t Baseline 3rd Column gamma Subscript 1 comma 3 comma t Baseline 4th Column ellipsis 5th Column gamma Subscript 1 comma m comma t Baseline 2nd Row 1st Column gamma Subscript 2 comma 1 comma t Baseline 2nd Column 0 3rd Column gamma Subscript 2 comma 3 comma t Baseline 4th Column ellipsis 5th Column gamma Subscript 2 comma m comma t Baseline 3rd Row 1st Column vertical ellipsis 2nd Column vertical ellipsis 3rd Column vertical ellipsis 4th Column down right diagonal ellipsis 5th Column vertical ellipsis 4th Row 1st Column gamma Subscript m comma 1 comma t Baseline 2nd Column gamma Subscript m comma 2 comma t Baseline 3rd Column ellipsis 4th Column gamma Subscript m comma m minus 1 comma t Baseline 5th Column 0 EndMatrix

where gamma Subscript j comma h comma t Baseline equals 0 for each h not an element of normal s normal p left parenthesis j right parenthesis comma j equals 1 colon m.

Use these definitions to represent the distribution of bold y Subscript t as

where the matrices bold upper A Subscript t and bold upper Sigma Subscript t are given by

and where bold upper Omega Subscript t is the precision matrix.

Finally, PROC DYNAMICLINEAR integrates the local-level dynamic linear model with the simultaneous graphical feature, as follows:

where bold upper F Subscript j t Baseline equals left parenthesis 1 comma bold y Subscript normal s normal p left parenthesis j right parenthesis comma t Baseline right parenthesis prime for all j, and bold upper G Subscript j t Baseline equals bold upper I for all j and t.

Forward Filtering in SGDLMs

To understand the forward filtering process in an SGDLM, you need to understand the distributions at play. Here are the key distributions that this process uses. For more information about these distributions, see Gruber and West (2016) and Gruber and West (2017).

  • The prior at time t is given by

    normal p left parenthesis bold upper Theta Subscript t Baseline comma bold upper Lamda Subscript t Baseline vertical bar script upper I Subscript t Baseline comma script upper D Subscript t minus 1 Baseline right parenthesis equals product Underscript j equals 1 colon m Endscripts normal p Subscript j t Baseline left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper I Subscript t Baseline comma script upper D Subscript t minus 1 Baseline right parenthesis

    where script upper I Subscript t denotes the information set at time t; script upper D Subscript t minus 1 denotes the data set up to time t minus 1; and normal p left parenthesis dot vertical bar dot right parenthesis denotes the density function of

    left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper I Subscript t Baseline comma script upper D Subscript t minus 1 Baseline right parenthesis tilde normal upper N normal upper G left parenthesis bold italic a Subscript j t Baseline comma bold upper R Subscript j t Baseline comma r Subscript j t Baseline comma c Subscript j t Baseline right parenthesis

    and

    StartLayout 1st Row 1st Column left parenthesis bold italic theta Subscript j t Baseline vertical bar lamda Subscript j t Baseline comma script upper I Subscript t Baseline comma script upper D Subscript t minus 1 Baseline right parenthesis 2nd Column tilde normal upper N left parenthesis bold italic a Subscript j t Baseline comma bold upper R Subscript j t Baseline divided by left parenthesis c Subscript j t Baseline lamda Subscript j t Baseline right parenthesis right parenthesis 2nd Row 1st Column left parenthesis lamda Subscript j t Baseline vertical bar script upper I Subscript t Baseline comma script upper D Subscript t minus 1 Baseline right parenthesis 2nd Column tilde normal upper G normal a normal m normal m normal a left parenthesis r Subscript j t Baseline divided by 2 comma r Subscript j t Baseline c Subscript j t Baseline divided by 2 right parenthesis EndLayout

    where normal upper N left parenthesis bold italic a comma bold upper A right parenthesis is a multivariate normal distribution that has the mean bold italic a and the covariance matrix bold upper A, and normal upper G normal a normal m normal m normal a left parenthesis alpha comma beta right parenthesis is a gamma distribution that has the shape parameter alpha and the rate parameter beta.

  • The naive posterior approximation at time t to the true posterior normal p left parenthesis bold upper Theta Subscript t Baseline comma bold upper Lamda Subscript t Baseline vertical bar script upper D Subscript t Baseline right parenthesis is given by

    ModifyingAbove normal p With tilde left parenthesis bold upper Theta Subscript t Baseline comma bold upper Lamda Subscript t Baseline vertical bar script upper D Subscript t Baseline right parenthesis equals product Underscript j equals 1 colon m Endscripts ModifyingAbove normal p With tilde Subscript j t Baseline left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper D Subscript t Baseline right parenthesis

    where

    left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper D Subscript t Baseline right parenthesis tilde normal upper N normal upper G left parenthesis bold italic m overtilde Subscript j t Baseline comma bold upper C overtilde Subscript j t Baseline comma n overtilde Subscript j t Baseline comma s overtilde Subscript j t Baseline right parenthesis
  • The true posterior at time t is

    normal p left parenthesis bold upper Theta Subscript t Baseline comma bold upper Lamda Subscript t Baseline vertical bar script upper D Subscript t Baseline right parenthesis proportional to StartAbsoluteValue upper I minus bold upper Gamma Subscript t Baseline EndAbsoluteValue product Underscript j equals 1 colon m Endscripts ModifyingAbove normal p With tilde Subscript j t Baseline left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper D Subscript t Baseline right parenthesis proportional to product Underscript j equals 1 colon m Endscripts normal p Subscript j t Baseline left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper D Subscript t Baseline right parenthesis

    with the components

    left parenthesis bold italic theta Subscript j t Baseline comma lamda Subscript j t Baseline vertical bar script upper D Subscript t Baseline right parenthesis tilde normal upper N normal upper G left parenthesis bold italic m Subscript j t Baseline comma bold upper C Subscript j t Baseline comma n Subscript j t Baseline comma s Subscript j t Baseline right parenthesis
Last updated: July 09, 2026