nnbeta <- function(p, k) nbeta(c(1,p),k)
dgp <- midas_pl_sim(250, m = 12, theta = nnbeta(c(2, 4), 24),
gfun = function(x) 0.25*x^3,
ar.x = 0.9, ar.y = 0.5, n.start = 100)
Run the code above in your browser using DataLab