处理违反比例风险假设的变量,可以同时使用阶跃函数和时间变换函数吗?

Dealing with variables violating proportional hazard assumption, is it okay to use both step function and time-transform function simultaneously?

提问人:JY J 提问时间:8/2/2023 最后编辑:JY J 更新时间:8/23/2023 访问量:30

问:

我尝试用 4 个违反比例风险假设的变量拟合生存模型。我阅读了以下小插曲(https://cran.r-project.org/web/packages/survival/vignettes/timedep.pdf)并在网上彻底烙印,但找不到类似的案例。.

是否可以对违反比例风险假设的变量使用阶跃函数? 例如

vet2 <- survSplit(Surv(time, status) ~ ., data= veteran, cut=c(90, 180),
episode= "tgroup", id="id")
fit <- coxph(Surv(tstart, time, status) ~ var1:strata(tgroup) + var2:strata(tgroup) + var3:strata(tgroup) + var4 + var5 + var6, data=vet) 

还是同时使用阶跃函数和时间变换函数? 如

fit1<-coxph(Surv(tstart, time, status) ~ var1:strata(timegroup) + var2:strata(timegroup)+ tt(var3) + var4 + var5 + var6, 
                                         data=vet,
                                         tt=function(x, t, ...) x*log(t+20))

期待您的帮助!!

R 生存分析 回归 生存 Cox

评论

0赞 Mike 8/2/2023
这个问题更适合统计交换,因为它是一个统计问题,而不是一个编程问题
0赞 IRTFM 8/14/2023
我在这两个公式中都没有看到“阶跃函数”。
0赞 JY J 8/23/2023
感谢您的评论。我会问这个问题来统计交换。
0赞 JY J 8/23/2023
@IRTFM,我只是使用了下面的小插曲部分 4.1 Step 函数 (cran.r-project.org/web/packages/survival/vignettes/timedep.pdf) 中的示例,尽管我在开始时错过了 survSplit() 部分。我修改了代码。Pleast,如果还有其他不清楚的地方,请告诉我。

答: 暂无答案