Я хотел бы проверить все наблюдения, которые достигли некоторого узла в дереве решений rpart. Например, в следующем коде:

fit <- rpart(Kyphosis ~ Age + Start, data = kyphosis)
fit

n= 81 

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

 1) root 81 17 absent (0.79012346 0.20987654)  
   2) Start>=8.5 62  6 absent (0.90322581 0.09677419)  
     4) Start>=14.5 29  0 absent (1.00000000 0.00000000) *
     5) Start< 14.5 33  6 absent (0.81818182 0.18181818)  
      10) Age< 55 12  0 absent (1.00000000 0.00000000) *
      11) Age>=55 21  6 absent (0.71428571 0.28571429)  
        22) Age>=111 14  2 absent (0.85714286 0.14285714) *
        23) Age< 111 7  3 present (0.42857143 0.57142857) *
   3) Start< 8.5 19  8 present (0.42105263 0.57894737) *

Я хотел бы видеть все наблюдения в узле (5) (т.е. 33 наблюдения, для которых Start> = 8.5 и Start <14.5). Очевидно, я мог бы вручную добраться до них. Но я хотел бы иметь такую ​​функцию, как (скажем) "get_node_date". Для которого я мог бы просто запустить get_node_date (5) - и получить соответствующие наблюдения.

Любые предложения о том, как это сделать?