中位生存率

Question

我試圖用R中的生存軟件包對中位生存率進行對數秩檢驗 - 這可能嗎？我正在使用的代碼：

survdiff(surv(Age, Survival) ~ Treatment, data = GBM) 

其中年齡=診斷時的年齡;生存=他們倖存多少年;治療=是，他們是否接受過治療（編碼爲1或2）

如果這是正確的，爲什麼我收到錯誤：

Error in survdiff.fit(y, groups, strata.keep, rho) : 
    There is only 1 group 
In addition: Warning message: 
In Surv(GBM$Age, GBM$Survival) : Invalid status value, converted to NA 

不治療需要編碼爲0或1，而不是1或2？

Answer 1

pkg中沒有surv函數：生存。這是Surv。

?Surv 

協變量是編碼爲0/1還是1/2都無關緊要。估計值將描述依賴於一個單位差異的風險比率。如果您想要設置Surv電話來處理年齡和持續時間，那麼活動時間應該是年齡+生存。當然，你需要一個事件指示器。