Я пытаюсь подобрать модель с интервальной цензурой и усечением влево с помощью пакета R survival.
Контекст
В эпидемиологическом анализе время начала может не иметь значения (заполнение формы не является событием само по себе), поэтому люди могут быть включены в исследование в соответствии с их возрастом.Тогда ось x больше не time, а age.Поскольку люди сообщают о времени наблюдения только с их возраста при поступлении, данные тогда усекаются влево (у нас нет информации об их характеристиках до их возраста при поступлении).Вы можете найти больше об этом здесь .
Мои события собираются с помощью анкет, поэтому событие может происходить где угодно между анкетами n и n-1.Таким образом, они подвергаются цензуре с интервалом .Распространенным способом решения этой проблемы является рассмотрение события, прошедшего цензуру, в середине периода между анкетами n и n-1, но мне интересно, может ли быть лучший способ.
Код
AFAIK, можно провести анализ выживаемости по усеченной справа (по возрасту) цензуре, используя структуру:
db$event = ifelse(is.event, 0, 1)
Surv(age_at_inclusion, age_at_event_or_censoring, event, type="counting")
#here "counting" could be omitted as it is the default when time2 is used
Здесь используется та же логика, что и в справочном примере Surv с продольнымнабор данных heart: Surv(heart$start, heart$stop, heart$event), но в поперечном наборе данных.
Но, читая справку Surv() , кажется, что вы можете добиться интервальной цензуры, используя следующую структуру:
db$event = ifelse(is.event, 0, 3)
Surv(time_interval_left, time_interval_right, event, type="interval")
Поскольку они, кажется, используютта же структура, как я могу смоделировать событие как с левой, так и с правой цензурой?