Выполнение одностадийной регрессии Кокса PH в течение 4 временных интервалов в R

Question

У меня 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
....

Я делаю регрессию Кокса с пакетом «выживания» (я также использую cph в пакете «Дизайн» для вычислений C-статистики).

Мой вопрос:

Есть ли лучший способ выполнить левое усечение и правую цензуру?

Идеально будет:

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

Я посмотрел на справку, и время, время2 и тип, кажется, справляются с усечением, но я думаю, что это более сложная обстановка, когда субъекты входят в исследование через 22 дня, а не для разделения данных на интервалы.

Редактировать

Я нашел функцию SurvSplit () в пакете для выживания, но, хотя по описанию он кажется правильным, я не уверен, как его приручить - пример на самом деле мне не помогает. У кого-нибудь есть опыт?

Answer 1

1. Я согласен с цензурой, которая выглядит просто и понятно.
2. Я не уверен, что вы должны урезать влево. Я бы чувствовал себя более комфортно, оставляя короткие времена выживания неизменными и просто увеличивая верхний предел цензуры. Если n-й период времени намного длиннее (n-1) -го - это не будет иметь большого значения, и если он не намного дольше, чем более короткое время выживания, то его не следует усекать.