如果我有一个cox比例风险模型,作为事件的预测因子一个分类变量(和其他协变量,如年龄等),让我们说例如我感兴趣的分类变量是肿瘤大小,可以是0-10例如,10-20,20-30,我看到随着肿瘤大小的增加,存在死亡率更高的趋势,我如何计算r并得到ap?
答案 0 :(得分:2)
接受不完整的规格存在一些危险,因为你可能会有一个大小分类会弄乱它,例如:0-10,10-20,20-30,30-80,80-120 ,120-240。 除非你仔细构建了你的因子才能正确地按顺序升序,否则我将建议的捷径将会失败。
survmdl <- coxph(Surv(tie,event) ~ as.numeric(fact), data=dat)
你会得到一个“趋势测试”,它将被解释为对于上升大小的日志(危险)的每类别增加,它将是单个系数。如果您想要更好的考虑回复,请发布您的实际因素水平。