require(dlmodeler)
# create the following model:
# deterministic level + quarterly seasonal + disturbance
mod1 <- dlmodeler.build.polynomial(0,sigmaH=.1) +
4* dlmodeler.build.dseasonal(4,sigmaH=0)
# create a multivariate model by binding the previous model
# with a stochastic trend model
mod2 <- mod1 %% dlmodeler.build.polynomial(1,sigmaH=0)
Run the code above in your browser using DataLab