如何在 h2o.automl 的排行榜中打印所有模型的可变重要性
how to print variable importance of all the models in the leaderboard of h2o.automl in r
我正在使用 R 中 H2o
包的 automl()
函数进行回归。
考虑一下我使用 "aml" 这个名字来构建模型。
aml <- h2o.automl(x=x, y=y, training_frame = train_set,
max_models = 20, seed = 1,
keep_cross_validation_predictions = TRUE)
automl()
的排行榜显示了表现最好的模型。我能够通过 h2o.varimp() 函数打印预测变量的重要性,并使用 h2o.varimp_plot()
函数仅针对领导者模型(automl 函数给出的最佳模型)绘制相同的图表。
h2o.varimp(aml@leader)
h2o.varimp_plot(aml@leader)
有没有办法打印排行榜中所有模型的预测变量的变量重要性,并使用上述两个函数绘制图表?
Stacked Ensembles(通常是领导者模型)尚不支持变量重要性 (JIRA here)。然而,其余模型的变量重要性可以在排行榜中的模型 id 循环中检索。请参阅下面的 R 代码。
library(h2o)
h2o.init()
# Import a sample binary outcome train/test set into H2O
train <- h2o.importFile("https://s3.amazonaws.com/erin-data/higgs/higgs_train_10k.csv")
# Identify predictors and response
y <- "response"
x <- setdiff(names(train), y)
# For binary classification, response should be a factor
train[,y] <- as.factor(train[,y])
# Run AutoML for 10 models
aml <- h2o.automl(x = x, y = y,
training_frame = train,
max_models = 10,
seed = 1)
# View the AutoML Leaderboard
lb <- aml@leaderboard
print(lb, n = nrow(lb))
# Get model ids for all models in the AutoML Leaderboard
model_ids <- as.data.frame(lb$model_id)[,1]
# View variable importance for all the models (besides Stacked Ensemble)
for (model_id in model_ids) {
print(model_id)
m <- h2o.getModel(model_id)
h2o.varimp(m)
h2o.varimp_plot(m)
}
我正在使用 R 中 H2o
包的 automl()
函数进行回归。
考虑一下我使用 "aml" 这个名字来构建模型。
aml <- h2o.automl(x=x, y=y, training_frame = train_set,
max_models = 20, seed = 1,
keep_cross_validation_predictions = TRUE)
automl()
的排行榜显示了表现最好的模型。我能够通过 h2o.varimp() 函数打印预测变量的重要性,并使用 h2o.varimp_plot()
函数仅针对领导者模型(automl 函数给出的最佳模型)绘制相同的图表。
h2o.varimp(aml@leader)
h2o.varimp_plot(aml@leader)
有没有办法打印排行榜中所有模型的预测变量的变量重要性,并使用上述两个函数绘制图表?
Stacked Ensembles(通常是领导者模型)尚不支持变量重要性 (JIRA here)。然而,其余模型的变量重要性可以在排行榜中的模型 id 循环中检索。请参阅下面的 R 代码。
library(h2o)
h2o.init()
# Import a sample binary outcome train/test set into H2O
train <- h2o.importFile("https://s3.amazonaws.com/erin-data/higgs/higgs_train_10k.csv")
# Identify predictors and response
y <- "response"
x <- setdiff(names(train), y)
# For binary classification, response should be a factor
train[,y] <- as.factor(train[,y])
# Run AutoML for 10 models
aml <- h2o.automl(x = x, y = y,
training_frame = train,
max_models = 10,
seed = 1)
# View the AutoML Leaderboard
lb <- aml@leaderboard
print(lb, n = nrow(lb))
# Get model ids for all models in the AutoML Leaderboard
model_ids <- as.data.frame(lb$model_id)[,1]
# View variable importance for all the models (besides Stacked Ensemble)
for (model_id in model_ids) {
print(model_id)
m <- h2o.getModel(model_id)
h2o.varimp(m)
h2o.varimp_plot(m)
}