Cox回歸中分層變量的參數估計和方差（分層/生存包）

Question

我使用survival計划運行Cox回歸來計算暴露A的死亡率風險比。我發現age變量違反比例風險假設（使用cox.zph ）並使用strata(age)來對其他模型中的strata(age)進行分層。

我需要一個age變量的參數估計，以及方差和協方差矩陣（計算速率提升期）......我不知道在哪里可以找到它們！

我錯過了什么，或者我誤解了什么strata在做什么？

---

這是一個可重復的例子，使用來自survival包的lung數據。

library(survival)    

我創建了生存對象並使用非分層年齡變量進行第一次Cox回歸。

lung$SurvObj <- with(lung, Surv(time, status == 2))
coxreg1 <- coxph(SurvObj ~ age + sex, data =  lung)  

因此，我得到參數估計的系數，方差和協方差矩陣。

>   coxreg1$coefficients
        age         sex 
 0.01704533 -0.51321852 

> vcov(coxreg1)
             age          sex
age 8.506877e-05 8.510634e-05
sex 8.510634e-05 2.804217e-02

現在，如果使用分層年齡變量進行第二次回歸，我不會得到任何系數估計，方差或協方差。

coxreg2 <- coxph(SurvObj ~ strata(age) + sex, data =  lung)

> coxreg2$coefficients
     sex 
-0.64471 

> vcov(coxreg2)
          sex
sex 0.0449369

謝謝您的幫助！

Answer 1

當您使用變量進行分層時，您不會得到任何系數估計值。 相反，估計不同年齡組的單獨基線危害。 分層Cox回歸的本質是擬合在每個層中具有不同基線危險的模型。