14.医学统计学生存分析.pptxVIP

生存分析(2)由于生存曲线只是对样本生存过程的统计描述样本生存率或生存曲线不同也可能是抽样误差所致，在生存分析中一个重要的问题是进一步对总体的生存曲线进行假设检验log-rank检验属于非参数方法该法并不指定生存时间服从某种特定的分布所比较的是整个生存时间的分布，而不是仅仅比较某个特定时间点的生存率例23-3 就例23-1数据，比较甲乙两疗法组脑瘤患者的生存率有无差别？两种治疗方式的脑瘤患者生存曲线相同两种治疗方式的脑瘤患者生存曲线不同两组在不同时点的期初观察例数列于第3、第7列，其合计列于第11列不同时点两组的死亡人数，分列于第4、第8列，其合计 列于第12列注意到不同时点期初观察例数等于上一时点期初观察例数减去上一时点的死亡数与删失数之后按下式计算各组期望死亡人数并分列于表中第5、第9列由于假设成立时，两组的生存分布相同，故可把两组的数据合并计算合并的死亡概率，以此计算相应的期望死亡人数，故将两组的完全生存时间混合排序列在表23-5第2列相同生存时间只列1次k=1,2 表示在每个时点，当两组的死亡率相等且均为该时点的总死亡数除以该时点的总观察数时，按照各组期初例数计算的期望死亡人数分别将两组各时点期望死亡人数相加列在第5、第9列合计处，而两组实际总死亡数为第4、第8列合计如果两组各时点生存率都相等，那么两组总的期望死亡数和总的实际死亡数相差不大注意事项（1）以上介绍的是log-rank检验的近似法，计算简便，但其结果较精确法（一般统计软件中输出精确法计算结果）保守。log-rank检验精确法 统计量计算公式为：wi为权重，对log-rank检验，wi=1，即该检验给任意时间点处两组间死亡的差别相同的权重。当比较的两总体生存曲线呈比例时，检验效能最大；wi=ni 则对应Gehan检验（1965）或Wilcoxon检验，该检验给两组间死亡的早期差别更大的权重（2）对于大样本频数表形式的生存曲线比较，基本方法与上述相同。另外，该法很容易推广到多个组的比较，在此不赘述。需要强调的是，生存曲线的比较也和前面均数、率的比较一样，要求组间具有可比性，最好是按照比较因素进行随机化分配之后再比较。如果是未经随机化分配的观察对比资料，要考虑是否有混杂因素干扰，若存在混杂因素，可进行分层分析或采用多因素分析方法（见本章第四节）。（3）当假设检验发现组间生存曲线有差别时，可通过中位生存期、相对危险度RR（relative risk）等指标评价其差别。相对危险度是两个对比组相对死亡比的比值，而相对死亡比是实际死亡数与期望死亡数之比，于是，第i组相对于第j组的相对危险度为：（4）log-rank检验用于整条生存曲线的比较，若比较两组某时间点处的生存率，则按下式计算：如比较多个时间点处生存率，检验水准可取Bonferroni校正，以保证总的I型错误概率不超过 趋势检验多组生存率比较时，若分组变量是等级变量，如肿瘤分期为Ⅰ期、Ⅱ期、Ⅲ期，或连续变量等级化分组，如年龄（岁）＜30、30～、40～、≥50，在log-rank检验组间生存率差别有统计学意义后，还可作趋势检验（trend test），分析危险率是否有随分组等级变化而变化的趋势。即是否有肿瘤分期越高，预后越差，或年龄越大（或越小），预后越差的情况。（1）按某种因素影响大小将病人分组 可采用临床上该因素的自然分组，如疾病的分期等，组数一般取奇数，如3组或5组。（2）计算每组的实际死亡数 与期望死亡数 。（3）进行趋势检验。 Cox比例风险回归模型log-rank检验属于生存时间比较的单因素分析方法，应该注意的是生存分析中的单因素比较在实验设计方面和均数、率的比较一样，要求对比组之间在非处理因素方面具有可比性。一般而言，经过随机化分配处理的实验设计数据在处理组之间可比性较好在实际工作中多见的观察对比资料其可比性通常不能满足或者研究者关心的影响生存时间的因素不止一个，此时应采用适当的多因素分析方法由于生存分析问题中反应变量比较特殊，是事件结局以及出现这一结局所经历的时间，普通的线性回归和logistic回归通常并不适用。如果仅考虑生存时间作为反应变量进行线性回归，由于生存时间通常并不是正态分布，不满足线性回归的模型要求；仅考虑某一时点事件结局作为反应变量进行logistic回归，生存时间长短的信息又未能充分利用；生存时间资料中还有删失数据的问题，上述两种模型都不能够利用这种不完全数据提供的信息模型结构式中x的 表示研究者认为可能影响生存的诸因素，也称协变量（covariates），这些变量在随访期间的取值不随时间变化而变化，例如根据研究目的可以是随访对象的年龄、性别、接受的不同治疗方式等。 t表示生存时间，好h(t,x) 称为具有协变量x 的个体在t 时刻的风险函数（hazard function），