data(llynx) plotts.wge(llynx)
fore.arma.wge(llynx, phi = c(0.7, 0.1, -0.2, -0.3), theta = -0.6, n.ahead = 20, limits = F)