powered by
This function is adapted from the plot.regsubsets function of the leaps package, and its main use is to plot the output from the armasubsets function.
# S3 method for armasubsets plot(x, labels = obj$xnames, main = NULL, scale = c("BIC", "AICc", "AIC", "Cp", "adjR2", "R2"), col = gray(c(seq(0.4, 0.7, length = 10), 0.9)), draw.grid = TRUE, axis.at.3 = TRUE, ...)
an object of class armasubsets
variable names
title for plot
which summary statistic to use for ordering plots
the last color should be close to but distinct from white
a logical argument; if it is true (default), gray grid lines are superimposed on the graph.
a logical argument; if if it is true (default), the x-labels are drawn on the upper horizontal axis.
other arguments
Plot the few best subset ARMA models.
armasubsets
# NOT RUN { set.seed(53331) test=arima.sim(model=list(ar=c(rep(0,11),.8),ma=c(rep(0,11),0.7)),n=120) res=armasubsets(y=test,nar=14,nma=14,y.name='test',ar.method='ols') plot(res) # }
Run the code above in your browser using DataLab