医学统计学考题按章节第6题05分随访资料的生存分析

1、五、其它30分（35道题目,每题610分）随访资料的生存分析：【06真题】九、某医生从2002年1月1日起对某医院收治的 6名急性心肌梗塞病人进行跟踪观察，2002年3月25日结束观察，共12周。记录的资料如下：（5分）-1始规豪时眉时:!l!坝義忖冋12002 1 1元于0桎22002 1.1 2002.3 J 532002.1.232002.2.1吒干车届42002.1 2920023 J 5仍右活5200X124&2002 2.142002J.101、上述资料随访时间单位以（日）、（月）、（年）哪个较合适？为什么?2、判断上述随访时间哪些属截尾值？写出观察对象编号。【05真题、04真题、

2、03真题】四、16例某癌症病人在不同时期经随机化分配到A、B两治疗组，并继续进行随访至1974年5月31日结束。资料如下表：（8分）16例某种癌症病人随访资料病人号治疗组分组日期终止日期是否该病死亡截尾值1A68. 05. 1268. 05. 30Y2B70. 10. 1871. 04. 16Y3B69. 02. 1270. 11. 06Y4A72. 01. 3074. 05. 31仍存活5A73. 11. 1174. 01. 02Y6B68. 03. 1273. 03. 30车祸死亡7A69. 01. 0669. 01. 04Y8A69. 02. 0870. 02. 08迁出9B71 . 0

3、5. 0271. 11 . 13Y10B68. 03. 0868. 05. 23Y11B73. 12. 1274. 02. 20Y12A74. 05. 0174. 05. 09Y13B72. 07. 0272. 07. 15Y14B68. 12. 1874. 04. 31失访15A69. 01. 0174. 05. 31仍存活16B73. 09. 0273. 09. 20Y1. 上述资料随访时间单位以（日）、（月）、（年）哪个较合适？为什么？2. 判断上述随访时间哪些属截尾值，写出观察对象编号。3. 要比较A、B疗法对该种癌症病人的疗效，宜选用何种统计检验方法？4. A B治疗组随访资料生存时

4、间的特征量（代表值）一般用何指标表示？【答案】jszb0、本资料中，第7号观察对象数据， 终止日期 竟然早于 分组日期，是典型的错误数据，应该排除。1、本资料并未按时间分组，实际上是A、B两个治疗组的未分组资料。一般情况下 较细的时间单位准确性较 高，当随访时间可以作较细的量化时，则应考虑用较细的时间单位。 但研究目的不同，时间单位不同，使用恰当的时间单位。本资料的目的是 比较A、B疗法对该种癌症病人的疗效，癌症病人的生存时间测度单位 如果以（日）太小,因此，本资料随访时间单位以月较合适？2、产生截尾数据的原因：包括 中途失访、研究结束时仍然存活、死于与研究疾病无关的原因。 因此，编号4、6、

5、8、14、15观察对象属截尾值。3、 两样本生存曲线比较， 宜选用 对数秩检验 (log-rank test) 。表示。4、A、B治疗组随访资料生存时间的特征量(代表值)一般用 中位生存时间 median survival time【知识点解析】 Jszb应变量 有a. 生存分析是将事件的结果和出现这一结果所经历的时间，结合起来分析的一种统计分析方法，所以它的 两个，即 生存时间和结局 。b. 随访资料做生存分析的条件为 有一定的死亡数 。一、随访资料的生存分析( Survival analysis)(一)定义是将事件的 结局 和出现这一结局所经历的时间 结合起来分析的统计分析方法。临床科研中

6、，评价某种药物的疗 效或者某种疾病的愈后，不仅可以对某各时点治愈率(或生存率)进行评价，还可以考察生存时间。利用生存分析的原理和方法，把“出生” 、“死”的涵义稍加变通，在医学研究中可以有广泛的应用。 【意义】1、分析放置的心脏起搏器有效持续时间。2、研究抛弃性隐形眼镜使用寿命。3、比较某因素暴露组与非暴露组人群某种疾病发生或死亡的威胁性。 一般而言，能够明确划分起止点的现象均可以用生存分析的原理和方法来进行研究。(二)特点 能同时分析有结局的完全数据和没有结局的截尾数据，充分利用信息。完全数据( complete data )：是确切的生存时间。若以死亡为结局事件，则完全数据是指从研究起点至

