r - 目前有没有办法使用生存包中的 tt() 函数来评估比例风险？

Question

我知道在生存包中，您可以定义一个函数，如果您有一个违反比例风险的协变量，您可以在该函数中转换时间。

但是，在我看到的文档中，您应该能够使用 cox.zph 运行比例风险测试，但是每次我尝试运行它时，它都会说“没有为具有 tt() 术语的模型定义函数。”

这是最近发生的变化，还是我执行不正确？代码如下。

cox0 <- coxph(Surv(time=st.time, time2=end.time, event) ~ black + tt(black) + christian + 
                cluster(id), data =df,
              tt = function(x, t, ...) {
                mtrx <- model.matrix(~x)[,-1]
                mtrx * log(t)
              })

test <- cox.zph(cox0, 
                transform=function(t)log(t))

标签： rsurvival-analysiscox-regressionsurvival