n <- 200
d <- 4
valinit <- 100
n2 <- n + d + 1
omega <- 0.05
alpha <- 0.05
beta <- 0.05
delta <- 2
pi <- rep(0.05, d)
e<-rnorm(n2+valinit)
Y<-e
for (t in 2:n2)
Y[t]<- 0.2*Y[t-1]+e[t]
x<-exp(Y)
X <- matrix(0, nrow = (n+valinit), ncol = length(pi))
for(j in 1:d)
X[, j] <- x[(d+2-j):(n+d+1-j+valinit)]
data <- simulate(n, omega, alpha, beta, delta, X, pi, valinit = valinit)
Run the code above in your browser using DataLab