##The parameter function
theta.h0 <- function(p, dk) {
i <- (1:dk-1)/100
pol <- p[3]*i + p[4]*i^2
(p[1] + p[2]*i)*exp(pol)
}
##Generate coefficients
theta0 <- theta.h0(c(-0.1,10,-10,-10),4*12)
##Plot the coefficients
plot(theta0)
##Generate the predictor variable
x <- simplearma.sim(list(ar=0.6),1500*12,1,12)
##Simulate the response variable
y <- midas.sim(500,theta0,x,1)
Run the code above in your browser using DataLab