wle (version 0.9-91)

plot.mle.cp: Plot the Mallows Cp

Description

Plot the Mallows Cp.

Usage

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

Arguments

x

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.

verbose

if TRUE warnings are printed.

...

graphical parameters can be given as arguments.

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)