如何在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-80y，而不是40-50yo怎么办？

我如何告诉Coxph将N'th群集作为默认值？

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)