医学统计学考题第6题

1、五、其它30分（35道题目，每题610分）随访资料的生存分析：【06真题】九、某医生从2002年1月1日起对某医院收治的6名急性心肌梗塞病人进行跟踪观察，2002年3月25日结束 观察，共12周。记录的资料如下：（5分）观察对象开始观察时间坤束观察时间12002.1.1M0E.2.4死于心F22002.1.12002.3.25仍存活.=2OT2.2.16死于车病42002.1.232002.3.25仍存活20022.132002.2.24失防6200?2.142002.3.10死于心槌1、上述资料随访时间单位以（日）、（月）、（年）哪个较合适？为什么？2、判断上述随访时间哪些属截尾值？写出观察

2、对象编号。【05真题、04真题、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.0

3、8迁出9B71. 05. 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. 20Y上述资料随访时间单位以（日）、（月）、（年）哪个较合适？为什么？判断上述随访时间哪些属截尾值，写出观察对象编号。要比较A、B疗法对该种癌症病人的疗效，宜选用何种统计检验方法？A、B治疗组随访资料生存时间的特征量（代表值）一般用

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

5、存曲线比较，宜选用对数秩检验(log-rank test)4、A、B治疗组随访资料生存时间的特征量(代表值)一般用中位生存时间median survival time表示。【知识点解析】Jszb生存分析是将事件的结果和出现这一结果所经历的时间，结合起来分析的一种统计分析方法，所以它的应变量有 两个，即生存时间和结局。随访资料做生存分析的条件为有一定的死亡数一、随访资料的生存分析(Survival analysis)定义是将事件的结局和出现这一结局所经历的时间结合起来分析的统计分析方法。临床科研中，评价某种药物的疗 效或者某种疾病的愈后，不仅可以对某各时点治愈率(或生存率)进行评价，还可以考察生

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

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

8、ate)1、“死亡”事件或称失败事件(failure event)反映处理因素失败或失效的特征。如乳腺癌病人手术后复发、肾移植病人肾功能衰竭、白血病患者化疗后的复发等。2、截尾值(censored value)有的观察对象终止随访不是由于失败事件发生，而是由于中途失访、死于其它原因、随访截止。由于不知道这 些观察对象发生失败事件的时间，他们的资料不能提供完全的信息，这些对象的观察值称为截尾值，常用符号+” 表示。3、生存时间(survival time)即随访观察持续的实足时间，按天、周、月、年等时间单位记录，常用符号t表示。一般情况下较细的时间单位准确性较高，当随访时间可以作较细的量化时，则

9、应考虑用较细的时间单位。但对 非住院病人常难以知道准确的死亡或失访时间。生存分析中最基本的问题就是计算生存时间，它有三个要素：时间的起点、终点、时间的测度单位(时、日、月、年)3、生存时间(survival time)如某病人1990年2月1日进入随访，1992年4月间发生失败事件，他的生存时间为t=26月；又如某白血病患者化疗3月后失去联系，他的随访结果为一截尾值，生存时间记为t=3月；安放心脏起搏器患者术后2年因意外事故死亡，他的随访结果也为一截尾值，生存时间记为t=2+年。4、生存率(survival rate)指某个观察对象活过t时刻的概率，常用P(Xt)表示。如P(X10)表示某对象

10、活过10天(或10月、10年)的概 率。根据不同随访资料的失败事件，生存率可以是缓解率、有效率等。死亡概率：q【对于未分组资料】某时刻的死亡概率=某时刻的死亡人数/该时刻以前的期初观察人数未分组资料生存率的估计主要是用乘积极限法(product-limit estimates)；又称Kaplan-Meier法，主要适用于观 察例数不多时，是一种非参数方法。【对于 分组资料】某年死亡概率=某年内死亡数/某年年初观察人数若有截尾，则分母用校正的年初观察人数(年初观察人数一1/2截尾人数)生存概率：p=1 q生存率：S(ti)S(ti)=ti时刻仍存活的例数/观察总例数如果含有截尾数据，分母需要校正

11、。生存率为累积生存概率，是各时段上的生存概率乘积。S(ti) =p1*p2* . *pi生存率的标准误生存曲线以时间t为横轴，生存率S(ti)为纵轴，表达了生存率随时间的变化趋势。常为阶梯形。(四)生存分析研究的内容描述生存过程。估计生存率及平均生存时间，绘制生存曲线。未分组资料，采用KaplanMeier法(简称KM法)。分组资料，采用寿命表法(Life table )o比较组间生存曲线是否有差别。对于未分组资料，采用对数秩检验(logrank test)。影响生存时间的因素分析。COX模型。【总体生存率可信区间的估计】用正态近似原理估计某时点总体生存率的可信区间。np 土 u s如本例3个

12、月时生存率的95%可信区间为：0.2724 -1.96 x 0.0457 = 0.18280.2724 +1.96 x 0.0457 = 0.3620下限： 上限：K-M法描述生存过程:Case Pi ocessiiHi SuimniyI1- Censored:删失，截尾Total N QfEM8nt NT T s Hl 28E5 I事件发生截尾截尾比例 的例歉 例数Survival TableTimeStatusCUEul白tl白 Propcrtion Surviving artth? Timen or Cumulative EventsN of Remaining CasesEstimat

