R语言 Cox回归 数据分层

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位大神，想问下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

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：Cox R语言 concordance Likelihood possible number events status

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

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

请问，楼主解决了嘛，我现在也有同样的疑惑

立妍Z 发表于 2021-5-24 10:56
请问，楼主解决了嘛，我现在也有同样的疑惑

可以看下，分层cox回归的定义，里面有个假设，是假设所有层的系数(即beta)的值是相等的，因此，这里算出来的系数并不是看哪一层的。可以说是，所有层的系数是相等的。这里分层cox回归主要是假设，不同层之间的基。准危险函数是不同的（即违背了风险成比例这一个假设）。希望可以对你们有所帮助、