# NOT RUN {
# set dimensions (p=covariates, n=individuals, T=time points)
p <- 3; n <- 12; T <- 10
# set model parameters
SigmaE <- diag(p)/4
Ax <- createA(3, "chain")
# generate time-varying covariate data
X <- dataVAR1(n, T, Ax, SigmaE)
# autoregression parameter matrices of VARX(1) model
A <- createA(p, topology="clique", nonzeroA=0.1, nClique=1)
B <- createA(p, topology="hub", nonzeroA=0.1, nHubs=1)
# generate data
Y <- dataVARX1(X, A, B, SigmaE, lagX=0)
# }
Run the code above in your browser using DataLab