R语言 Cox回归数据分层

12552

收藏 2016-03-29

各位大神，想问下Cox回归中设定了分层变量（level=2），也就是说会将数据按照不同的两层进行拟合，coxph()得出的其他变量的参数应该有两套啊，为什么只有一套呢？代码和输出结果如下：
>cox1<-coxph(Surv(days,status==1)~log(thick)+sex+strata(ulc))#这是模型，包含了一个分层变量ulc；
> summary(cox1) Call:coxph(formula = Surv(days, status == 1) ~ log(thick) + sex + strata(ulc)) n= 205, number of events= 57 coef exp(coef) se(coef) z Pr(>|z|) log(thick) 0.5599 1.7505 0.1784 3.139 0.0017 **sex 0.3600 1.4333 0.2702 1.332 0.1828 #sex和log(thick)的系数估计值只一个，是针对哪一层的呢---Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 exp(coef) exp(-coef) lower .95 upper .95log(thick) 1.750 0.5713 1.234 2.483sex 1.433 0.6977 0.844 2.434Concordance= 0.673 (se = 0.058 )Rsquare= 0.063 (max possible= 0.9 )Likelihood ratio test= 13.3 on 2 df, p=0.001296Wald test = 12.88 on 2 df, p=0.001598Score (logrank) test = 12.98 on 2 df, p=0.00152

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

seanfang1992

2016-6-23 20:10:08

我刚刚也做了，我觉得他是按ulc分层后，ulc各个曾经做的分析？求高手指点咯

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

xjh_109

2016-8-31 09:39:27

我也弄清楚了，确实是针对strata(ulc)中的ulc分类进行的，我修改了下，是这样进行的！
cox1<-coxph(Surv(days,status==1)~log(thick)+sex,data=mydata[ mydata$gender==1 ,])
这样的结果是针对gender赋值为1的分组进行HR分析！