问题描述
关于如何使用 survival
软件包在R中进行生存分析的信息(互联网和教科书)很多。但是,当您左审查数据时,我没有找到有关如何执行此操作的任何信息。
there is quite a lot of information (internet and textbooks) on how to do survival analysis in R with the survival
package. But I don't find any information on how to do this when you have left censored data.
问题背景:
我有一个带有公开生存数据的自建数据集。通常会给出事件的时间和最近一次随访的日期(正确的检查)。但是,有一项研究仅指出该事件发生在360天之前。所以我左删掉了此数据。
I have a self constructed data set with published survival data. Usually the event time and the date of the last follow-up (right censoring) is given. There is however one study that only states that the event happened before day 360. So I left censored this data.
想要做的事情:
我想分析带有左截断,事件和右截断的完整数据集。我想按性别绘制Kaplan-Meier曲线,然后
I want to analyse the complete data set with left truncation, events, and right truncation. I want to plot the Kaplan-Meier curve by gender and then
- 进行对数秩检验
- 进行Cox回归
我需要什么:
我能够使用 type = interval2
创建一个Surv对象。但是,这不允许计算生存包的 survdiff
或 coxph
。
I am able to create a Surv object with type = interval2
. But this does neither allow to calculate survdiff
, nor coxph
of the survival package.
intcox
软件包已从CRAN中删除,我在 icenReg $ c $中找不到搜索的内容c>或
间隔
软件包。
The intcox
package was removed from CRAN and I don't find what I search in the icenReg
or interval
packages.
有人可以给我吗后如何解决我的问题或在哪里可以找到有关此方面的实用信息?我已经在这上面花了几天时间。
Can anyone please give me a hind how to solve my problem or where to find practical information on this? I am already spending days on this one.
非常感谢!
推荐答案
通过使用 ic_sp
函数,您可以在 icenReg
中进行左右审查的Cox-PH模型。您可以使用标准的 Surv
响应变量来拟合它,即
You can fit a Cox-PH model with both right and left censoring in icenReg
by using the ic_sp
function. You can fit this using the standard Surv
response variable, i.e.
fit <- ic_sp(Surv(L, R, type = 'interval2') ~ treatment, data = myData)
或更简洁
fit <- ic_sp(cbind(L, R) ~ treatment, data = myData)
在 icenReg $ c中不支持日志等级测试$ c>,但可以在
interval
包中找到。
Log-rank tests are not available in icenReg
, but can be found in the interval
package.
这篇关于间隔检查数据:R的Cox比例风险和生存差异的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!