如何在R中的Cox模型中選擇默認群集？

Question

我正在研究Cox模型，並且有一個在群集中剪切的變量。

這是一些工作示例：

lung2=lung
lung2$age.cl = cut(lung2$age, c(40,50,60,70,80,90))
cox = coxph(Surv(time, status) ~ age.cl, data =  lung2)

它返回此表：

                coef exp(coef) se(coef)    z     p
age.cl(50,60]  0.184     1.202    0.297 0.62 0.536
age.cl(60,70]  0.237     1.267    0.288 0.82 0.411
age.cl(70,80]  0.532     1.703    0.309 1.73 0.085
age.cl(80,90]  3.926    50.707    0.825 4.76 2e-06

默認情況下， coxph使用第一個集群作為默認值來計算所有內容。 但是，如果臨床默認值為70-80歲而不是40-50歲怎么辦？

如何告訴coxph將第n個簇設為默認簇？

Answer 1

lung2=lung
lung2$age.cl = cut(lung2$age, c(40,50,60,70,80,90))
lung2$age.cl_relevel <- relevel(lung2$age.cl, '(70,80]')
cox = coxph(Surv(time, status) ~ age.cl_relevel, data =  lung2)