library(lars)
ans <- cv.lars(data, K = 10, index, trace = FALSE, plot.it = TRUE,type = "lasso")
dev.new()
plot(ans)
不过,不用lars呢~~~
library(glmnet)
ans <- cv.glmnet(data, return[,1], standardize = TRUE)
par(mfrow=c(1, 2))
plot(ans$glmnet.fit, "norm", label=TRUE)
plot(ans$glmnet.fit, "lambda", label=TRUE)
dev.new()
plot(ans)