7、死亡 所经历的时间。截尾数据( censored data )：由于某种原因，在研究结束时未能观察到病人的明确结局，所以不知道该病人 的确切生存时间。是从起点至研究结束所经历的时间。产生截尾数据的原因：1. 中途失访2. 在研究结束时仍然存活3. 死于与研究疾病无关的原因 生存分析常有两个效应变量：结局，生存时间1) “死亡”：即处理失败，终止随访时间为“死亡”时间。2) 中途失访：如失去联系、中途退出等，终止随访时间为最后一次访问时间为准。3) 死于其它与研究疾病无关的原因，终止随访时间为死亡时间。4) 随访截止 : 随访研究结束时观察对象仍存活，终止随访时间为研究结束时间。【基本术语】1、

8、“死亡”事件或称失败事件 (failure event)2、截尾值 (censored value)3、生存时间 (survival time)4、生存率 (survival rate)1、“死亡”事件或称失败事件 (failure event)反映处理因素失败或失效的特征。 如乳腺癌病人手术后复发、 肾移植病人肾功能衰竭、 白血病患者化疗后的复发等。2、截尾值 (censored value) 有的观察对象终止随访不是由于失败事件发生，而是由于中途失访、死于其它原因、随访截止。由于不知道这 些观察对象发生失败事件的时间，他们的资料不能提供完全的信息，这些对象的观察值称为截尾值，常用符号“+”

9、表示。3、生存时间 (survival time) 即随访观察持续的实足时间，按天、周、月、年等时间单位记录，常用符号 t 表示。 一般情况下 较细的时间单位准确性较 高，当随访时间可以作较细的量化时，则应考虑用较细的时间单位。但对 非住院病人常难以知道准确的死亡或失访时间。生存分析中最基本的问题就是计算生存时间，它有三个要素： 时间的起点、终点、时间的测度单位(时、日、月、年)3、生存时间 (survival time)如某病人 1990 年 2 月 1 日进入随访， 1992 年 4 月间发生失败事件，他的生存时间为 t=26 月；又如某白血病患者化疗 3月后失去联系，他的随访结果为一截尾

10、值，生存时间记为t=3月;t=2+ 年。安放心脏起搏器患者术后2年因意外事故死亡，他的随访结果也为一截尾值，生存时间记为 4、生存率(survival rate)指某个观察对象活过 t时刻的概率，常用 P(Xt)表示。如P(X10)表示某对象活过10天(或10月、10年)的概 率。根据不同随访资料的失败事件，生存率可以是缓解率、有效率等。1. 死亡概率：q【对于未分组资料】某时刻的死亡概率=某时刻的死亡人数/该时刻以前的期初观察人数未分组资料生存率的估计主要是用乘积极限法(product-limit estimates);又称Kaplan-Meier法，主要适用于观察例数不多时，是一种非参数方

