使用朝阳视图绘制 rpart 决策树模型

Question

我找到了一种使用 sunburstR 包从 rpart 绘制决策树解决方案的方法。要绘制旭日形图，必须有 data.frame 表示一个序列。我将决策树结果修改为如下序列

决策树的结果

rpart(Species~.,data=iris)

n= 150 

node), split, n, loss, yval, (yprob)
      * denotes terminal node

1) root 150 100 setosa (0.33333333 0.33333333 0.33333333)  
  2) Petal.Length< 2.45 50   0 setosa (1.00000000 0.00000000 0.00000000) *
  3) Petal.Length>=2.45 100  50 versicolor (0.00000000 0.50000000 0.50000000)  
    6) Petal.Width< 1.75 54   5 versicolor (0.00000000 0.90740741 0.09259259) *
    7) Petal.Width>=1.75 46   1 virginica (0.00000000 0.02173913 0.97826087) *

旭日序列：

sequences_1<-1
sequences_1<-data.frame(sequences_1)
sequences_1[1:3,]<-1
sequences_1$V1[1]<-"Petal.Length<_2.45-setosa"
sequences_1$V1[2]<-"Petal.Length>=2.45-Petal.Width<_1.75_54_5-versicolor"
sequences_1$V1[3]<-"Petal.Length>=2.45-Petal.Width>=1.75_46_1-virginica"
sequences_1$V2[1]<-50
sequences_1$V2[2]<-54
sequences_1$V2[3]<-46
sequences_1$sequences_1<-NULL

绘制旭日图：

library(sunburstR)
sunburst(sequences_1,count=TRUE)

旭日图的序列，我手动做的。有人知道如何从 rpart 决策树的结果中像上面那样自动构建序列吗？

Answer 1

d3r提供函数d3_party将rpart/partykit转换成d3层级。 sunburst 可以使用 d3_party 的结果，稍加修改即可将 "rule" 更改为 "name"。这并不理想，但在大多数情况下可以完美运行。

library(rpart)
library(d3r)
# d3_party requires partykit
# install.packages("partykit")
library(sunburstR)

rp <- rpart(Species~.,data=iris)
rp_d3 <- d3_party(rp)

# one trick/hack required since sunburst expects
#   name but d3_party gives rule
#   this is ugly but let's replace all "rule" with "name"
#   with gsub
rp_d3 <- gsub(
  x = rp_d3,
  pattern = '"rule":',
  replacement = '"name":'
)

sunburst(
  rp_d3,
  valueField = "n",
  sumNodes = FALSE,
  count = TRUE,
  legend = FALSE
)

使用朝阳视图绘制 rpart 决策树模型

Plotting a rpart decision tree modell with the sunburst view

r

sunburst-diagram