### R code from vignette source 'partsm.Rnw'

###################################################
### code chunk number 1: partsm.Rnw:164-167
###################################################
library(partsm)
data("gergnp")
lgergnp <- log(gergnp, base=exp(1))


###################################################
### code chunk number 2: partsm.Rnw:178-188
###################################################
detcomp <- list(regular=c(0,0,0), seasonal=c(1,0), regvar=0)
aic <- bic <- Fnextp <- Fpval <- rep(NA, 4)
for(p in 1:4){
  lmpar <- fit.ar.par(wts=lgergnp, detcomp=detcomp, type="PAR", p=p)
  aic[p] <- AIC(lmpar@lm.par, k=2)
  bic[p] <- AIC(lmpar@lm.par, k=log(length(residuals(lmpar@lm.par))))
  Fout <- Fnextp.test(wts=lgergnp, detcomp=detcomp, p=p, type="PAR")
  Fnextp[p] <- Fout@Fstat
  Fpval[p] <- Fout@pval
}


###################################################
### code chunk number 3: partsm.Rnw:217-223
###################################################
dcsi <- list(regular=c(0,0,0), seasonal=c(1,0), regvar=0)
out.Fparsi  <- Fpar.test(wts=lgergnp, detcomp=dcsi, p=2)
show(out.Fparsi)
dcsit <- list(regular=c(0,0,0), seasonal=c(1,1), regvar=0)
out.Fparsit <- Fpar.test(wts=lgergnp, detcomp=dcsit, p=2)
show(out.Fparsit)


###################################################
### code chunk number 4: partsm.Rnw:232-235
###################################################
par2 <- fit.ar.par(wts=lgergnp, type="PAR", p=2, detcomp=detcomp)
Fsh.out <- Fsh.test(res=residuals(par2@lm.par), s=frequency(lgergnp))
show(Fsh.out)


###################################################
### code chunk number 5: partsm.Rnw:247-250
###################################################
out.par <- fit.ar.par(wts=lgergnp, type="PAR", detcomp=detcomp, p=2)
out.MV <- PAR.MVrepr(out.par)
out.MV


###################################################
### code chunk number 6: partsm.Rnw:272-274
###################################################
out.LR <- LRurpar.test(wts=lgergnp, detcomp=detcomp, p=2)
show(out.LR)


###################################################
### code chunk number 7: partsm.Rnw:279-281
###################################################
Fpari1.out <- Fpari.piar.test(wts=lgergnp, detcomp=detcomp, p=2, type="PARI1")
show(Fpari1.out)


###################################################
### code chunk number 8: partsm.Rnw:337-340
###################################################
out.piar <- fit.piar(wts=lgergnp, detcomp=detcomp, p=2)
out.MV <- PAR.MVrepr(out.piar)
out.MV


###################################################
### code chunk number 9: pdiffplot
###################################################
plotpdiff(out.piar)


###################################################
### code chunk number 10: partsm.Rnw:362-364
###################################################
out.pred <- predictpiar(wts=lgergnp, p=2, hpred=24)
show(out.pred)


###################################################
### code chunk number 11: partsm.Rnw:369-374
###################################################
out.pred@wts <- exp(1)^out.pred@wts
out.pred@fcast <- exp(1)^out.pred@fcast
out.pred@ucb <- exp(1)^out.pred@ucb
out.pred@lcb <- exp(1)^out.pred@lcb
plotpredpiar(out.pred)


###################################################
### code chunk number 12: predplot
###################################################
plotpredpiar(out.pred)