11、法。【对于 分组资料】某年死亡概率=某年内死亡数/某年年初观察人数若有截尾，则分母用校正的年初观察人数(年初观察人数-1/2截尾人数)2. 生存概率：p = 1 q3. 生存率：S(ti)S(ti)=ti时刻仍存活的例数/观察总例数如果含有截尾数据，分母需要校正。生存率为累积生存概率，是各时段上的生存概率乘积。S(ti) = p1*p2* . *pi4. 生存率的标准误5ES(r;) = 5(OxE-V (-心5. 生存曲线以时间t为横轴，生存率 S(ti)为纵轴，表达了生存率随时间的变化趋势。常为阶梯形。(四)生存分析研究的内容1. 描述生存过程。估计生存率及平均生存时间，绘制生存曲线。 未

12、分组资料，采用Kaplan Meier法(简称K M法)0分组资料，采用寿命表法(Life table )。2. 比较组间生存曲线是否有差别。对于未分组资料，采用对数秩检验(log rank test )。3. 影响生存时间的因素分析。COX莫型。【总体生存率可信区间的估计】 用正态近似原理估计某时点总体生存率的可信区间。nP0 U Snpo如本例3个月时生存率的95%可信区间为：下限:0.2724 1.96 0.04570.1828上限：0.2724 1.96 0.04570.3620(一) K-M法描述生存过程：CasePiocessing SiiininaiyTotal NN of Ev

13、ents715 iCensoredPercent | | 28.6% |Censored:删失，截尾事件发生截尾截尾比例 的例数 例数Survival TableTimeStatusCumulative Proportion Surviving at the TimeNof Cumulative EventsNof Remaining CasesEstmateStd Error185.0001.857.132162125.0001253125 0001.571.187344125.0000335183.0001.199426258 0001.190.168517433.000050各时刻的生存

14、率生存率的标准误Cumulative proportion surviving at the time:累积生存比例，简称“生存率 JMeans and Medians foi SurviViil TimeMedianEstimateStd. Error95% Confidence IntervalEstimateStd. Error95% Confidence IntervalLcwer BoundUpper BoundLower B oundUpper Bound214.33348.429119413309.254I 183.00060.69964.030301.970a. Estimat

15、ion is limited tt the largest survival time if it is censored.生存时间的均数.标准误及95%生存时间的中位数.标准误及95% 置信区间置信区间Petcenlles25.0%60.0%75.0%EstimateStd. ErrorEstimateStd ErrorEstimateStd Error 1258.00065.975183.00060.699125 00026.166 I生存时间的P25及生存时间的P75及标准误标准误Survi irul r urictiDn1 QD 0生存+6 d0 oIAlAJnIcnUJlloppui

16、rvalFunction+ Cerso? largast survive I lime if ilis censored.OvierdI Comparisons由图可见，随时间的延拴，两组生存率逐渐下降， 第一组（A药）生存率低于第二组（B药），但无统 计学基异（可能和样本量太小有关）.【五、两样本生存率曲线的比较】1、两时点生存率的比较对于大样本资料，根据正态近似原理，可用u检验。 该法是正态近似法，当样本的生存率很小或很大，如小于0.05，或大于0.99，或各组死亡人数、生存人数少于5时，该方法误差较大； 由于该法只检验某时点两样本生存率的差别有无统计学意义，而不能对整个生存期作出评价。2

17、、 两样本生存曲线比较对数秩检验（log-ra nk test）其基本思想是如果无效假设成立，即两总体生存曲线无差别时，则两种处理在各个时期的实际死 亡数与理论死亡数不会相差太大，否则应认为无效假设不可能成立，两条生存率曲线差异有统计学意 义。Log-rank检验：时序检验，属于非参数方法，比较整个生存时间的分布 大样本卡方检验2条或多条生存曲线比较对数秩检验（log-rank test） 公式：2 (A T)2x T 组数1当有T v 5时，用校正公式。H0:单放疗组与放化疗组的生存率曲线分布相同 H:单放疗组与放化疗组的生存率曲线分布不同=0.052计算检验统计量：X = 4.65， 得P

18、=0.03按 =0.05水准，拒绝H0，接受H，故认为单放疗组与放化疗组的生存率曲线分布不同。【生存率比较注意事项】当比较不同治疗方案的治疗效果时，仅比较某个时点的生存率，有时可能出现不正确的结论。事实上， 每一种防治措施的效果都是一条 时间效应曲线，所以评价该措施的效果应对这条时间效应曲线进行 全 面评价，而不是仅对某个时点的生存率进行评价。1. 生存分析：生存分析(survival analysis )是将事件的结果和出现这一结果所经历的时间，结合起来分析的一种统计分析方法，它不仅可以从事件结局的好坏，如疾病的治愈(成功)和死亡(失败)，而且可以从事件的持续时间，如某病经治疗后存活的时间长

19、短进行分析比较，因而能够更全面、更精确地反映该治疗的效果。2. 生存时间：生存时间(survival time)指观察到的存活时间。3完全数据：完全数据(complete data )指从起点至死亡所经历的时间，即死者的存活时间。4截尾数据：由于失访、改变防治方案、研究时间结束时事件尚未发生等情况，使得部分病人不能随访到底，称之为截尾。从起点至 截尾所经历的时间，称为截尾数据(censored data )。5. 死亡率：某年内死亡例数与年中观察例数之比称为死亡率( mortality rate )。6. 死亡概率：死亡概率(mortality probability )是指某年内死亡例数与年

20、初观察例数之比，若年内有截尾，分母用校正人口数。7. 生存率：生存率(survival rate )指病人经历tk个单位时间后仍存活的概率，即tk时刻仍存活的例数与观察总例数之比。8. 生存概率：生存概率(survival probability)表示在某单位时段开始时存活的个体到该时段结束时仍存活的机会大小，它是某年 活满一年人数与年初观察例数之比，若年内有截尾，分母用校正人口数。中位生存时间 median survival time :半数生存时间/平均生存时间，恰好由50%个体存活s(t)=0.5 的时间生存时间一般为非正态，故用中位生存期。如果生存期最长的死亡对象的生存率高于50%则无

21、法估计中位生存期。（三）寿命表法描述生存过程:lit色 Nible9inlivdl 51hf1 TinsNurnber EnleririQ inirn 创Number Wthdravmg dunngi iririNJ nub erEjipcsed to RiskriunrhErOT Temmal tvyntsProportisnTEnri ridingProportitnSunMliijCunulaW PccortiDni StinirlnQ ai End tl irileivaiTd Eiroror CuinulatN Pre ip 匚 rticn Sur/rving a!End Jf im&rvdii0255202562.