学习资料生存数据分析.pptVIP

寿命表法估计生存率步骤如下: 1.将观察例数按时间段(年)0-,1-,2-,划分,分别计数期初例数,死亡,失访, 截尾例数列入表22.2的1-5列。事实上,从第二个时间段开始,期初人数ni 系由下式算得: ni=ni-1-di-ui-wi 例如第二行,即时间段1-,有 n2=126-47-4-15=60 2.计算各时间段期初实际观察例数,(亦称有效例数)ni ni=ni-ui/2-wi/2 上式表明该时间段期初例数中的失访,及截尾例数只计其半时,即得有效例数。 如第一行,n1=126-4/2-15/2=116.5 3.分别用(22.5)(22.6)(22.7)式计算条件生存率^S(ti/ti-1),累积生存率s(ti)及其标准误。 计算结果已列于表22.2中,第7,8,11列,表中9,10二列系用于第11列的计算。 例如时间段0--中 ^S(ti/ti-1)=(116.5-47)/116.5=0.5966 ^S(ti)=1×0.5966=0.5966 SE(S(ti))=0.5966×√5.805×10-3=0.0455 故一年生存率的估计为0.5966±0.0455 同样二年生存率的估计为0.5386±0.0479 由于寿命表法与积限法的累积生存率及其标准误的计算公式完全相同,所以,当分组资料中每一个分组区间中最多只有1个观察值时,寿命表法就是积限法。 课件 第四节 生存率的比较 当有两个或两个以上的生存分布时,我们常需比较它们是否来自同一生存分布,此时的假设检验为: H0:样本所来自的总体生存分布相同。 H1:样本所来自的总体生存分布不相同。 可选用的检验方法有:Logrank法, 广义Wilcoxon法,和Cox-Mantel法等。当拒绝H0时,认为几个生存分布不相同。 课件 当不需要整体比较，而只要比较个别时间点上几组生存率时可用下面方法： （1）两个生存率比较 生存率S1和S2，其方差为V1和V2 用卡方检验： ?2 =（S1-S2）2 / （V1+V2） df=1 课件 （2）两个以上两个生存率比较 生存率S1，S2和S3，方差为V1，V2和V3 用卡方检验： 权重W1=1/V1, W2=1/V2，W3=1/V3 加权平均生存率： S=（W1*S1+W2*S2+W3*S3） /（W1+W2+W3） ?2 =W1*(S1-S)2 ++W2*(S2-S)2 +W3*(S3-S)2 df=3-1 课件 课件 §3.1 Logrank检验(Log Rank Test) 当比较的几个样本生存分布,全部为完全数据时,本检验又称为Savage检验。 Logrank检验的计算步骤如下: 1.将两样本的生存数据混合,由小到大排列,并给以秩次i1, 当截尾数据与完全数据数值相同时,截尾数据排列在后。并设两样本含量分别为m1,m2,总例数n=m1+m2。 例22.1中药组与对照组生存数据排列结果见表22.3中第1,2列。 2.列出所比较的两组中任一个组的序号i2(本处选用中药组),记入表22.3中第3列。 3.列出死亡例的序号i3(见表22.3中第4列)。 4.计算非截尾数据(完全数据)各时间点处于危险状态的例数r,它表示该时刻时还剩下多少例数。r系由与i3相应的i1值计算而得。 r=n-i1+1 (22.10) 例如与生存期7(月)相应的r值系由 r=26-9+1=18 算得,见表中第5列。 课件 5.对秩次i1作logrank变换,即计算logrank变换值W,其算法为 秩次为i1 的序号为i3非截尾数据的W值为 i3 W=∑ 1/rj-1 (22.11) j=1 秩次为i1的截尾数据,首先判断它在那二个非截尾数据之间,如果它在序号i3与i3+1之间则W为 i3 W=∑ 1/rj (22.12) j=1 特别地,当截尾数据在第一个非截尾数据之前时,取W=0,几个截尾数据落在同样序号的非截尾数据之间时,它们具有相同的W值。 例如表22.3中第1个数据为非截尾(已死亡)则由(22.11)式得 W=1/26-1=-0.96 第2,3数据都是截尾数据(存活),