如何获取ctree创建的生存对象？

我想得到由ctree函数创建的生存对象。 原因是我想得到描述树的每片叶子中曲线的向量

有人知道怎么做吗

提前谢谢！

---

Alley

您可以使用

treeresponse

方法获得曲线。也许有更好的方法，但这就是我想到的

以下是使用

？ctree

中的生存树示例的图示：

require(party)
data("GBSG2", package = "ipred")
GBSG2ct <- ctree(Surv(time, cens) ~ .,data = GBSG2)
plot(GBSG2ct)

每个组件out

out

都是生存对象，类

“survfit”

对于这棵树，我们有四个终端节点，因此只有四个唯一的生存对象。您可以使用

where

方法查看观测值所在的节点

wnode <- where(GBSG2ct)

其中，除了轴上的范围之外，是先前绘制的树上节点3的面板中使用的绘图

---

对于本例，您可以对节点4、6和7重复此操作，但您需要根据数据和拟合树调整使用的节点。

ctree

不是基本的R函数。你指的是什么包裹？一些示例代码也将帮助您获得答案。@Prasad，它在package

party

中。加载包后，请参见

？ctree

，以获得一个简单的生存树示例。非常感谢Parsad！我试试看。你知道我怎么用生存物体的“surv”参数吗？我想得到生成曲线的概率向量。再次感谢-Alley@Alley您需要查看生存软件包的帮助，并阅读

？survfit.object

谢谢Gavin，我会查看的。您帮了大忙，很抱歉上次没认出是您。好的，使用：out[[1]]$surv确实有效！：）顺便问一下，你为什么要写：out[[1]]而不是out[1]？这是什么意思？再次谢谢你去学点什么

[

当应用于这样的列表时，

out[1]

返回一个包含单个组件的列表，即所选组件。

[

将返回该列表中的对象。

> class(out[[1]])
[1] "survfit"

wnode <- where(GBSG2ct)

> n3 <- which(wnode == 3 & !duplicated(wnode))
> n3
[1] 1
> out[[n3]]
> out[[n3]]
records   n.max n.start  events  median 0.95LCL 0.95UCL 
    686     248     248      88    2093    1814      NA

plot(out[[n3]], conf.int = FALSE, mark.time = FALSE)