首页 > 编程知识 正文

r语言经典实例,coxa

时间:2023-05-06 15:51:28 阅读:176233 作者:3810

与部分共变量不满足PHA的情况相比,extendedcoxmodelextendedcoxmodel通过在模型中添加它们的共变量和时间的interaction项,可以考虑共变量随时间变化的情况。

接下来考虑的是sex不满足PHA的情况。

library(survival ) Anderson=read.table (choose.files )、col.names=c )、' time '、' censor '、' sx、' logWbBBBI cut=8,end='time ',event='censor ', start='start ' ) #以t=8做cut,分割后的数据Anderson2$ gt=(Anderson2$ start==8)0# heavisidefunctiondetach (Anderson ) attach ) Anderson2) s=surv ) staaach censor (PHM _2=Cox ph (s ) sxtrtlogwbcsx :

call 3360 Cox ph (formula=s至sxtrtlogwbcsx : gt ) n=66,numberofevents=30coefexp(coef ) se ) coef ) zpr(|z| ) sx0.9813 )2. 6679.5412.813.0698.TRT1. 1583.1847 ) 450.0155 log WBC1. 5834.87130.34274.6203.83 e-06 * * * lower.95 upper.95sx2. 6679.37480.92367.7059 TRT3. 18470.31401.26292.0307 log WBC4. 8713.2053.488539.535356 ssssssss300000000 Likelihood ratio test=52.16 on 4 df、p=1e-10Wald test=38.25 on 4 df、p=1e-07score(logrank ) test=57.08 on 4 df、 p=1e-11 stratifiedcoxmodelstratifiedcoxmodel也针对不满足PHA而提出,这里的思想是对不满足PHA的合作变量层次化

协变量的效果在一层(部分)中是相同的。 也就是说,层内没有interaction,效果是常数。 这就是非交互辅助。

对于“no interaction”模型,虽然每层的基本函数不同,但指数项系数一致;

考虑interaction的模型baseline function h 0 h_0 h0和系数在各层不同。

phm_si=Coxph(s~# ) strata(sx ) trtlogwbcsx : logwbcsx : TRT ) strata函数对sx层次summary (PHM _ si ) call3360coxph ) FHM_si ) Cox ph trtlogwbcsx numberofevents=30coefexp(coef ) se ) coef ) zpr(|z| ) TRT0.3113 ) 1.3652 ) 0.5636.552(5807logWBC1.2061 ) 3.3406 ) 5035.396.0000 log WBC : sx 0.53661.71020.73520.73520 ign if.codes :‘* * * 0.001’0.05‘.’0.1‘’1 exp (coef ) exp (coef )-coef ) ) )。 lower.95 upper.95 TRT 1.365.73250.4524.120 log WBC 3.3410.29931.2452.962 log WBC 3360 sx 1.710.5847.4048.226 ) Likelihood ratio test=35.83 on 4 df,p=3e-07Wald test=21.69 on 4 df,p=2e-04scorrdance ) )。

版权声明:该文观点仅代表作者本人。处理文章:请发送邮件至 三1五14八八95#扣扣.com 举报,一经查实,本站将立刻删除。