lmSubsets (version 0.5-1)

image: Variable Selection Heatmaps

Description

Visualization of variable subsets.

Usage

# S3 method for lmSubsets
image(x, size = NULL, best = 1, which = NULL,
      hilite, hilite_penalty, main, sub, xlab = NULL, ylab,
      ann = par("ann"), axes = TRUE, col = c("gray40", "gray90"),
      lab = "lab", col_hilite = cbind("red", "pink"),
      lab_hilite = "lab", pad_size = 3, pad_best = 1,
      pad_which = 3, axis_pos = -4, axis_tck = -4,
      axis_lab = -10, ...)
# S3 method for lmSelect
image(x, best = NULL, which = NULL, hilite,
      hilite_penalty, main, sub = NULL, xlab = NULL, ylab,
      ann = par("ann"), axes = TRUE, col = c("gray40", "gray90"),
      lab = "lab",  col_hilite = cbind("red", "pink"),
      lab_hilite = "lab", pad_best = 2, pad_which = 2,
      axis_pos = -4, axis_tck = -4, axis_lab = -10, ...)

Arguments

x

An object of class "lmSubsets" or "lmSelect".

main, sub, xlab, ylab

Main, sub- and axis titles.

size, best

Submodels to be plotted.

which

Regressors to be plotted.

col, lab

Color and label style.

hilite, hilite_penalty

Submodels to be highlighted.

col_hilite, lab_hilite

Highlighting style.

pad_size, pad_best, pad_which

Padding.

axis_pos, axis_tck, axis_lab

Position of axes, tick length, and position of labels

…

Ignored.

axes

Plot axes.

ann

Annotate plot.

See Also

lmSubsets, lmSelect.

Examples

# NOT RUN {
## data
data("AirPollution", package = "lmSubsets")


#################
##  lmSubsets  ##
#################

lm_all <- lmSubsets(mortality ~ ., data = AirPollution, nbest = 20)

## heatmap
image(lm_all, best = 1:3)

## highlight 5 best (BIC)
image(lm_all, best = 1:3, hilite = 1:5, hilite_penalty = "BIC")


################
##  lmSelect  ##
################

## default criterion: BIC
lm_best <- lmSelect(lm_all)

## highlight 5 best (AIC)
image(lm_best, hilite = 1:5, hilite_penalty = "AIC")

## axis labels
image(lm_best, lab = c("bold(lab)", "lab"), hilite = 1,
      lab_hilite = "underline(lab)")
# }

Run the code above in your browser using DataCamp Workspace