作者热门文章
- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我确信在运行 cva.glmnet 后有一种优雅的方法可以提取最好的 alpha 和 lambda,但不知何故我找不到它。
这是我在此期间使用的代码。
谢谢
library(data.table);library(glmnetUtils);library(useful)
# make some dummy data
data(iris)
x <- useful::build.x(data = iris,formula = Sepal.Length ~ .)
y <- iris$Sepal.Length
# run cv for alpha in c(0,0.5,1)
output.of.cva.glmnet <- cva.glmnet(x=x,y=y,alpha = c(0,0.5,1))
# extract the best parameters
number.of.alphas.tested <- length(output.of.cva.glmnet$alpha)
cv.glmnet.dt <- data.table()
for (i in 1:number.of.alphas.tested){
glmnet.model <- output.of.cva.glmnet$modlist[[i]]
min.mse <- min(glmnet.model$cvm)
min.lambda <- glmnet.model$lambda.min
alpha.value <- output.of.cva.glmnet$alpha[i]
new.cv.glmnet.dt <- data.table(alpha=alpha.value,min_mse=min.mse,min_lambda=min.lambda)
cv.glmnet.dt <- rbind(cv.glmnet.dt,new.cv.glmnet.dt)
}
best.params <- cv.glmnet.dt[which.min(cv.glmnet.dt$min_mse)]
最佳答案
基于我在 GitHub 上读到的一个线程,作者希望人们使用 plot(fit)
而不是仅仅输出最好的参数。但是,这并不总是可行的,尤其是在涉及交叉验证时。这些辅助函数可能是一个很好的解决方法。
# Train model.
fit <- cva.glmnet(X, y)
# Get alpha.
get_alpha <- function(fit) {
alpha <- fit$alpha
error <- sapply(fit$modlist, function(mod) {min(mod$cvm)})
alpha[which.min(error)]
}
# Get all parameters.
get_model_params <- function(fit) {
alpha <- fit$alpha
lambdaMin <- sapply(fit$modlist, `[[`, "lambda.min")
lambdaSE <- sapply(fit$modlist, `[[`, "lambda.1se")
error <- sapply(fit$modlist, function(mod) {min(mod$cvm)})
best <- which.min(error)
data.frame(alpha = alpha[best], lambdaMin = lambdaMin[best],
lambdaSE = lambdaSE[best], eror = error[best])
}
关于r - 从 cva.glmnet 对象中提取最佳参数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/54803990/
我想在弹性网回归模型中强制一些自变量。使用时出现错误惩罚因子在 cva.glmnet() 函数中。任何建议和意见表示赞赏。 # --------------- import data # outcom
我确信在运行 cva.glmnet 后有一种优雅的方法可以提取最好的 alpha 和 lambda,但不知何故我找不到它。 这是我在此期间使用的代码。 谢谢 library(data.table);l
我是一名优秀的程序员,十分优秀!