regsubsets function in the leaps package finds
  optimal subsets of predictors. This function plots a measure of fit
  (see the statistic argument below) against subset size).subsets(object, ...)
## S3 method for class 'regsubsets':
subsets(object, 
    names=abbreviate(object$xnames, minlength = abbrev), 
    abbrev=1, min.size=1, max.size=length(names), legend, 
    statistic=c("bic", "cp", "adjr2", "rsq", "rss"), 
    las=par('las'), cex.subsets=1, ...)regsubsets object produced by the regsubsets function
    in the leaps package.object."bic", Bayes Information Criterion; 
    "cp", Mallows's $C_{p}$;
    "adjr2", $R^{2}$ adjusted for degrees of freedom;
    "rsq", unadjusted0, ticks labels are drawn parallel to the
    axis; set to 1 for horizontal labels (see par).1.subsets.regsubsets and plot.NULL. This function is used for its side effect --
  to create a plot.regsubsetslibrary(leaps)
subsets(regsubsets(undercount ~ ., data=Ericksen))Run the code above in your browser using DataLab