我想知道如何从cv.glmnet对象中提取特定lambda的拟合值,残差和摘要统计信息>(例如“ lambda.1se”)。假设只有我可以访问cv.glmnet对象,而不能直接访问训练数据。这是一个示例:
data(iris) Inx <- sample(nrow(iris),100) iris$Species <- factor(iris$Species) train_data <- iris[Inx, ] test_data <- iris[-Inx,] Formula <- "Sepal.Length ~ Sepal.Width + Petal.Length + Species:Petal.Width + Sepal.Width:Petal.Length + Species + splines::bs(Petal.Width, df = 2, degree = 2)"
--------- glm交叉验证--------------
ModelMatrix <- predict(caret::dummyVars(Formula, train_data, fullRank = T, sep = ""), train_data) cvglm <- glmnet::cv.glmnet(x = ModelMatrix, y = train_data$Sepal.Length, nfolds = 4, keep = TRUE, alpha = 1, parallel = F, type.measure = 'mse')
谢谢您的任何建议。
[我想知道如何从cv.glmnet对象中提取特定lambda(例如“ lambda.1se”)的拟合值,残差和摘要统计信息。假设只有我可以访问cv.glmnet对象...
拟合值: