间隔检查数据:R的Cox比例风险和生存差异 [英] Interval censored data: Cox proportional hazard and surival difference in R

查看:161
本文介绍了间隔检查数据:R的Cox比例风险和生存差异的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!

问题描述

关于如何使用 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


  1. 进行对数秩检验

  2. 进行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 间隔软件包。

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 ,但可以在 interval 包中找到。

Log-rank tests are not available in icenReg, but can be found in the interval package.

这篇关于间隔检查数据:R的Cox比例风险和生存差异的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!

查看全文
登录 关闭
扫码关注1秒登录
发送“验证码”获取 | 15天全站免登陆