基于多元回归下的白血病不同疗法白细胞生存数目对比

长春工业大学实习报告题目基于多元回归下的白血病不同疗法白细胞生存数目对比学院基础科学学院班级120802班姓名王春光实习时间2015年11301225实习地点林园图书馆7楼机房指导教师李纯静2015年12月25日实习题目2015年11月30日2014年12月25日班级120802学号20112986姓名王春光成绩_指导教师签字_【实习目的】我们通过系统的学习，在最后的时光里利用这次宝贵的实习机会，通过对数据的处理，老师的讲授学习处理实际问题。在论文中我们主要采用口服抗生素和输入红细胞，采用LOGRANK检验，根据检验的结果可以知道有显著意义的因素。用非参数的COX比例风险回归模型对肺结核治愈率进行单因素和多因素分析,当数据个体较少时，为充分利用每个数据的信息，必须采用更为精确的估计方法，所以我们采用KAPLANMEIER（也叫乘积限估计）方法来估计生存函数，然后画出生存函数图。收集到符合入组条件的30名急性白血病患者，用其生存资料来分析这两种治疗方案的生存情况。生存分析方法适用于急性白血病疗效指标的动态描述及其影响因素的分析。【数据来源】来自吉大一院的肿瘤科【使用软件】SAS【实习内容】1；首先是数据的采集处理2；利用统计知识建立模型3；进行模型拟合4；检验拟合结果5。导出模型陈述结果【研究过程】1查找数据我们将30名急性白血病患者分别以口服抗生素和输入红细胞两种方法进行治疗，治疗后继续对这些人进行随访，至两年半后结束，数据录入见附录所示，试进行COX模型分析。所采用的数据是右删失数据，根据所做的调查研究中，要考察和比较不同的治疗方法，不同病理类型，病人的某些特征对疾病预后的影响，可运用非参数方法进行组间比较，这类方法虽然使用简单，但在多因素共存条件下，单一因素的比较会受到其他因素干扰和混杂，组间难以达到均衡，而且不能分析和考察因素间的关系（如交互作用）和进行定量评价。因此，我们采用COX回归模型。2模型建立生存分析还可以估计生存率函数STEKTM、危险率函数KTMKKTM1,拟合生存率曲线。常用的方法是采用非参数的KAPLANMEIER估计生存率。生存分析回归模型由基本部分H0T和修改部分XX为预后因素X的函数两部分组成。这里协变量XX1,X2,XP可以是不同的治疗方法，不同病理类型或病人的性别、年龄等特征，也可以是复合变量（如交互作用）。当基本部分H0T与修改部分为乘积关系时，即危险函数由基本危险X率函数乘上一个协变量的常数因子组成。HT,XH0TX3程序输出结果分析在SAS程序中，我们把生存月数中的删失数据记为0，未删失数据记为1；采用口服抗生素的治疗方法记为A，采用输入红细胞的治疗方法为B；白细胞未损伤的状态记为1（表现为正号），白细胞已损伤的状态记为0（表现为负号）。程序调用PHREG过程分析，PHREG为PROPORTIONALHAZARDREGRESSION（比例风险回归）。PHREG和MODEL为COX模型分析的必须语句，MODEL语句左边为生存时间和截尾指示，右边为需分析的解释变量。【模型结果】生存分析还可以估计生存率函数STEKTM、危险率函数KTMKKTM1,拟合生存率曲线。常用的方法是采用非参数的KAPLANMEIER估计生存率。生存分析回归模型由基本部分H0T和修改部分XX为预后因素X的函数两部分组成。这里协变量XX1,X2,XP可以是不同的治疗方法，不同病理类型或病人的性别、年龄等特征，也可以是复合变量（如交互作用）。当基本部分H0T与修改部分为乘积关系时，即危险函数由基本危险X率函数乘上一个协变量的常数因子组成。HT,XH0TX称此模型为成比例危险模型（PROPORTIONALHAZARDSMODEL）。带协变量的修改部分最常见的形式为EXPX,则模型为XHT,XH0TEXPXH0TEXP（1X12X2PXP）其中（1，2，P）是P个未知待估计的回归系数，它是描述各个因素对生存期影响大小的参数向量；XX1,X2,XP是一个P维协变量向量，XP可以是影响预后的治疗方法、病人的某些特征，也可以是这些因素的交互作用项；H0T是未知的、任意分布的基准危险函数。在SAS程序中，我们把生存月数中的删失数据记为0，未删失数据记为1；采用口服抗生素的治疗方法记为A，采用输入红细胞的治疗方法为B；白细胞未损伤的状态记为1（表现为正号），白细胞已损伤的状态记为0（表现为负号）。程序调用PHREG过程分析，PHREG为PROPORTIONALHAZARDREGRESSION（比例风险回归）。PHREG和MODEL为COX模型分析的必须语句，MODEL语句左边为生存时间和截尾指示，右边为需分析的解释变量。【实证结论】COX模型程序结果输出如下PHREG过程模型信息数据集WORKCOX因变量T删失变量A删失值0结值处理BRESLOW读取的观测数使用的观测数3030事件和删失值个数汇总总计事件删失删失百分比302551667第1步，效应“B”已输入。模型包含以下效应收敛状态满足收敛准则GCONV1E8。模型拟合统计量准则没有协变量具有协变量2LOGL131800126436AIC131800128436SBC131800129654检验全局零假设BETA0检验卡方自由度PR卡方似然比73649100005评分76067100079WALD72641100018NOTENOADDITIONALEFFECTSMETTHE005LEVELFORENTRYINTOTHEMODEL最大似然估计值分析参数自由度参数估计值标准误差卡方PR卡方危险比B110151904929342415003942760C100407204934100068093421042程序第一步引入B变量，经过以上三种检验，发现都有差别，P599,有2205,统计学意义。两变量引进后不能剔除，建立COX回归模型。分别列出了变量名、自由度、参数估计值、标准误、WALD单变量值、P值和危险率比，两变量2WALD的与384比较，P皆小于005，所建立COX回归模型为20,EXP1590472HTBC【参考文献】1何晓群，刘文卿应用回归分析M北京中国人民大学出版社，20112高惠璇统计计算M北京大学出版社，20083汪海波，罗莉，汪海灵，徐俊杰。SAS统计分析与应用从入门到精通。北京人民邮电出版社，20064中国上海第二医科大生物统计教研室第十八章生存分析和COX回归2010附录数据编号生存月数疗法白细胞损伤12A224B36A411B54B65B717A820A98B104A111B1230A1316B1411B159B166B1715A1826A1913B203A2112B2214A2317B2418A256B2627A2712A2821A294B3013BSAS程序1、COX模型DATACOXINPUTTABCCARDS211124120611011121412051211701020111812041101120301111612011121912060201511126010131203110121201411117011181106120270111211121110412013120PROCPHREGMODELTA0BCRUN2、LOGRANK检验DATASAINPUTTABCIFA1THENCENSOR1ELSECENSOR0CARDS21112412061101112141205121170102011181204110112030111161