[英]Doing a one step Cox PH regression for 4 time intervals in R
我有4个兴趣区间:
现在,我像这样子集数据集:
# 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()函数,但尽管从描述上看它似乎是正确的,但我不确定如何驯服它-该示例并不能真正帮助我。 任何人都有经验吗?
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.