plot.mle.cp: Plot the Mallows Cp

In wle: Weighted Likelihood Estimation

Description

Plot the Mallows Cp.

Usage

## S3 method for class 'mle.cp'
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.

Author(s)

Claudio Agostinelli

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)

wle documentation built on May 29, 2017, 11:48 a.m.