艾滋病疗法的评价及疗效的预测模型
摘 要
本论文通过对影响艾滋病疗法的疗效的主要因素进行了综合评价,通过统计方法得到已知数据的特征变化趋势，利用灰色预测法对疗效做了分析预测,由此得出最佳治疗终止时间和最佳疗法．
针对问题一，对缺省数据，采用了插值和取截断值的方法进行恢复，按时间分段后统计得出病人体内CD4和HIV浓度的变化趋势（见文中图一、图二），并对这些数据进行分析和预测得出最佳治疗终止时间在25周左右．
针对问题二，将数据按整周点作为时间点划分，根据数据的特点将CD4浓度分段，统计出按时间和CD4分段后的人数百分比（文中表1至表4），分析数据对四种疗法进行综合评价，可得结论为：疗法四为最佳。
对此最优疗法进行预测，由于有很多不确定的因素，运用灰色预测法建立GM（1，1）
dx(1)?ax(1)??可得CD4浓度与时间的预测关系： 微分方程模型：dtx(k?1)?[x(0)(1)?^(1)?a]e?ak??a （k=0,1,2…n）由于残差检验观察相对误差序列有的
CD4
浓度与时间的预测关系：
误差很大，修正后
x(k?1)?[x(0)(1)?^(1)?a]e?ak??a??(k?1)(?ae)[e(0)(1)??eae]e?aek由此得到停止治疗时间为
第40周。
对于问题三，在考虑疗效的同时还应考虑费用问题，因此主要利用问题二得出的结论选择一个疗效比较好且费用较低的疗法，得出结论为：对于考虑经济效益的情况下最好选择疗法三。
关键词
灰色预测法、最小二乘法、关联度检验、后验差检验
1
问题重述
艾滋病是当前人类社会最严重的瘟疫之一，从1981年发现以来的20多年间，它已经吞噬了近3000万人的生命。
导致艾滋病的主要是一种叫做HIV的病毒破坏人的免疫系统，使人体丧失抵抗各种疾病的能力，从而严重危害人的生命。人类免疫系统的CD4细胞在抵御HIV的入侵中起着重要作用，当CD4被HIV感染而裂解时，其数量会急剧减少，HIV将迅速增加，导致AIDS发作。据统计，凡CD4+T淋巴细胞＜200/mm3或CD4+T淋巴细胞的百分比＜14%的HIV感染者可归入艾滋病。
艾滋病治疗的目的，是尽量减少人体内HIV的数量，同时产生更多的CD4，至少要有效地降低CD4减少的速度，以提高人体免疫能力。
据很多资料表示，成功的抗病毒治疗应该达到以下效果：血浆中HIV-RNA（HIV核酸法检验的病毒载量）的水平4周内应下降1个Log以上， 4～6个月内病毒降至检测不到的水平(HIV RNA<50拷贝/ml)，CD4+T淋巴细胞计数应逐渐上升。如果达不到相应指标，则应根据实验室检测结果调整治疗方案。
迄今为止还没有找到能根治AIDS的疗法，目前的一些AIDS疗法不仅对人体有副作用，而且成本也很高。本问题就是要确定一个较好的治疗方案。
现在得到了美国艾滋病医疗试验机构ACTG公布的两组数据。一组数据是服用了相同药物的300多名病人每隔几周测试的CD4和HIV的浓度（每毫升血液里的数量）。另一组数据是将1300多名病人随机地分为4组，每组按4种不同的疗法中的一种服药，大约每隔8周测试的CD4浓度（这组数据缺HIV浓度，它的测试成本很高）。根据提供数据完成以下问题：
（1）利用附件1的数据，预测继续治疗的效果，或者确定最佳治疗终止时间（继续治疗指在测试终止后继续服药，如果认为继续服药效果不好，则可选择提前终止治疗）。
（2）利用附件2的数据，评价4种疗法的优劣（仅以CD4为标准），并对较优的疗法预测继续治疗的效果，或者确定最佳治疗终止时间。
(3) 艾滋病药品的主要供给商对不发达国家提供的药品价格如下：600mg zidovudine 1.60美元，400mg didanosine 0.85美元，2.25 mg zalcitabine 1.85美元，400 mg nevirapine 1.20美元。如果病人需要考虑4种疗法的费用，对（2）中的评价和预测（或者提前终止）有什么改变。
模型建立与求解
问题(1):
将附录1中数据统计得出，可将时间分段，分为五个阶段：表一和表二中I表示0—4周，II表示4—8周，III表示8—25周，IV表示25—40周，V表示40—57周。 Ⅰ 1．69% Ⅱ 7．63% Ⅲ 13．10% Ⅳ 14．79% Ⅴ 11．54% ???,?50? 2
??50,0? 14．04% 42．70% 26．40% 10．67% 2．53% 1．97% 20．34% 50．56% 13．84% 5．08% 1．13% 1．41% 24．11% 34．23% 16．37% 8．04% 2．98% 1．19% 表一
23．47% 34．69% 18．88% 5．10% 2．04% 1．53% 19．23% 50．00% 19．23% 0．00% 0．00% 0．00% ?0,50? ?50,100? ?100,150? ?150,200? ?200,??? 表一说明：左列表示CD4的值的变化范围，为负表示CD4减少的值，为正表示CD4增加的值，数值表示所占的比例。
图一
图一说明：1—7分别表示表一中CD4值变化的七个范围，图一是随着服药时间的延续CD4在各范围内的增加或减少的人数变化趋势,CD4减少则看做是药品的副作用，经过统合评价得出结论在第三阶段和第四阶段总体上CD4在增加，因此最佳治疗终止时间在这两个阶段内。
Ⅰ 0．00% 0．57% Ⅱ 4．65% 7.56% Ⅲ 6．85% 16.20% Ⅳ 5．26% 12.91% ?1.5,??? ?0.5,1.5? 3
?0,0.5? 4.57% 12.57% 74.29% 8.00% 10.76% 63.66% 12.21% 1.16% 表二
8.41% 53.89% 13.18% 1.56% 20.57% 45.45% 14.53% 5．26% ??1,0? ??3,?1? ???,?3? 表二说明：左列表示HIV的值的变化范围，为负表示HIV值减少，为正表示HIV增加的值，数值表示所占的比例。
图二
图二说明：1—6分别表示表二中HIV值变化的6个范围，图二是随着服药时间的延续HIV在各范围内的增加或减少的人数变化趋势,HIV增加则看做药品的副作用，经过综合评价得出结论在第二阶段和第三阶段总体上HIV在减少，因此最佳治疗终止时间在这两个阶段内。
综合以上两方面的因素，确定最佳治疗终止时间为第三阶段，第三阶段为8—25周，所以我们认为最佳终止时间为开始服药后的第25周。
问题(2):
将采用不同治疗方法的病人分组进行分析，最后比较优劣。
根据检测间隔时间大约都是8周，我们将40周平均分为五个区间，将CD4浓度也等间隔划分为七个区间（具体划分见下表）。利用统计方法统计出处于各个时间段不同CD4浓度的人数。
由于这1300多个病人是随机分给各组的，直接比较没有可比性，故我们将统计的各个时间段的人数求和，用统计得到的各个人数值除以相应时间段的人数和，得到各个时间段CD4浓度不同的人数比值（表中数据）。通过比较四种方法的这些比值来确定疗效的优劣。
4
（比值用%表示） 表1 Log(CD4 count+1) [0,1) 检测时间（周） [1,2) [2,3) [3,4) [4,5) [5,6) [0,8] (8,16] (16,24] (24,32] (32,40] 表2 Log(CD4 count+1) 检测时间（周） 3.07 5.96 5.71 5.70 10.43 8.11 13.30 15.71 15.02 17.79 35.75 33.95 40.95 45.08 41.10 44.52 40.37 30.95 29.53 28.83 8.11 6.42 5.71 4.66 1.86 0.44 0.95 [0,1) [1,2) [2,3) [3,4) [4,5) [5,6) [0,8] (8,16] (16,24] (24,32] (32,40]
4.38 5.15 5.26 10.11 9.38 10.07 12.45 13.88 15.96 14.38 30.85 31.33 41.63 33.51 34.38 45.30 35.62 30.14 32.98 32.50 8.53 14.16 9.09 7.45 8.75 0.88 1.29 0.63 5