sealasso: Standard Error Adjusted Adaptive Lasso

In sealasso: Standard Error Adjusted Adaptive Lasso

Standard Error Adjusted Adaptive Lasso

Description

This provides the solution path of some variants of the adaptive lasso (e.g. SEA-lasso, NSEA-lasso, OLS-adaptive lasso), and the corresponding optimal model selected by BIC criterion.

Usage

sealasso(x, y, method = c("nsealasso", "sealasso", "olsalasso", "lasso"))

Arguments

x	The model matrix.
y	The response.
method	One of "nsealasso", "sealasso", "olsalasso" and "lasso", which represent NSEA-lasso, SEA-lasso, OLS-adaptive lasso and the lasso, respectively. The default is "nsealasso".

Details

SEA-lasso and NSEA-lasso (Qian and Yang, 2010) are two versions of the adaptive lasso. They may be used for variable selection, especially in cases where condition index of the scaled model matrix is large (e.g. > 10) and collinearity is a concern. This function provides condition index, solution path and the suggested optimal model based on BIC. The estimated coefficients are also given for transition points of the path.

Value

method	The method used.
condition.index	Condition index of the scaled model matrix.
path	Solution path and corresponding BIC values at transition points.
beta	The estimated coefficients at transition points of solution path.
optim.beta	The estimated coefficients of the optimal model based on BIC criterion.

References

Qian, W. and Yang, Y. (2010) "Model Selection via Standard Error Adjusted Adaptive Lasso." Technical Report, University of Minnesota.

Examples

  # use the diabetes dataset from "lars" package
  data(diabetes)
  x <- diabetes$x
  y <- diabetes$y
  sealasso(x, y)

  # with quadratic terms
  x2 <- cbind(diabetes$x1,diabetes$x2)
  object <- sealasso(x2, y, "sealasso")
  object$condition.index
  object$optim.beta

sealasso documentation built on April 23, 2022, 1:05 a.m.