第三节 小样本病例随访资料统计分析
随访病例较少时,可按下法求不同时期的生存率(或缓解率)及其统计学意义分析。
一、资料统计方法和曲线描绘分析
例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
一、资料统计方法和曲线描绘分析
例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