在插入符号中获取交叉验证的 glmnet 模型的系数

Question

如何在 R 的 caret 包中获得交叉验证的 glmnet 模型的系数？

library(caret)
x = matrix(rnorm(500), ncol=10)
y = rnorm(100)
control = trainControl(
         method="LOOCV",
         allowParallel = TRUE,
         number = nrow(sample_df), 
         verboseIter = FALSE,
         returnData = FALSE
        )

  my_elasticnet <- train(x, y,
                     method = "glmnet",
                     preProc = c("center", "scale"),
                     trControl = control)

Answer 1

我想你忘了显示 sample_df 的代码。但是，假设以下情况，您可以按如下方式访问它：

library(caret)
x = matrix(rnorm(500), ncol=10)
y = rnorm(100)
sample_df = cbind.data.frame(y,x)
control = trainControl(
  method="LOOCV",
  allowParallel = TRUE,
  number = nrow(sample_df), 
  verboseIter = FALSE,
  returnData = FALSE
)

my_elasticnet <- train(sample_df[2:11], sample_df$y,
                       method = "glmnet",
                       preProc = c("center", "scale"),
                       trControl = control)

my_elasticnet$finalModel$beta

如果您只看名称，您将获得有关最终模型所需的一切：

> names(my_elasticnet$finalModel)
 [1] "a0"          "beta"        "df"          "dim"         "lambda"      "dev.ratio"   "nulldev"    
 [8] "npasses"     "jerr"        "offset"      "call"        "nobs"        "lambdaOpt"   "xNames"     
[15] "problemType" "tuneValue"   "obsLevels"   "param"    

编辑：回应评论

最终模型取决于您选择的 alpha 和 lambda 级别。有 66 个这样的参数。如果你想选择一台机器认为最好的，你可以这样做：

coef(my_elasticnet$finalModel, my_elasticnet$bestTune$lambda)

这将为您提供您正在寻找的 11x1 矢量。