## Load data and select the deterministic components.
data("gergnp")
lgergnp <- log(gergnp, base=exp(1))
## Fit a PIAR(2) model with seasonal intercepts.
out.pred <- predictpiar(wts=lgergnp, p=2, hpred=24)
plotpredpiar(out.pred)
Run the code above in your browser using DataLab