This function propagates the state trajectory based on the exogenous inputs only (without measurement update), and calculates the corresponding log-likelihood
propagate(theta, u, v, y, stdlik = TRUE)
A list of system parameters (A, B, C, D, Q, R)'
Input matrix for the state equation (m_u rows, T columns)
Input matrix for the output equation (m_v rows, T columns)
Observations
Boolean, whether the likelihood is divided by the number of observations. Standardizing the likelihood this way may speed up convergence in the case of long time series.
A list of predictions and log-likelihood (X, Y, V, lik)