Данные, подвергнутые цензуре с интервалом: пропорциональная опасность Кокса и разница в суривале в R

Question

Есть довольно много информации (в Интернете и учебниках) о том, как провести анализ выживания в R с survival пакет. Но я не нахожу никакой информации о том, как это сделать, когда вы оставили цензурированные данные.

---

Проблема фона:

У меня есть собственный набор данных с опубликованными данными о выживании. Обычно указывается время события и дата последнего наблюдения (правая цензура). Однако есть одно исследование, в котором только утверждается, что событие произошло до дня 360. Поэтому я оставил цензуру этих данных.

Что я хочу сделать:

Я хочу проанализировать полный набор данных с усечением слева, событиями и усечением справа. Я хочу построить кривую Каплана-Мейера по полу, а затем

1. сделать тест лог-ранга
2. сделать регрессию Кокса

Что мне нужно:

Я могу создать объект Surv с type = interval2, Но это не позволяет ни рассчитать survdiffни coxph пакета выживания.

intcox пакет был удален из CRAN, и я не могу найти то, что я ищу в icenReg или же interval пакеты.

---

Кто-нибудь может подсказать, пожалуйста, как решить мою проблему или где найти практическую информацию по этому вопросу? Я уже трачу дни на это.

Большое спасибо!

Answer 1

Вы можете установить модель Cox-PH с правой и левой цензурой в icenReg используя ic_sp функция. Вы можете соответствовать этому, используя стандартные Surv переменная ответа, т.е.

fit <- ic_sp(Surv(L, R, type = 'interval2') ~ treatment, data = myData)

или немного более кратко с

fit <- ic_sp(cbind(L, R) ~ treatment, data = myData)

Тесты лог-ранга недоступны в icenReg, но можно найти в interval пакет.