cox比例风险回归模型及其R程序课件.ppt

1、CoxCox回归分析回归分析刘瑞红1ppt课件 )exp(0000expexp,exp,XtttSdtXthdtXthXtSCox 模型不直接考察生存函数 tS与协变量的关系，而 是用风险率函数 th作为因变量，并假定：mmXXXexpthXexpthX,th221100 利用生存率函数S(t,X)与风险函数h(t,X)的关系可导出 较好地解决截尾值的问题 反映了协变量X与生存函数的关系 Cox模型的基本形式 2ppt课件所有危险因素为0时的基础风险率，它是未知的，但假定它与h(t,X)是呈比例的。右侧可分为两部分：h0(t)没有明确的定义，分布无明确的假定，参数无法估计，为非参数部分；另一部

2、分是参数部分，其参数可以通过样本的实际观察值来估计的，正因为Cox模型有非参数和参数两部分组成，故又称为半参数模型。3ppt课件 若若 i0，则，则RR1，该因素为危险因素；，该因素为危险因素；若若 i0，则，则RR0，说明该个体死亡风险高于平均水平；l若PIj0，说明该个体死亡风险低于平均水平。mmjjjxbxbxbPI12211ibiiiisbb ijx27ppt课件3通过估计生存率，对群体定量地进行预后评价。l由于生存率与基础生存率相关，故只要估计出基础生存率，再结合各因素的偏回归系数就可以估计出生存率，即 XtttSdtXthdtXthXtS expexpexp,exp,000028p

3、pt课件29ppt课件ti 时刻的基础时刻的基础生存率公式生存率公式 ti 时刻的基时刻的基础累计风险础累计风险函 数 公 式函 数 公 式在在tk时刻的死亡人数时刻的死亡人数 ti 时刻的生存时刻的生存率 计 算 公 式率 计 算 公 式 2.Breslow法30ppt课件5.Cox模型拟合优度的考察l可将研究对象按个体预后指数恰当地分组，用乘积极限法估计各组的生存率曲线，并与按Cox模型预后指数分类的生存率曲线在同一坐标系内进行比较，若两种曲线具有一致性，说明模型拟合较理想。也可以使用卡方检验的方法。31ppt课件比例风险假定（比例风险假定（PH假定）：假定）：lPH假定的检验方法大致可以

4、分为图法和正规的检验方法图法和正规的检验方法。l图法：即通过观察散点图中散点的分布或趋势是否满足既定模型基本假设下的形状来判断资料是否满足或近似满足模型假定，主要的方法有：比较COX-KM生存曲线、基于累计风险函数的图示法、Schoenfeld残差图、Score残差图；l正规的检验方法：通过构造满足既定模型基本假定下服从某一已知分布的统计量，利用P值来检验资料是否满足或近似满足模型假定。主要方法有：时协变量法、线性相关检验、加权残差Score检验、三次样条函数法l.Cox比例风险回归模型诊断及预测有关问题的研究_余红梅.caj.caj32ppt课件对数线性假定对数线性假定l主要的方法：多重法；

5、鞅残差法；改进的鞅残差法；构造变量法Cox回归影响点的分析鞅残差，剩余残差，得分残差，Schoenfeld残差，加权Schoenfeld残差33ppt课件R程序：程序：l在survival包中的函数：lcoxph(formula,data,weights,subset,na.action）l formula:反应变量要有Surv()函数定义lcox.zph(fit,transform=,global=T/F)检验比例风险假定lbasehaz(fit,centered=T/F)计算基础风险率的函数34ppt课件getwd()setwd(C:/Users/Administrator/Desktop

6、)my-read.csv(pharynx.csv,header=T)attach(my)my-myCOND!=9&GRADE!=9,#去掉两个缺失值 CONDCOND=3|COND=4-2CONDCOND=0-1aa-table(COND)aamydata-my,-c(1,11)l library(survival)lfit-coxph(Surv(TIME,as.numeric(STATUS).,data=mydata)l#cox比例风险模型函数lplot(survfit(fit)l summary(fit)35ppt课件lbasehaz(fit)#基础风险函数laa-cox.zph(fit)

7、#比例风险假定的检验lprint(aa)lplot(aa6)36ppt课件多重分数多项式模型多重分数多项式模型llibrary(mfp)lf-mfp(Surv(TIME,as.numeric(STATUS)fp(AGE,df=4,select=0.05)+INST+SEX+TX+GRADE+COND+SITE+T_STAGE+N_STAGE,select=0.15,family=cox,data=mydata)lprint(f)l（rsq=1-sum(f$residuals)2/sum(mydata$TIME-mean(mydata$TIME)2)#R237ppt课件SAS程序：程序：lpro

8、c phreg data=;lmodel=/选项；lstrata lfreq；lby；lProportionality_test:test;/*比例风险假定,时协变量法*/lrun;38ppt课件llibname ll F:R语言学习COX;ldata a;lset ll.pharynx;lrun;ldata a;lset a;lif cond=9 or grade=9 then delete;lif cond=3 or cond=4 then cond=2;lif cond=0 then cond=1;lrun;lproc phreg data=a;lmodel time*status(0)

9、=tx cond site t_stage n_stage age grade inst/selection=stepwise risklimits sle=0.15 sls=0.15;lbaseline out=cox survival=suv stderr=stderr /method=ch;lrun;39ppt课件四、Cox模型的适用范围及注意事项1.适用范围适用范围lCox模型适用于生存资料的统计分析，属半参数模型，对模型适用于生存资料的统计分析，属半参数模型，对资料没有特殊的要求，也可以估计各因素的参数，并能资料没有特殊的要求，也可以估计各因素的参数，并能做多因素的统计分析。该模型的

10、主要优点在于能从众多做多因素的统计分析。该模型的主要优点在于能从众多的影响因素中排除混杂因素的影响，找出影响生存时间的影响因素中排除混杂因素的影响，找出影响生存时间的因素，根据各因素的参数估计出个体的生存率。的因素，根据各因素的参数估计出个体的生存率。l另外，另外，Cox模型能分析具有截尾数值的生存时间。模型能分析具有截尾数值的生存时间。lCox模型可以分析多种因素对疾病预后的影响，使其更适模型可以分析多种因素对疾病预后的影响，使其更适合于临床的随访研究。合于临床的随访研究。Cox模型使临床观察的定性指标又模型使临床观察的定性指标又加上定量指标进行分析，提高了分析的效率。加上定量指标进行分析，

11、提高了分析的效率。40ppt课件2.注意事项注意事项(1)设计阶段应注意的问题设计阶段应注意的问题 样本的代表性要好。样本的代表性要好。样本的含量为观察协变量的5-20倍。所有危险因素要在设计时考虑全面，避免漏掉重要因素和所有危险因素要在设计时考虑全面，避免漏掉重要因素和加入无关因素。加入无关因素。生存时间的定义要明确。起始事件和终点事件要有明确规生存时间的定义要明确。起始事件和终点事件要有明确规定定,时间尽可能精确测量。时间尽可能精确测量。研究的协变量在研究对象中的分布要适中，否则会给参数研究的协变量在研究对象中的分布要适中，否则会给参数的估计带来困难。的估计带来困难。在设计时要注意影响时间

12、的效应因素。在设计时要注意影响时间的效应因素。如研究吸烟对肺癌患者生存率的影响，若本来吸烟的人因患肺癌而戒烟，则一般模型不易分析吸烟的作用，需考虑吸烟量随时间变化的趋势。如果研究的因素随时间而发生变化，必要时可以采用伴时协变量的Cox模型进行分析。41ppt课件 注意事项(2)模型拟合时应注意的问题模型拟合时应注意的问题 多重共线性问题：医学研究中的许多变量间并不是独立的，多重共线性问题：医学研究中的许多变量间并不是独立的，但通常不会影响分析的结果，如果变量间存在高度的相关，但通常不会影响分析的结果，如果变量间存在高度的相关，则会影响则会影响Cox模型的参数估计，此时可采用主成分分析法或模型的

13、参数估计，此时可采用主成分分析法或R型聚类分析法消除多元共线性的影响。型聚类分析法消除多元共线性的影响。应注意应注意Cox模型要求病人的风险函数与基础风险函数呈比例，模型要求病人的风险函数与基础风险函数呈比例，如果这一假定不成立，则不能用如果这一假定不成立，则不能用Cox模型进行分析模型进行分析。42ppt课件 注意事项(3)模型应用时应注意的问题模型应用时应注意的问题 结果的解释应结合专业知识。结果的解释应结合专业知识。Cox模型与其他回归分模型与其他回归分析一样，当进入模型中的因素有统计学意义时，该因析一样，当进入模型中的因素有统计学意义时，该因素与生存时间不一定有因果关系，其中有一部分因

14、素素与生存时间不一定有因果关系，其中有一部分因素与生存时间的关系为伴随关系。与生存时间的关系为伴随关系。Cox回归的生存率一般不宜用于不同资料之间的比较，回归的生存率一般不宜用于不同资料之间的比较，因为基准危险率函数只在同一份资料内保持相同，不因为基准危险率函数只在同一份资料内保持相同，不同资料的基准危险率往往不同；而且在多因素分析的同资料的基准危险率往往不同；而且在多因素分析的情况下，协变量组合也很难一致，可比性难以保证。情况下，协变量组合也很难一致，可比性难以保证。43ppt课件3.Cox模型的局限性模型的局限性 lCox模型估计参数时，首先要假定偏似然函数具有最模型估计参数时，首先要假定

15、偏似然函数具有最大似然的性质，这个问题在理论上尚不完善。大似然的性质，这个问题在理论上尚不完善。lCox模型对异常值较为敏感，所以在进行模型配合时模型对异常值较为敏感，所以在进行模型配合时要注意要注意Cox模型拟合优度的检验。模型拟合优度的检验。lCox模型估计参数时，不是利用精确的生存时间，而模型估计参数时，不是利用精确的生存时间，而是利用生存时间的顺序统计量，这损失了一定的样本是利用生存时间的顺序统计量，这损失了一定的样本信息。当引进的协变量随时间的变化剧烈时，偏似然信息。当引进的协变量随时间的变化剧烈时，偏似然函数损失的信息也增多。函数损失的信息也增多。44ppt课件五、扩展的五、扩展的Cox回归模型回归模型1.伴时协变量Cox回归模型2.分层Cox回归模型3.分组数据的Cox回归模型4.多状态Cox回归模型45ppt课件46ppt课件