partykitからctreeを使用して、各ターミナルノードのリーフからルートへのフルパスを取得する

Question

現在、Rパッケージ「partykit」からctreeを使用していますが、ターミナルノードからルートへのフルパスを取得する方法があるかどうか疑問に思っていました。各リーフについて、ノードIDを含むベクトルとして表されるルートへのフルパスが必要です。

library(partykit)
ct <- ctree(Species ~ ., data = iris) 

Model formula:
Species ~ Sepal.Length + Sepal.Width + Petal.Length + Petal.Width

Fitted party:
[1] root
|   [2] Petal.Length <= 1.9: setosa (n = 50, err = 0.0%)
|   [3] Petal.Length > 1.9
|   |   [4] Petal.Width <= 1.7
|   |   |   [5] Petal.Length <= 4.8: versicolor (n = 46, err = 2.2%)
|   |   |   [6] Petal.Length > 4.8: versicolor (n = 8, err = 50.0%)
|   |   [7] Petal.Width > 1.7: virginica (n = 46, err = 2.2%)

Number of inner nodes:    3
Number of terminal nodes: 4

木のプロット

これは基本的に私が必要なものです：

[[1]]
[1] 2 1

[[2]]
[1] 5 4 3 1

[[3]]
[1] 6 4 3 1

[[4]]
[1] 7 3 1

助けていただければ幸いです！ありがとう！

Answer 1

次の関数でうまくいくはずです。最初の行はノードごとの子のリストを抽出し、これからすべてのノードを再帰的に調べることができます。

get_path <- function(object) {
  ## list of kids per node (NULL if terminal)
  kids <- lapply(as.list(object$node), "[[", "kids")

  ## recursively add node IDs of children
  add_ids <- function(x) {
    ki <- kids[[x[1L]]]
    if(is.null(ki)) {
      return(list(x))
    } else {
      x <- lapply(ki, "c", x)
      return(do.call("c", lapply(x, add_ids)))
    }
  }
  add_ids(1L)
}

これは、任意のpartyオブジェクトに適用できます。

get_path(ct)
## [[1]]
## [1] 2 1
## 
## [[2]]
## [1] 5 4 3 1
## 
## [[3]]
## [1] 6 4 3 1
## 
## [[4]]
## [1] 7 3 1