R 目前是否有一种方法可以使用生存软件包中的tt()函数来评估比例风险?

R 目前是否有一种方法可以使用生存软件包中的tt()函数来评估比例风险?,r,survival-analysis,cox-regression,survival,R,Survival Analysis,Cox Regression,Survival,我知道在生存包中,你可以定义一个函数,当你有违反比例风险的协变量时,它允许你转换时间 然而,在我看到的文档中,您应该能够使用cox.zph运行一个比例风险测试,但是每次我尝试运行它时,它都会说“函数不是为带有tt()项的模型定义的。” 这是最近发生的变化,还是我实施得不正确?代码如下 cox0 <- coxph(Surv(time=st.time, time2=end.time, event) ~ black + tt(black) + christian +

我知道在生存包中,你可以定义一个函数,当你有违反比例风险的协变量时,它允许你转换时间

然而,在我看到的文档中,您应该能够使用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))

cox0在“问题标记[r]”标题中写道“对于与统计相关的问题,请使用”。您可能会在那里得到更好的帮助。我还担心在公式中同时使用未转换的协变量和同一协变量的时间转换的概念。似乎这可能会产生严重的多重共线性问题。在“问题标记[r]”标题中,它说“对于与统计相关的问题,请使用”。您可能会在那里得到更好的帮助。我还担心在公式中同时使用未转换的协变量和同一协变量的时间转换的概念。似乎这可能会造成严重的多重共线性问题。