# NOT RUN {
# runs longly
library(MTS)
set.seed(1)
yt <- VARMAsim(300,arlags=1,malags=1,
phi=matrix(c(0.2,-0.6,0.3,1.1),2,2),
theta=matrix(c(-0.5,0,0,-0.5),2,2),
sigma=diag(2))$series
M <- VARMA(yt,p=2,q=2,include.mean=TRUE)
M[-c(1,6,7,9,10)]
pMvarma<-MTS2pM(M)
print(pMvarma,style="broad",digits=3)
pMvarma
print(pMvarma,"matrix",digits=3)
# clean up
# rm(yt, M, pMvarma)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab