wle (version 0.5)

plot.mle.cp: Plot the Mallows Cp

Description

Plot the Mallows Cp.

Usage

plot.mle.cp(object, base.line=0, num.max=20, 
           plot.it=TRUE, log.scale=FALSE, 
           xlab="Number of Predictors", ylab=NULL)

Arguments

object

an object of class mle.cp.

base.line

the intercept of the line to split the submodels in acceptable (good) and not-acceptable (bad), (the slope is always one).

num.max

maximum number of submodels plotted.

plot.it

if TRUE the graph is plotted.

log.scale

if TRUE the y-axis as log10 scale.

xlab

a title for the x axis.

ylab

a title for the y axis.

Value

• num.good
  number of submodels below the base.line
• num.bad
  number of submodels above the base.line
• cp.good
  list of the submodels below the base.line with their Cp.
• cp.bad
  list of the submodels above the base.line with their Cp.

See Also

mle.cp a function to calculate the Mallows Cp.

Examples

library(wle)

data(hald)

result <- mle.cp(y.hald~x.hald)

plot(result,num.max=7)