tsdiagplot(x, p.max=2, q.max=p.max, model=c(p.max, 0, q.max), ## S-Plus order=c(p.max, 0, q.max), ## R lag.max=36, gof.lag=lag.max, armas=arma.loop(x, order=order, series=deparse(substitute(x)), ...), diags=diag.arma.loop(armas, x, lag.max=lag.max, gof.lag=gof.lag), ts.diag=rearrange.diag.arma.loop(diags), lag.units=ts.diag$tspar["frequency"], lag.lim=range(pretty(ts.diag$acf$lag))*lag.units, lag.x.at=pretty(ts.diag$acf$lag)*lag.units, lag.x.labels={tmp <- lag.x.at tmp[as.integer(tmp)!=tmp] <- "" tmp}, lag.0=TRUE, main, lwd=0, ...)
acfplot(rdal, type="acf", main=paste("ACF of std.resid:", rdal$series, " model:", rdal$model), lag.units=rdal$tspar["frequency"], lag.lim=range(pretty(rdal[[type]]$lag)*lag.units), lag.x.at=pretty(rdal[[type]]$lag)*lag.units, lag.x.labels={tmp <- lag.x.at tmp[as.integer(tmp)!=tmp] <- "" tmp}, lag.0=TRUE, xlim=xlim.function(lag.lim/lag.units), ...)
aicsigplot(z, z.name=deparse(substitute(z)), series.name="ts", model=NULL, xlab="", ylab=z.name, main=paste(z.name, series.name, model), layout=c(1,2), between=list(x=1,y=1), ...)
residplot(rdal, main=paste("std.resid:", rdal$series, " model:", rdal$model), ...)
gofplot(rdal, main=paste("P-value for gof:", rdal$series, " model:", rdal$model), lag.units=rdal$tspar["frequency"], lag.lim=range(pretty(rdal$gof$lag)*lag.units), lag.x.at=pretty(rdal$gof$lag)*lag.units, lag.x.labels={tmp <- lag.x.at tmp[as.integer(tmp)!=tmp] <- "" tmp}, xlim=xlim.function(lag.lim/lag.units), pch=16, ...)