x=rnorm(1000)*0.1
y=matrix(0,1000,1)
y[1]=2
lambda=matrix(0,1000,1)
for (i in 2:1000){
lambda[i]=2+0.2*y[i-1]/exp(x[i-1])+0.5*lambda[i-1]
y[i]=rpois(1,exp(x[i])*lambda[i])
}
ACMx(y,order=c(1,1),x,"po")
Run the code above in your browser using DataLab