繁体   English   中英

在R中以4个时间间隔进行一步Cox PH回归

[英]Doing a one step Cox PH regression for 4 time intervals in R

我有4个兴趣区间:

  • 0至30天
  • 30天-½年
  • ½-2年
  • 2年-10年

现在,我像这样子集数据集:

# Set time period
time_period.first <-  30/365.25
time_period.intermediate <- .5
....

# TREOP = Time in years
data.first = all_data
# Remove already censored data
data.intermediate = subset(data.first, data.first$TREOP > time_period.first)

# Set all outside as censored
data.first$RREOP[data.first$TREOP > time_period.first] = 0
data.first$TREOP[data.first$TREOP > time_period.first] = time_period.first

data.intermediate$RREOP[data.intermediate$TREOP > time_period.second] = 0
data.intermediate$TREOP[data.intermediate$TREOP > time_period.second] = time_period.second
....

我正在使用“生存”程序包进行cox回归(我也使用“设计”程序包中的cph进行C统计量计算)。

我的问题:

有没有更好的方法来执行这种左截断和右删截?

理想的是:

# TREOP - time in years
# RREOP - event
surv <- Surv(TREOP, RREOP, start=30/365.25, stop=.5)

我已经查看了帮助,时间,time2和类型似乎可以处理截断,但我认为这是为了让受试者在22天后进入研究环境,而不是将数据分成多个间隔。

编辑

我已经在生存程序包中找到了survSplit()函数,但尽管从描述上看它似乎是正确的,但我不确定如何驯服它-该示例并不能真正帮助我。 任何人都有经验吗?

  1. 我同意看起来简单明了的权利审查。
  2. 我不确定您是否应该截断。 我会更舒适地保留较短的生存时间,而只是增加检查上限。 如果第n个时间段比第(n-1)个时间段长得多-无关紧要,并且如果它不比较短的生存时间长很多,则不应缩短。

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM