ncvreg (version 3.1-0)

cv.ncvreg: Cross-validation for ncvreg

Description

Performs k-fold cross validation for MCP- or SCAD-penalized regression models over a grid of values for the regularization parameter lambda.

Usage

cv.ncvreg(X, y, ..., nfolds=10, seed, trace=FALSE)

Arguments

X

The design matrix, without an intercept, as in ncvreg.

y

The response vector, as in ncvreg.

...

Additional arguments to ncvreg.

nfolds

The number of cross-validation folds. Default is 10.

seed

You may set the seed of the random number generator in order to obtain reproducible results.

trace

If set to TRUE, cv.ncvreg will inform the user of its progress by announcing the beginning of each CV fold. Default is FALSE.

Value

• An object with S3 class "cv.ncvreg" containing:
• cve
  The error for each value of lambda, averaged across the cross-validation folds.
• cvse
  The estimated standard error associated with each value of for cve.
• lambda
  The sequence of regularization parameter values along which the cross-validation error was calculated.
• fit
  The fitted ncvreg object for the whole data.
• min
  The index of lambda corresponding to lambda.min.
• lambda.min
  The value of lambda with the minimum cross-validation error.
• null.dev
  The deviance for the intercept-only model.
• pe
  If family="binomial", the cross-validation prediction error for each value of lambda.

Details

The function calls ncvreg nfolds times, each time leaving out 1/nfolds of the data. The cross-validation error is based on the residual sum of squares when family="gaussian" and the binomial deviance when family="binomial".

References

Breheny, P. and Huang, J. (2011) Coordinate descent algorithms for nonconvex penalized regression, with applications to biological feature selection. Ann. Appl. Statist., 5: 232-253.

See Also

ncvreg, plot.cv.ncvreg, summary.cv.ncvreg

Examples

data(prostate)
X <- as.matrix(prostate[,1:8])
y <- prostate$lpsa

cvfit <- cv.ncvreg(X,y)
plot(cvfit)
summary(cvfit)

fit <- cvfit$fit
plot(fit)
beta <- fit$beta[,cvfit$min]