関心のある 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回帰を行っています（C統計計算には「デザイン」パッケージのcphも使用しています）。
私の質問：
この左切り捨てと右打ち切りを実行するより良い方法はありますか?
理想は次のとおりです。
# TREOP - time in years
# RREOP - event
surv <- Surv(TREOP, RREOP, start=30/365.25, stop=.5)
ヘルプと time、time2 & type は切り捨てを処理しているようですが、被験者が 22 日後に試験に参加するより複雑な設定のためであり、データを間隔に分割するためのものではないと思います。
編集
SurvSplit() 関数をサバイバル パッケージで見つけましたが、説明からすると正しいように見えますが、それを飼いならす方法がわかりません。この例は実際には役に立ちません。誰でも経験がありますか？