區間刪失數據：Cox 比例風險和 R 中的生存差異

Question

有很多關於如何使用survival包在 R 中進行生存分析的信息（互聯網和教科書）。 但是，當您離開審查數據時，我沒有找到有關如何執行此操作的任何信息。

---

問題背景：

我有一個自建的數據集，其中包含已發布的生存數據。 通常會給出事件時間和最后一次隨訪（右刪失）的日期。 然而，有一項研究僅表明該事件發生在 360 天之前。所以我對這些數據進行了刪失。

我想做的事：

我想分析帶有左截斷、事件和右截斷的完整數據集。 我想按性別繪制 Kaplan-Meier 曲線，然后

1. 做一個對數秩檢驗
2. 做一個 Cox 回歸

我需要的：

我能夠創建一個type = interval2的 Surv 對象。 但這既不允許計算survdiff ，也不允許計算生存包的coxph 。

intcox包已從 CRAN 中刪除，我在icenReg或interval包中找不到我搜索的icenReg 。

---

任何人都可以告訴我如何解決我的問題或在哪里可以找到有關此的實用信息？ 我已經在這個上花了幾天時間。

非常感謝！

Answer 1

您可以使用ic_sp函數在icenReg擬合具有左右icenReg的 Cox-PH 模型。 您可以使用標准的Surv響應變量進行擬合，即

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

或者更簡潔一點

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

對數秩測試在icenReg中不可用，但可以在interval包中找到。