研究生医学统计学-生存分析 教程教案.ppt

19 生存分析; 生存分析的目的;;第一节　生存时间资料的特点 一、数据结构 二、统计描述指标 三、资料要求 ;一、数据结构;医学例子：起始事件 随访时间 终点事件;随访（follow-up）资料的记录;月 份（2000年） 图12-1 生存时间原始记录示意 (+截尾);生存时间资料的特点 ; 2个效应变量（1）生存时间(天数)，（2）结局(死亡与否、是否阳性等) 错误1：忽略生存时间，采用Logistic回归分析死亡率 错误2：忽略结局，采用t检验、线性回归分析生存时间;　　得不到确切的生存时间，但它们提供的生存时间长于观察期的时间，这种数据为不完全数据。或截尾数据、删失数据或终检值。 （如有确切的生存时间，则这种数据称为完全数据。） 两种错误的做法： 错误1：丢弃截尾数据，只考虑确切数据。（损失了信息） 错误2：将截尾数据当作确切数据处理。（低估了生存时间的平均水平）。 ;　　 两种错误的做法： 错误1：采用平均生存时间而不是采用中位生存时间来表示生存时间的平均水平。 错误2：采用常规t检验或方差分析进行组间比较。（应采用log-rank检验比较几组生存时间 ） ;二、生存分析的统计描述指标;[例1] 手术治疗60例肺癌病人，术后每年死亡10例，无删失。试求基本生存分析指标。N=60 ； ; 1 . 死亡率、死亡概率、生存概率 (1) 死亡率 (mortality rate，death rate) 表示某单位时间内的死亡强度。 年平均人口数=(年初人口数+年末人 口数)/2 ;(2) 死亡概率 ( mortality probability ) 是指单位时段开始时存活的个体到该时段结束时死亡的可能性，即死于某时段t～t+n内的概率。 ＝ d/n0 　　　　　　　　　　　 年内有删失，分母用校正人口数： 校正人口数= 年初人口数—删失例数 / 2 ＝n0－c/2 ;(3) 生存概率 ( survival probability ) 指某单位时段开始时存活的个体到该时段结束时仍存活的可能性的大小。 　　　　　　　　　　　　　　　　　　　　　　　　＝ (n0-d)/n0 分子为年末尚存人数， 若年内有删失，分母用校正人口数。;生存率 (survival rate) ( 累积生存概率 cumulative probability of survival ) 指病人经历t个单位时间后仍存活的概率。 生存概率指单个时段的概率， 生存率指从0～t多个时段的积累概率。 ;;如数据中无删失，生存率可用下式计算 　　　　　　　　　　　　　　　　　 如果有删失数据，则需分时段计算不同单位时间的生存概率Pi(i=1,2,…,t),然后利用概率乘法原理将Pi相乘得到t时刻生存率。 　S ( t ) = P (T ≥ tk ) = p1 p2 … pk ;例如　手术治疗50例肺癌病人，术后1，2，3年的 死亡数分别为10，10，10例，无截尾数据。 试求各年的 生存概率和3年生存率。 解： 各年生存概率 p1 = ( 50 – 10 ) / 50, p2 = ( 40 – 10 ) / 40, p3 = ( 30 – 10 ) / 30 3 年生存率 S(3) = P(T ≥3) ＝(n0-d)/N = ( 50 – 30 ) / 50 = 0.4000 　或 S(3) = p1 p2 p3 = 0.4000 ;生存率的标准误 生存率置信区间