bvartools (version 0.2.4)

kalman_dk: Durbin and Koopman Simulation Smoother

Description

An implementation of the Kalman filter and backward smoothing algorithm proposed by Durbin and Koopman (2002).

Usage

kalman_dk(y, z, sigma_u, sigma_v, B, a_init, P_init)

Value

A \(M \times T+1\) matrix of state vector draws.

Arguments

y

a \(K \times T\) matrix of endogenous variables.

z

a \(KT \times M\) matrix of explanatory variables.

sigma_u

the constant \(K \times K\) error variance-covariance matrix. For time varying variance-covariance matrices a \(KT \times K\) can be specified.

sigma_v

the constant \(M \times M\) coefficient variance-covariance matrix. For time varying variance-covariance matrices a \(MT \times M\) can be specified.

B

an \(M \times M\) autocorrelation matrix of the transition equation.

a_init

an M-dimensional vector of initial states.

P_init

an \(M \times M\) variance-covariance matrix of the initial states.

Details

The function uses algorithm 2 from Durbin and Koopman (2002) to produce a draw of the state vector \(a_t\) for \(t = 1,...,T\) for a state space model with measurement equation $$y_t = Z_t a_t + u_t$$ and transition equation $$a_{t + 1} = B_t a_{t} + v_t,$$ where \(u_t \sim N(0, \Sigma_{u,t})\) and \(v_t \sim N(0, \Sigma_{v,t})\). \(y_t\) is a K-dimensional vector of endogenous variables and \(Z_t = z_t^{\prime} \otimes I_K\) is a \(K \times M\) matrix of regressors with \(z_t\) as a vector of regressors.

The algorithm takes into account Jarociński (2015), where a possible missunderstanding in the implementation of the algorithm of Durbin and Koopman (2002) is pointed out. Following that note the function sets the mean of the initial state to zero in the first step of the algorithm.

References

Durbin, J., & Koopman, S. J. (2002). A simple and efficient simulation smoother for state space time series analysis. Biometrika, 89(3), 603--615.

Jarociński, M. (2015). A note on implementing the Durbin and Koopman simulation smoother. Computational Statistics and Data Analysis, 91, 1--3. tools:::Rd_expr_doi("10.1016/j.csda.2015.05.001")

Examples

# Load data
data("e1")
data <- diff(log(e1))

# Generate model data
temp <- gen_var(data, p = 2, deterministic = "const")
y <- t(temp$data$Y)
z <- temp$data$SUR
k <- nrow(y)
tt <- ncol(y)
m <- ncol(z)

# Priors
a_mu_prior <- matrix(0, m)
a_v_i_prior <- diag(0.1, m)

a_Q <- diag(.0001, m)

# Initial value of Sigma
sigma <- tcrossprod(y) / tt
sigma_i <- solve(sigma)

# Initial values for Kalman filter
y_init <- y * 0
a_filter <- matrix(0, m, tt + 1)

# Initialise the Kalman filter
for (i in 1:tt) {
  y_init[, i] <- y[, i] - z[(i - 1) * k + 1:k,] %*% a_filter[, i]
}
a_init <- post_normal_sur(y = y_init, z = z, sigma_i = sigma_i,
                          a_prior = a_mu_prior, v_i_prior = a_v_i_prior)
y_filter <- matrix(y) - z %*% a_init
y_filter <- matrix(y_filter, k) # Reshape

# Kalman filter and backward smoother
a_filter <- kalman_dk(y = y_filter, z = z, sigma_u = sigma,
                      sigma_v = a_Q, B = diag(1, m),
                      a_init = matrix(0, m), P_init = a_Q)
                      
a <- a_filter + matrix(a_init, m, tt + 1)