资源描述
肝手术病人生存时间的预测模型
第三队
白 璇 罗基莉 张晓辉
(女) (女) (男)
肝手术病人生存时间的预测模型
摘 要 本文针对肝手术病人四项指标与生存时间的关系数据进行了分析计算. 运用多项式拟合的方法,建立多元线性回归模型,并经过多次改进,最终得出病人生存时间的预测模型.
首先,利用题中所给数据,分别做出凝血值,预后指数,酵素化验值,肝功化验值与病人生存时间的散点图,由图像拟合得出此四项指标与生存时间的线性回归模型. 在此基础上,考虑交叉项的影响,建立多元非线性回归模型一,利用MATLAB计算出模型一的自变量系数,置信区间和检验统计量,依此判断模型中各参量对生存时间的影响程度,得出交叉项对生存时间影响不太显著,没有必要引入模型.
其次,不考虑交叉项的影响做出以生存时间为因变量,四项指标为自变量的多元线性回归模型二,利用与模型一同样的求解方法得出肝功化验值对生存时间的影响不显著,故将其在模型中剔除,建立改进后的模型三. 利用同样的方法计算出结果并剔除异常数据后,最终得出预测模型四.
再次,利用表数据拟合模型四得到预测生存时间,与真实生存时间进行对比分析,从而验证了模型四的合理性和精确度. 得出凝血值,预后指数,酵素化验值对生存时间的灵敏度依次增大的结论.
最后,本文对模型进行了分析评价, 然后将该模型推广到其它预测问题中.
关键词 线性回归 检验统计量 预测模型
一.问题重述
为了对做过某一类型肝手术病人的生存时间做预测,某医院外科随机的选取了位需要做此手术的病人为研究对象.对每位病人手术前考察了下列四个指标:凝血值;预后指数;酵素化验值;肝功化验值.根据题目中表 给出的数据,建立合适的数学模型,以预测病人手术后的生存时间.并解答以下问题:
1.考虑自变量交互项对模型结果的影响,并判断有无将其引入模型的必要.
2.在此分析的基础上建立模型.模型建立后,需说明模型的合理性.
3.使用残差图分析的方法,分析剔除肝功能化验值对模型结果的影响,以改进模型.
4. 评估关于误差正态性假定的合理性.
5.引入其他有用的数据进行分析.
根据题目中表所给的数据,首先拟合问题1.2中所选择的模型,得出相关数据,并与原来所得数据进行比较,以评价所选择模型的稳定性和预测能力.
最后,综合表和表中的数据,对所选择模型的精细和确认性分析,对所得结果进行评价.
二. 问题分析
为较为准确地预测做过某一类肝手术病人的生存时间,需运用曲线拟合的方法,建立以病人手术前的四项指标(凝血值,预后指数,酵素化验值,肝功化验值)为回归变量,手术后生存时间对数变换值为因变量的回归模型.
针对问题1.1,要确定最终的回归模型中是否引入自变量交叉项,需考虑引入自变量交叉项对因变量的影响程度,如果影响不显著,则不需引入;反之,则需要引入.
针对问题1.2,建立合理模型不仅需考虑交叉项对模型结果的影响,还需考虑某项指标对模型结果的影响,是否引入该项亦须考虑其影响的显著程度. 验证模型合理性只需将自变量值代入预测模型当中,比较预测值与实际值间的差值,即可验证.
针对问题1.3,根据问题1.2的求解过程,可画出剔除肝功化验值前后模型的残差图. 对比两个残差图,若变化明显,则不可剔除肝功化验值;反之,则可剔除.
针对问题1.4,可对问题1.2中的残差值画出残差图和残差正态分布图,观察图像即可评估误差正态性假定的合理性.
针对问题1.5,肝手术后生存时间还与其他指标相关,搜索相关指标的数据,并将其作为新的自变量进行曲线拟合,对模型进行修正.
针对问题2.1,首先运用曲线拟合的方法对问题一中所选模型进行拟合,并与模型一中各自变量的回归系数以及检验统计量进行比较,以此判断所选模型的合理性.
针对问题2.2,要评价模型的预报能力,须将表2中的自变量代入问题一中的预测模型中,计算出预测值与真实之间的误差,以此为依据评价所选模型的准确性.
针对问题2.3,要对模型做出精细性分析,可以考虑自变量的灵敏度. 分别剔除一个自变量后,计算出模型的检验统计值,与模型三作比较,以此为参考得出各自变量的灵敏度.
三. 基本假设
1. 仅考虑题中所给指标对手术后生存时间的影响;
2. 患者无其他病变和意外情况发生;
3. 残差大致符合均值为零的正态分布.
四. 符号表示
病人编号
第位病人的凝血值
第位病人的预后指数
第位病人的酵素化验值
第位病人的肝功化验值
病人生存时间对数变换值
病人生存时间对数变换预测值
模型拟合产生的随机误差
普通残差. 即实际观测值与回归值的偏差
模型中自变量的系数
回归模型的决定系数
模型的统计量值
统计量值对应的概率值
模型的置信水平
随机误差的方差
方差的无偏估计值
五. 模型建立与求解
根据题目所给表格,确定四项指标与病人生存时间对数值的关系,然后分析模型的合理性并给出预测.
首先考察各自变量对因变量的影响,为了大致地分析与,,,的关系,利用表1的数据分别作出对,,,的散点图(见图1).
0
5
10
15
1.5
2
2.5
3
0
50
100
1.5
2
2.5
3
0
50
100
150
1.5
2
2.5
3
0
2
4
6
8
1.5
2
2.5
3
图1 对,,,的散点图
从图可以发现,随着的增加,值有较明显的线性增长趋势,图中的直线是用线性模型拟合得到,综合四个图可建立线性回归模型如下:
式右端的称为回归变量,是给定凝血值,预后指数,酵素化验值,肝功化验值时,生存时间对数的平均值,其中的参数为回归系数,由表的数据估计,影响的其它因素作用都包含在随机误差中,如果模型选择的合适,应大致服从均值为零的正态分布.
5.1 问题1.1的求解:
为考察自变量交互作用对病人生存时间的影响,不妨建立包括各交互项的回归模型,通过对交互项回归系数及置信区间的观察,判断有无必要将这些交叉项引入到模型中.
模型一建立:
从图可以发现,随着,,,的增加,的值有比较明显的线性增长趋势,故建立如下的回归模型:
在这个模型中,的均值与及它们的交叉项有关,由系数确定.
模型一求解:
直接利用MATLAB统计工具箱中的命令regress求解,得到模型一的回归系数及其置信区间(置信水平)、检验统计量,,的结果见表:
表 模型一的计算结果
参数
参数估计值
参数置信区间
表显示,指因变量(生存时间的对数变换值)的可由模型确定,值远远超过检验的临界值(),值远小于,因而模型一从整体来看是可用的.
表的回归系数给出了模型一中的到的估计值.检查它们的置信区间发现,的置信区间均包含零点(且区间两端端点距零点很近),表明回归变量的各交叉项对因变量的影响不显著,故各交叉项可不引入到模型中.
5.2 问题1.2的求解:
由问题1.1中的分析可得,在建立回归模型时,可不引入各交叉项.即只需建立与回归变量线性相关的回归模型.
模型二建立:
回归变量,,,对因变量的影响是相互独立的,即病人生存时间的对数平均值与凝血值,预后指数,酵素化验值,肝功化验值的关系分别由回归系数,,,确定,得到:
在这个模型中,的均值与给定凝血值,预后指数,酵素化验值,肝功化验值有关.
模型二求解:
直接利用MATLAB统计工具箱中的命令regress求解,得到模型二的回归系数及其置信区间(置信水平)、检验统计量,,的结果见表:
表 模型二的计算结果
参数
参数估计值
参数置信区间
表显示,指因变量(生存时间的对数变换值)的可由模型确定,值远远超过检验的临界值(),值远小于,因而模型二从整体来看是可用的.
表的回归系数给出了模型一中的,,,,的估计值,即,,,,.检查它们的置信区间发现,只有的置信区间包含零点(但区间左端点距零点很近),表明回归变量(对因变量的影响)不太显著,模型需进行改进.
对模型二各回归变量的回归系数以及置信区间进行考察可以看出,(肝功化验值)对影响不太显著. 故在改进模型时,可将剔除.
模型三建立:
如果剔除肝功能化验值,模型变为:
在这个模型中,的均值仅与给定凝血值,预后指数,酵素化验值有关,由系数确定.
模型三求解:
利用MATLAB统计工具箱中的命令regress求解,得到模型三的回归系数及其置信区间(置信水平)、检验统计量,,的结果见表:
表 模型三的计算结果
参数
参数估计值
参数置信区间
对比表和表可以看出,剔除后的值变化微小,而且的值增大,所有参数的置信区间不包含零点,由此判断,剔除肝功能化验值对模型的预测影响不大,但相比模型二更合理.
为了对残差进行分析,图给出模型三的残差图.
图 模型三的残差图
从图可发现两个异常点,为了使个别的数据不影响整个模型,应该将这个异常数据去掉,对模型三重新估计回归系数,得到的结果如表,残差图见图,可以看出,去掉异常点数据后结果有所改善.
图 模型三去掉异常数据后的残差图
表 模型三去掉异常数据后的计算结果
参数
参数估计值
参数置信区间
对比表和表 可以看出,剔除两个异常数据后,得到的和均有所增大,所有参数的置信区间不包含零点,由此可以得出该问题的最优预测模型,称为模型四:
将组数据代入预测模型得到的预测值与真实值进行比较(比较结果见附录表),得出平均偏差为,所以,模型三去除两个异常点后建立的预测模型四相对符合实际,更合理.
5.3 问题1.3的求解:
模型四中已经剔除掉肝功化验值,为了对残差进行分析,下面给出了模型二和模型三的残差图:
图 模型二的残差图 图 模型三的残差图
对比模型二和模型四的残差图可看出,模型二和模型四的残差图差别很小,而且图中的残差分布相对集中,说明肝功化验值对预测模型的影响很小,可以不考虑.
5.4 问题1.4的求解:
对模型三的残差从小到大排序后,利用MATLAB统计工具箱中提供的图检验其是否符合正态分布,模型四的残差图如下图:
图 模型三残差图检验
由图中可以看到,模型三的残差点近似分布在一条直线上,则可认为它符合正态分布,则误差正态性假定是合理的,下面给出模型三的残差正态分布图:
图 模型三残差正态分布图
5.5 问题1.5的求解:
做过肝手术病人的生存时间除了和凝血值,预后指数,酵素化验值,肝功化验值有关外,还与抗凝指标值,纤溶指标值等有关系.
5.6 问题2.1求解:
问题一中所选模型为
根据表中的数据,运用MATLAB统计工具箱中的命令regress重新求解,得到模型三的回归系数及其置信区间(置信水平),检验统计量,,的结果见表:
表 用表中的数据拟合模型三的计算结果
参数
参数估计值
参数置信区间
得到预测模型如下:
利用式对表中所给的指标进行预测,所得预测值与真实值进行比较(见附录表),平均偏差为,表明所选模型比较稳定,未出现较大程度的偏差. 再将表与表中的各项数据进行比较,发现各项数据的变化较小,说明模型三的精确度高,稳定性好.
对表数据的残差图(图)进行分析可见数据分布集中,这也说明模型三精确度高,能够较准确的预测病人生存时间.
图 用表数据拟合模型三所得残差图
5.7 问题2.2的求解:
对一个实际问题建立起多元线性回归方程后,一个重要应用就是利用方程去预测. 除了要知道预测值外,还希望知道预测精度,即更希望能给出的一个预测范围(预测区间),使这个区间包含的概率达到指定的.
随机误差的方差的无偏估计值为:
则的置信水平达到和的近似预测区间分别为:
和
由上式可知,近似预测范围较大,故模型精确度高,即模型准确性高,预报能力强.
5.8 问题2.3的求解:
为了对模型做出更精细的分析,可分别考虑模型三中凝血值,预后指数,酵素化验值的灵敏度.
在模型三中剔除凝血值后,模型变为:
运用MATLAB统计工具箱中的命令regress求解,得到该模型的回归系数及其置信区间(置信水平),检验统计量,,的结果见表:
表 模型三中剔除的计算结果
参数
参数估计值
参数置信区间
对比表和表可以发现,去除后和都有明显的下降,参数置信区间也增大,从而可以得出对模型的灵敏度较大,即凝血值对病人的生存时间影响较大.
在模型三中剔除预后指数后,模型变为:
运用MATLAB统计工具箱中的命令regress求解,得到该模型的回归系数及其置信区间(置信水平),检验统计量,,的结果见表:
表 模型三中剔除的计算结果
参数
参数估计值
参数置信区间
对比表和表可以发现,去除后和都有明显的下降,参数置信区间也增大,从而可以得出对模型的灵敏度较大,即预后指数对病人的生存时间影响较大.
在模型三中剔除酵素化验值后,模型变为:
运用MATLAB统计工具箱中的命令regress求解,得到该模型的回归系数及其置信区间(置信水平)、检验统计量,,的结果见表:
表 模型三中剔除的计算结果
参数
参数估计值
参数置信区间
对比表和表可以发现,去除后和都有很大程度的下降,参数置信区间也明显增大,从而可以得出对模型的灵敏度很大,即酵素化验值对病人的生存时间影响很大.
综上分析可以得出,凝血值,预后指数,酵素化验值三项指标对病人生存时间的影响都很显著,且灵敏度依次增大,即酵素化验值的灵敏度最高,对病人生存时间的影响最大.
六.结果分析
对肝病患者进行手术前,需先对其相关身体指标进行化验检查,以确定患者是否适合进行肝脏手术以及有无进行手术的必要.本题中给出病人的四项指标数据和病人生存时间,通过模型的建立和分析,得出肝功化验值对生存时间的影响很小,最终得出预测模型. 该预测模型精确度较高,可以通过其它三项指标数值计算出病人的生存时间,给医学上的预测带来很大的方便.
七.模型改进
本题所选模型分析过程相对简单,在消除交叉项影响的分析过程中,仅考虑了全部引入的情况,没有进行逐步引进交叉项时回归系数变化情况的分析,现使用MATLAB工具箱中的stepwise工具,通过对各自变量以及交叉项引入前后回归系数和偏差变化的分析,最后确定是否引入这些项.
使用stepwise工具所得结果如图和图,由它们的对比可以发现,仅考虑,,时,模型结果与考虑全部自变量及交叉项时的模型结果相近.
其中,图中Stepwise Table窗口中的数据为:
,,;
图中Stepwise Table窗口中的数据为:
,,.
检验统计量的比较也显示只考虑,,时,所得模型合理.
图 仅考虑,,时模型的stepwise图
图 考虑,,,及交叉项时模型的stepwise图
八.模型评价与推广
本文建立的模型和算法具有一般性,适用于预测肝病人手术后的生存时间,模型运用多项式拟合方法,经过逐步分析和不断改进,最终建立多元线性回归模型,并计算出自变量的系数,从而得到预测模型. 经检验,该模型的精确度较高,可以较为准确地预测肝病人手术后的生存时间.
题目中只给出四项指标对生存时间的影响,但在实际生活中,仍有其它因素影响病人的生存时间. 本文建立的模型可以推广到其它多元非线性回归问题上,并做出相关预测,同时也可在增加其它数据和条件的基础上广泛应用于预测问题中.
参考文献
[1]刘振航,数学建模[M],北京:中国人民大学出版社,2004
[2]刘新平 魏启恩,概率论与数理统计[M],西安:西安出版社,2002
[3]李宏艳 王雅芝,数学实验[M],北京:清华大学出版社,2007
[4]任玉杰,数值分析及其MATLAB实现[M],北京:高等教育出版社,2008
[5]肖树铁,数学实验[M],北京:高等教育出版社,1999
附录
表1 54位肝手术病人的观测数据
病人编号
凝血值
预后指数
酵素化验值
肝功化验值
生存时间
对数变换值
1
6.7
62
81
2.59
200
2.3010
2
5.1
59
66
1.70
101
2.0043
3
7.4
57
83
2.16
204
2.3096
4
6.5
73
41
2.01
101
2.0043
5
7.8
65
115
4.30
590
2.7067
6
5.8
38
72
1.42
80
1.9031
7
5.7
46
63
1.91
80
1.9031
8
3.7
68
81
2.57
127
2.1038
9
6.0
67
93
2.50
202
2.3054
10
3.7
76
94
2.40
203
2.3075
11
6.3
84
83
4.13
329
2.5172
12
6.7
51
43
1.86
65
1.8129
13
5.8
96
114
3.95
830
2.9191
14
5.8
83
88
3.95
330
2.5185
15
7.7
62
67
3.40
168
2.2253
16
7.4
74
68
2.40
217
2.3365
17
6.0
85
28
2.98
87
1.9395
18
3.7
51
41
1.55
34
1.5315
19
7.3
68
74
3.56
215
2.3324
20
5.6
57
87
3.02
172
2.2355
21
5.2
52
76
2.85
109
2.0374
22
3.4
83
53
1.12
136
2.1335
23
6.7
26
68
2.10
70
1.8451
24
5.8
67
86
3.40
220
2.3424
25
6.3
59
100
2.95
276
2.4409
26
5.8
61
73
3.50
144
2.1584
27
5.2
52
86
2.45
181
2.2577
28
11.2
76
90
5.59
574
2.7589
29
5.2
54
56
2.71
72
1.8573
30
5.8
76
59
2.58
178
2.2504
31
3.2
64
65
0.74
71
1.8513
32
8.7
45
23
2.52
58
1.7634
33
5.0
59
73
3.50
116
2.0645
34
5.8
72
93
3.30
295
2.4698
35
5.4
58
70
2.64
115
2.0607
36
5.3
51
99
2.60
184
2.2648
37
2.6
74
86
2.05
118
2.0719
38
4.3
8
119
2.85
120
2.0792
39
4.8
61
76
2.45
151
2.1790
40
5.4
52
88
1.81
148
2.1703
41
5.2
49
72
1.84
95
1.9777
续表
42
3.6
28
99
1.30
75
1.8751
43
8.8
86
88
6.40
483
2.6840
44
6.5
56
77
2.85
153
2.1847
45
3.4
77
93
1.48
191
2.2810
46
6.5
40
84
3.00
123
2.0899
47
4.5
73
106
3.05
311
2.4928
48
4.8
86
101
4.10
398
2.5999
49
5.1
67
77
2.86
158
2.1987
50
3.9
82
103
4.55
310
2.4914
51
6.6
77
46
1.95
124
2.0934
52
6.4
85
40
1.21
125
2.0969
53
6.4
59
85
2.33
198
2.2967
54
8.8
78
72
3.20
313
2.4955
表2 肝手术病人生存时间的模型确认数据
病人编号
凝血值
预后指数
霉素化验值
肝功化验值
对数变换值
1
7.1
23
78
1.93
2.0326
2
4.9
66
91
3.05
2.4086
3
6.4
90
35
1.06
2.2177
4
5.7
35
70
2.13
1.9078
5
6.1
42
69
2.25
2.0035
6
8
27
83
2.03
2.0945
7
6.8
34
51
1.27
1.7652
8
4.7
63
36
1.71
1.7925
9
7
47
67
1.6
2.1295
10
6.7
69
65
2.91
2.2295
11
6.7
46
78
3.26
2.1524
12
5.8
60
86
3.11
2.3188
13
6.7
56
32
1.53
1.9039
14
6.8
51
58
2.18
2.0508
15
7.2
95
82
4.68
2.6525
16
7.4
52
67
3.28
2.2053
17
5.3
53
62
2.42
1.9246
18
3.5
58
84
1.74
2.1541
19
6.8
74
79
2.25
2.497
20
4.4
47
49
2.42
1.7237
21
7
66
118
4.69
2.8339
22
6.7
61
57
3.87
2.1282
23
5.6
75
103
3.11
2.6884
24
6.9
58
88
3.46
2.4284
25
6.2
62
57
1.25
2.0261
26
4.7
97
27
1.77
2.0843
27
6.8
69
60
2.9
2.2826
28
6
73
58
1.22
2.2073
续表
29
5.9
50
62
3.19
2.0443
30
5.5
88
74
3.21
2.4863
31
3.8
55
52
1.41
1.9037
32
4.3
99
83
3.93
2.6647
33
6.6
48
54
2.94
1.9071
34
6.2
42
63
1.85
1.9093
35
5
60
105
3.17
2.4389
36
5.8
62
82
3.18
2.3343
37
4.7
42
10
0.28
1.3379
38
5.7
70
59
2.28
2.1996
39
4.7
64
48
1.3
1.8795
40
7.8
74
40
2.58
2.1504
41
2.9
43
32
0.94
1.433
42
4.9
72
90
3.51
2.4381
43
4.6
73
57
2.82
2.1075
44
5.9
78
70
4.28
2.2843
45
4.6
69
70
3.17
2.1615
46
6.1
53
52
1.84
2.0558
47
5.9
88
98
3.33
2.7249
48
4.7
66
68
1.8
2.052
49
10.4
62
85
4.65
2.681
50
5.8
70
64
2.52
2.2604
51
5.4
64
81
1.36
2.2553
52
6.9
90
33
2.78
2.1745
53
7.9
45
55
2.46
2.0224
54
4.5
68
60
2.07
2.1413
表11 模型四预测值与真实值的比较
1
2.301
2.3002
0.0008
29
1.8573
1.8731
0.0158
2
2.0043
2.0093
0.005
30
2.2504
2.1464
0.104
3
2.3096
2.3254
0.0158
31
1.8513
1.9065
0.0552
4
2.0043
1.9938
0.0105
32
1.7634
1.7229
0.0405
5
2.7067
2.7409
0.0342
33
2.0645
2.0706
0.0061
6
1.9031
1.928
0.0249
34
2.4698
2.4432
0.0266
7
1.9031
1.9053
0.0022
35
2.0607
2.0613
0.0006
8
2.1038
2.1361
0.0323
36
2.2648
2.2745
0.0097
10
2.3075
2.3363
0.0288
37
2.0719
2.1595
0.0876
11
2.5172
2.4909
0.0263
38
2.0792
2.0063
0.0729
12
1.8129
1.8277
0.0148
39
2.179
2.1036
0.0754
13
2.9191
2.8674
0.0517
40
2.1703
2.1831
0.0128
14
2.5185
2.4943
0.0242
41
1.9777
1.9844
0.0067
15
2.2253
2.2359
0.0106
42
1.8751
1.9412
0.0661
16
2.3365
2.3331
0.0034
43
2.684
2.7403
0.0563
续表
17
1.9395
1.9391
0.0004
44
2.1847
2.1919
0.0072
18
1.5315
1.5894
0.0579
45
2.281
2.3138
0.0328
19
2.3324
2.33
0.0024
46
2.0899
2.1149
0.025
20
2.2355
2.2333
0.0022
47
2.4928
2.4849
0.0079
21
2.1335
2.0509
0.0826
48
2.5999
2.5761
0.0238
23
1.8451
1.8452
0.0001
49
2.1987
2.1899
0.0088
24
2.3424
2.3291
0.0133
50
2.4914
2.4937
0.0023
25
2.4409
2.43
0.0109
51
2.0934
2.0864
0.007
26
2.1584
2.1471
0.0113
52
2.0969
2.0859
0.011
27
2.2577
2.1489
0.1088
53
2.2967
2.2903
0.0064
28
2.7589
2.8439
0.085
54
2.4955
2.5107
0.0152
表12 用表2中数据拟合模型四预测值与真实值比较
编号
真实值
预测值
偏差
编号
真实值
预测值
偏差
1
2.0326
1.9422
0.0904
29
2.0443
1.9455
0.0988
2
2.4086
2.2998
0.1088
30
2.4863
2.3787
0.1076
3
2.2177
2.0839
0.1338
31
1.9037
1.7416
0.1621
4
1.9078
1.8721
0.0357
32
2.6647
2.479
0.1857
5
2.0035
1.9552
0.0483
33
1.9071
1.9006
0.0065
6
2.0945
2.0924
0.0021
34
1.9093
1.9042
0.0051
7
1.7652
1.7586
0.0066
35
2.4389
2.3883
0.0506
8
1.7925
1.7245
0.068
36
2.3343
2.2415
0.0928
9
2.1295
2.0466
0.0829
37
1.3379
1.2812
0.0567
10
2.2295
2.2056
0.0239
38
2.1996
2.0839
0.1157
11
2.1524
2.1224
0.03
39
1.8795
1.85
0.0295
12
2.3188
2.2621
0.0567
40
2.1504
2.0885
0.0619
13
1.9039
1.7672
0.1367
41
1.433
1.373
0.06
14
2.0508
1.9812
0.0696
42
2.4381
2.3447
0.0934
15
2.6525
2.6434
0.0091
43
2.1075
2.012
0.0955
16
2.2053
2.1211
0.0842
44
2.2843
2.2779
0.0064
17
1.9246
1.9293
-0.0047
45
2.1615
2.1017
0.0598
18
2.1541
2.0575
0.0966
46
2.0558
1.8904
0.1654
19
2.497
2.3942
0.1028
47
2.7249
2.6405
0.0844
20
1.7237
1.6832
0.0405
48
2.052
2.0622
-0.0102
21
2.8339
2.7142
0.1197
49
2.681
2.6045
0.0765
22
2.1282
2.0552
0.073
50
2.2604
2.1397
0.1207
23
2.6884
2.549
0.1394
51
2.2553
2.2209
0.0344
24
2.4284
2.3431
0.0853
52
2.1745
2.1008
0.0737
25
2.0261
2.028
-0.0019
53
2.0224
1.9773
0.0451
26
2.0843
1.9466
0.1377
54
2.1413
1.9883
0.153
27
2.2826
2.1644
0.1182
平均值
0.074854
28
2.2073
2.1233
0.084
内容总结
22
展开阅读全文