Создайте объект выживания с помощью наблюдений как правой, так и левой усеченной + правой.

Question

У меня проблема с созданием моих объектов выживания в R. Я хочу смоделировать выживание от клиентов по подписке (ежемесячные данные), но для создания объектов выживания мне нужно включить оба типа цензуры:

• классическая правая цензура: некоторые клиенты вошли в периоды наблюдения и еще не отписались ("событие смерти")

• усечение влево и цензура вправо: другие клиенты вводили ДО периодов наблюдения, но неизвестно, когда, потому что не было таблиц хронологии из ранее

Конечно, у меня есть случаи событий (отписались). Но проблема для меня заключается в том, как создать объекты выживания обоих типов в одном наборе данных для последующего моделирования. Я думаю, что если не принимать во внимание "усеченную" ситуацию, это смещет и переоценит длительность пребывания некоторых клиентов, поэтому я не хочу отказываться от этих дел.

Таким образом, я знаю время начала для тех, кто вступил в период наблюдения. Но для тех, кто вошел ранее, у меня просто время начала "0", а не реальные (неизвестные).

Пока я пробовал эти коды:

1) survobj <- Surv(TIME, EVENT)
## i loose information of the truncated ones.

2) survobj <- Surv(ifelse(T0==0,NA,T0), T1, EVENT) 
## will create "interval-censored" objects, not my case

3) survobj <- Surv(T0, data$T1, EVENT, type='counting')  
## all objects will be "left truncated", not just the T0=0

Заранее благодарю за любую помощь.

Answer 1

Одним из способов включения как левой, так и правой цензуры в длительность является создание Surv объект с type = interval2, Как survival Документация гласит:

[...] думайте о каждом наблюдении как о временном интервале с (-infinity, t) для левой цензуры, (t, бесконечность) для правой цензуры, (t, t) для точного и (t1, t2) для интервала. Этот подход используется для type = interval2. Бесконечные значения могут быть представлены либо фактической бесконечностью (Inf), либо NA.

Следовательно, вам нужно назначить известные длительности в векторе для time (где будут оставлены случаи цензуры NA) и в другом для time2 (где будут правильные случаи цензуры NA).

Подобная проблема обсуждается здесь.