13、eSid. Error105.0001.857132162125.0001253125.0001.671187344125.0000335183.0001.36119942658.0001.190.108517433.000050各时刻的生存率生存率的标准误Cumulative proportion surviving at the time:累积生存比例，简称“生存率：aiidModnn&rof SiifwalTimeMean*MedianEstimateStd. Error95% Confidence IntervalE 乱 rnateStd. Error35% Confidence In

14、tervalLower BoundUpperBoundUwer BoundUpper Bound21 4 33348.420119.413309.254163.00060.69S64.030301.970a. Estimation is limited ttllie lareet survive I time ifitiB censored.生存时间的均数、标准误及95%生存时间的中位数、标准误及95% 置信区间置信区间Pefcentlles25.0%75.0%EstimateStd ErrarEstimate8tdl. ErrorEstimateStd. Error259D0CG5.9751

15、33.000so.sag125.0C02B.1S611生存时间的P&及生存时间的P75及标准误标准误生存率Survival FjncIbiDnpjSurvivalI Function-j- Censored1D0 2OJ 3Q0400| hniE I生存时间由图可见，随时间的延长，生存率逐渐下降。图中“ + ”表示载尾敷据.（二）比较组间生存曲线是否有差别：分别给出两组在各时刻的生存率估计值及标准误。Case Processing SuiniihiiygroupTotal NN ofEntsCensor&dTTPercent175128.6%776111.3%OverallU11321/%Su

16、ivvtil Tableg rou pTimeStaluscumulate rreport on Surviving at the TimeMof Cumulative EventsFJ OfRemaining CasesEstimateEtd. Error1185.0001.顷/32162125.0JD12521 25.0001.571.1373441 25.000a33b1 83.0001.381.1994I6268.0001.19D.168617433.000a60217.0301.857.1321621 38.00017U.171253I 82.000i.571.1ST344236.0

17、001.423.187435235.030a42fion mn1.211.178517477.0001.ODD.00060两组生存时间的均数，中位数及相应的标准误、95%置信区间；hi-SquaredfSig.Leg Rank (Mantel-Cox)1一6晶ITest of equelitf of suivival cistribufions forthe differentle/eI group.对数秩检验11 D-D.B-DED.4-D-2-必值SuFYival Functionsgroup_n=+ 1-ceiwred+ 2-censcredD1lM2D030D400500由图可见随时间

18、的延长，两组生存率遂渐下降，第一组（A药）生存率低于第二组（B药），但无统计学差异（可能和样本量大小有关）【五、两样本生存率曲线的比较】1、两时点生存率的比较对于大样本资料，根据正态近似原理，可用u检验。该法是正态近似法，当样本的生存率很小或很大，如小于0.05，或大于0.99，或各组死亡人数、 生存人数少于5时，该方法误差较大；由于该法只检验某时点两样本生存率的差别有无统计学意义，而不能对整个生存期作出评价。2、两样本生存曲线比较对数秩检验（log-rank test）其基本思想是如果无效假设成立，即两总体生存曲线无差别时，则两种处理在各个时期的实际死 亡数与理论死亡数不会相差太大，否则应认

19、为无效假设不可能成立，两条生存率曲线差异有统计学意 义。 Log-rank检验：时序检验，属于非参数方法，比较整个生存时间的分布大样本卡方检验 2条或多条生存曲线比较对数秩检验（l og-rank test） 公式：Means and Medians for Survival Timegroup . a MeanMedianEstimateStc. Error95% ConfidencE IntervalEstimateSid. Error&5% Confidence InternalLower BoundUrper BoundLort，er BoundJiper Bound1214.3334

20、S.429H 19.41 3309.2541S3.00060.69964.030301.9702284.00063.955158.844409.156235.00069.393gs.geg371.011Overall250.58343.268165.779335.3881E3.00046.11792.610273.390a. Estimation is limited to the largest sun/iy/al time if it is censored.Overall Compel isonsV 二组数-1当有TV5时，用校正公式。H：单放疗组与放化疗组的生存率曲线分布相同H0：单放

21、疗组与放化疗组的生存率曲线分布不同a = 0.05计算检验统计量：X2 = 4.65,得P=0.03按a = 0.05水准，拒绝H。，接受H1，故认为单放疗组与放化疗组的生存率曲线分布不同。【生存率比较注意事项】当比较不同治疗方案的治疗效果时，仅比较某个时点的生存率，有时可能出现不正确的结论。事实上， 每一种防治措施的效果都是一条时间效应曲线，所以评价该措施的效果应对这条时间效应曲线进行全 面评价，而不是仅对某个时点的生存率进行评价。生存分析：生存分析(survival analysis )是将事件的结果和出现这一结果所经历的时间，结合起来分析的一种统计分析方法， 它不仅可以从事件结局的好坏，

22、如疾病的治愈(成功)和死亡(失败)，而且可以从事件的持续时间，如某病经治疗后存活的时间长短 进行分析比较，因而能够更全面、更精确地反映该治疗的效果。生存时间：生存时间(survival time)指观察到的存活时间。完全数据：完全数据(complete data)指从起点至死亡所经历的时间，即死者的存活时间。截尾数据：由于失访、改变防治方案、研究时间结束时事件尚未发生等情况，使得部分病人不能随访到底，称之为截尾。从起点至 截尾所经历的时间，称为截尾数据(censored data)。死亡率：某年内死亡例数与年中观察例数之比称为死亡率(mortality rate)。死亡概率：死亡概率(mortality probability )是指某年内死亡例数与年初观察例数之比，若年内有截尾，分母用校正人口数。生存率：生存率(survival rate)指病人经历tk个单位时间