Hi,

I am trying to use the simulate function in the dse1 package to 
generate VAR (multivariate) models. I *think* the first example gives a 
simple VAR(2) time series, here is the code.....

 > AR   <- array(c(1, .5, .3, 0, .2, .1, 0, .2, .05, 1, .5, .3) 
,c(3,2,2))
 > VAR  <- ARMA(A=AR, B=diag(1,2))
 > print(VAR)

A(L) =
1+0.5L1+0.3L2    0+0.2L1+0.05L2
0+0.2L1+0.1L2    1+0.5L1+0.3L2

B(L) =
1    0
0    1

 > simData <- simulate(VAR)

 > simData$output[1:4,]
            [,1]       [,2]
[1,]  0.3153194 -1.3748001
[2,] -0.6794556  1.0139211
[3,] -0.9385841 -0.4116173
[4,]  1.3702363  0.3295771

However, I really do not understand what is going on here. I can see 
that simulate has returned two output time series y_{1,t} and y_{2,t}.

What univariate models are y_{1,t} and y_{2,t} and how is the 
relationship (i.e. multivariate) between the two series defined?

Any help would be greatly appreciated.

Best regards,

Sam.
	[[alternative text/enriched version deleted]]


