簡體   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