[英]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.