预防医学/小样本病例随访资料统计分析
随访病例较少时,可按下法求不同时期的生存率(或缓解率)及其统计学意义分析。
一、资料统计方法和曲线描绘分析
例23.3某单位用甲、乙两法治疗何杰金病。甲法治疗15例中已复发9例;乙法治疗14例,有4例复发。两组随访情况如表23-3。
先以甲疗法为例说明不同随访时期的缓解率及其标准误。演算结果如表23-4。
表23-4 甲、乙两法治疗何杰金病随访天数
甲疗法 | 乙疗法 | ||
已复发者 | 尚未复发者 | 已复发者 | 尚未复发者 |
141 | 1446+ | 505 | 615+ |
364 | 836+ | 296 | 570+ |
950 | 498+ | 1375 | 1205+ |
570 | 173+ | 688 | 1726+ |
312 | 1540+ | 1190+ | |
570 | 836+ | 822+ | |
173 | 1408+ | ||
401 | 1493+ | ||
86 | 1645+ | ||
1570+ |
尚未复发者随访天数后加“+”号,表明缓解天数至少多于随访天数
表23-4 甲疗法治疗何杰金病不同时期缓解率计算
病序(1) | 随访天数n(2) | 复发例数r(3) | 期初病例数R(4) | 复发概率qx(5) | 缓解概率px(6) | 累计缓解概率np0(7) | 标准误snp0(8) |
1 | 86 | 1 | 15 | 0.0667 | 0.9333 | 0.933 | 0.064 |
2 | 141 | 1 | 14 | 0.0714 | 0.9286 | 0.867 | 0.088 |
3 | 173 | 1 | 13 | 0.0769 | 0.9231 | 0.800 | 0.103 |
4 | 173 | … | 12 | 0.0000 | 1.0000 | 0.800 | - |
5 | 312 | 1 | 11 | 0.0909 | 0.9091 | 0.727 | 0.117 |
6 | 364 | 1 | 10 | 0.1000 | 0.9000 | 0.654 | 0.126 |
7 | 401 | 1 | 9 | 0.1111 | 0.8889 | 0.581 | 0.131 |
8 | 498+ | … | 8 | 0.0000 | 1.0000 | 0.581 | — |
9 | 570 570 |
2 | 7 | 0.2857 | 0.7143 | 0.415 | 0.136 |
10 | |||||||
11 | 836 836 |
… | 5 | 0.0000 | 1.0000 | 0.415 | — |
12 | |||||||
13 | 950 | 1 | 3 | 0.3333 | 0.6667 | 0.277 | 0.145 |
14 | 1446+ | … | 2 | 0.0000 | 1.0000 | 0.277 | — |
15 | 1540+ | … | 1 | 0.0000 | 1.0000 | 0.277 | - |
1.按随访天数从小到大依次排列,如遇复发者天数和未复发者随访天数相同时,以复发者排在前面。
2.填写不同随访天数的复发例数及期初病例数如表23-4的(3)、(4)栏。
3.求出不同随访天数的复发概率qx(复发例数÷期安病例数)和缓解概率px(1-qx)如(5)、(6)栏。
4.根据公式(23.6)求出累计缓解概率np0如(7)栏。
5.按下式求不同时点累计缓解率的标准误。
公式(23.8)
本例173天时点累计缓解率的标准误:
同法可以求得乙疗法的累计缓解率及其标准误,学者试自演算求解。
6.缓解率曲线描绘以横轴为随访天数(n),纵轴为累计缓解率(np0),将两疗法的演算结果各点的坐标准确标出,然后将各点向右连成与横轴平行的阶梯形,得出两组缓解曲线如图23-1。可以看出乙疗法累计缓解率水平始终在甲法之上。
图23-1 甲、乙疗法累计缓解率的比较
二、两疗法差异的统计学意义分析
如果要分析两疗法差异有无统计学意义,可用时序检验法(log rank test)。假定两组疗法效果相同,求各时点预期复发数,再进一步作x2检验。演算如表23-5。
表23-5按检验假设算得甲、乙两组的预期复发数(即理论值)和实际数,分别为:
A甲=9,T甲=5.138;A乙=4,T乙=7.817
代入x2检验公式
查x2值表,x20.05(1)=3.84,今x2>4.675,P<0.05,表明两法累计缓解率曲线的差别有统计学意义。
表23-5 甲、乙两疗法预期复发数计算表
疗法分组(1) | 观察天数(2) | 复发例数 | 期初病例数 | 预期复发数 | |||||
甲组(3) | 乙组(4) | 合计(5)=(3)+(4) | 甲组(6) | 乙组(7) | 合计(8)=(6)+(7) | 甲组(9)=(5)(6)/(8) | 乙组(10)=(5)(7)/(8) | ||
甲 | 86 | 1 | 1 | 15 | 14 | 29 | 0.517 | 0.483 | |
甲 | 141 | 1 | 1 | 14 | 14 | 28 | 0.500 | 0.500 | |
甲 | 173 | 1 | 1 | 13 | 14 | 27 | 0.481 | 0.519 | |
甲 | 173+ | … | 12 | 14 | 26 | … | … | ||
乙 | 296 | 1 | 1 | 11 | 14 | 25 | 0.440 | 0.560 | |
甲 | 812 | 1 | 1 | 11 | 13 | 24 | 0.458 | 0.542 | |
甲 | 364 | 1 | 1 | 10 | 13 | 23 | 0.435 | 0.565 | |
甲 | 401 | 1 | 1 | 9 | 13 | 22 | 0.409 | 0.591 | |
甲 | 498+ | … | 8 | 13 | 21 | … | … | ||
乙 | 505 | 1 | 1 | 7 | 13 | 20 | 0.350 | 0.650 | |
甲 甲 |
570 >570 |
1 >2 1 |
1 > 1 |
7 | 12 | 19 | 0.737 | 1.263 | |
乙 | 570+ | … | 5 | 12 | 17 | … | |||
乙 | 615+ | … | 5 | 11 | 16 | … | |||
乙 | 688 | 1 | 1 | 5 | 10 | 15 | 0.333 | 0.667 | |
乙 | 822+ | … | 5 | 9 | 14 | … | … | ||
甲 | 836+ > 836+ |
… … > … |
5 | 8 | 13 | … | … | ||
甲 | |||||||||
甲 | 950 | 1 | 1 | 3 | 8 | 11 | 0.273 | 0.727 | |
乙 | 1190+ | … | 2 | 8 | 10 | … | … | ||
乙 | 1205+ | … | 2 | 7 | 9 | … | … | ||
乙 | 1375 | 1 | 1 | 2 | 6 | 8 | 0.250 | 0.750 | |
乙 | 1408+ | … | 2 | 5 | 7 | … | … | ||
甲 | 1446+ | … | 2 | 4 | 6 | … | … | ||
乙 | 1493+ | … | 1 | 4 | 5 | … | … | ||
甲 | 1540+ | … | 1 | 3 | 4 | … | … | ||
乙 | 1570+ | … | 0 | 3 | 3 | … | … | ||
乙 | 1645+ | … | 0 | 2 | 2 | … | … | ||
乙 | 1726+ | … | 0 | 1 | 1 | … | … | ||
总和 | (A)9 | (A)4 | 13 | 15 | 14 | 29 | (T)5.183 | (T)7.817 |