从tidymodel绘制决策树结果

问题描述

我已经成功地使用tidymodels软件包构建了决策树模型，但是我不确定如何提取结果并绘制树形图。我知道我可以使用rpart和rpart.plot软件包来实现相同的目的，但是我宁愿使用tidymodels，因为这是我正在学习的内容。以下是使用mtcars数据的示例。

library(tidymodels)
library(rpart)
library(rpart.plot)
library(dplyr) #contains mtcars

#data
df <- mtcars %>%
    mutate(gear = factor(gear))


#train/test
set.seed(1234)

df_split <- initial_split(df)
df_train <- training(df_split)
df_test <- testing(df_split)


df_recipe <- recipe(gear~ .,data = df) %>%
  step_normalize(all_numeric())


#building model
tree <- decision_tree() %>%
   set_engine("rpart") %>%
   set_mode("classification")

#workflow
 tree_wf <- workflow() %>%
   add_recipe(df_recipe) %>%
   add_model(tree) %>%
   fit(df_train) #results are found here 

rpart.plot(tree_wf$fit$fit) #error is here

我得到的错误是Error in rpart.plot(tree_wf$fit$fit) : Not an rpart object，这是有道理的，但我不知道是否缺少将结果转换为rpart.plot可以打印的格式的软件包或步骤。这可能是不可能的，但将不胜感激。

解决方法

您还可以使用workflows::pull_workflow_fit()功能。它使代码更加优雅。

tree_fit <- tree_wf %>% 
  pull_workflow_fit()
rpart.plot(tree_fit$fit)

以下作品（请注意额外的$fit）：

rpart.plot(tree_wf$fit$fit$fit)

这不是一个非常优雅的解决方案，但是它确实绘制了树。

使用防风草0.1.3和rpart.plot 3.0.8进行了测试。

decision-tree r r tidymodels