艾滋病疗法的评价及疗效的预测但数学建模竞赛优秀.doc
《艾滋病疗法的评价及疗效的预测但数学建模竞赛优秀.doc》由会员分享,可在线阅读,更多相关《艾滋病疗法的评价及疗效的预测但数学建模竞赛优秀.doc(26页珍藏版)》请在三一文库上搜索。
1、高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。我们参赛的题目是: 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 参
2、赛队员 (打印并签名) :1. 2. 3. 指导教师或指导教师组负责人 (打印并签名): 日期: 年 月 日赛区评阅编号(由赛区组委会评阅前进行编号):1高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):艾滋病疗法的评价及疗效的预测摘要针对问题一我们建立了两个模型:曲线拟合模型和灰色预测模型。对曲线拟合模型,我们采用了分区间分类讨论的思想,使预测精度有了很大提高;在建立灰色预测模型的的过程中,我们对经典的GM(1,1
3、)模型进行了拓展,建立了不等时距的拓GM(1,1)模型,并对其算法进行了简化。最终得出结果:在未来一段时间内CD4浓度将持续增长,HIV指标将缓慢低速下降,应继续服药。问题二可以分解为两个子问题:对四种疗法的评价和对最优疗法的预测。我们的总体思想是:先将年龄分段,再对每个年龄段进行评价和预测。在评价中,我们将疗效定义为:期望与方差的加权组合。最终结果为:适合青、中、老年的最优疗法分别为:疗法四、二、一。在预测中,我们又引进了指数平滑模型,对不同年龄段对应的最优疗法作出了合理的预测,并进确定出个疗法的最佳治疗终止时间为:青年试验结束后135.68周;中年试验开始后816周;老年试验结束后0.47
4、2周。问题三是对问题二的进一步延伸,我们在问题二的评价基础上,引入花费效应,利用双指标评价模型,得到了不同权重下各疗法的优劣。通过对数据的趋势分析,结合实际情况,引入了比较可靠的权重作为评价参数,得出了各年龄段的最优疗法。另外我们在确定最佳治疗终止时间时,通对过金钱与疗效之间相关系数k的引入,将问题简化为单目标优化问题。在对数据的处理上,我们采用了化区间为点的思想,提高了计算效率;在模型的进一步讨论中,我们对模型的遗留问题与拓展方向进行了研究,使之更具有普遍性。 关键字 曲线拟合 灰色预测 不等时距的GM(1,1)模型 权变理论 指数平滑 双指标评价 问题重述艾滋病是当前人类社会最严重的瘟疫之
5、一,它是由艾滋病毒(医学全名为“人体免疫缺损病毒”, 英文简称HIV)引起的。人类免疫系统的CD4细胞在抵御HIV的入侵中起着重要作用,当CD4被HIV感染而裂解时,其数量会急剧减少,HIV将迅速增加,导致AIDS发作。 艾滋病治疗的目的,是尽量减少人体内HIV的数量,同时产生更多的CD4,至少要有效地降低CD4减少的速度,以提高人体免疫能力。现在得到了美国艾滋病医疗试验机构ACTG公布的两组数据。 ACTG320(见附件1)是同时服用zidovudine(齐多夫定),lamivudine(拉美夫定)和indinavir(茚地那韦)3种药物的300多名病人每隔几周测试的CD4和HIV的浓度(每
6、毫升血液里的数量)。193A(见附件2)是将1300多名病人随机地分为4组,每组按下述4种疗法中的一种服药,大约每隔8周测试的CD4浓度(这组数据缺HIV浓度,它的测试成本很高)。4种疗法的日用药分别为:600mg zidovudine或400mg didanosine(去羟基苷),这两种药按月轮换使用;600 mg zidovudine加2.25 mg zalcitabine(扎西他滨);600 mg zidovudine加400 mg didanosine;600 mg zidovudine加400 mg didanosine,再加400 mg nevirapine(奈韦拉平)。问题:(1
7、)利用附件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)中的评价和预测(或者提前终
8、止)有什么改变。条件假设与符号说明(一) 条件假设(1) 附件中的数据都是真实、可靠的;(2) 在测试期间,患者接受治疗是连续的;(3) 只考虑药物的正面影响,不考虑副作用;(4) CD4与HIV是相对独立的,作为两个指标;(5) 在测试时刻间隙,患者的病情无急剧变化现象;(6) 参与试验的患者相对于所有艾滋病患者是随机的。(二) 符号说明 Pi: 第i种疗法 yi: 第i个年龄区间 Ti: 第i个测试时间段 Ci: 第i个测试时间段CD4 的累加值 Hi: 第i个测试段的HIVRi: 第i种疗法的综合效应Ti: 第i种疗法的最佳停止时间ci: 第i个测试段的CD4 hi: 第i个测试时间段H
9、IV si: 第i个年龄段的药效的方差 li: 第i个疗法的加权效应(不考虑花费)pi: 第i种疗法的平均治疗费用 ni 第i周检测者的个数问题分析问题一是一个典型的预测问题,我们应在对已知数据进行充分分析的基础上,建立疗效与时间的关系,进行预测。首先,由于疗效的评价标准有两个方面:CD4及HIV。由前面的假设,这两者是相互独立的,这就从总体上讲该问题分解为两个子问题。题目给出的数据量庞大的且看起来似乎无严格的规律,因此必须对数据进行预处理,可以采用加权平均或累积求和的方法进行数据的重生成。问题二也可以分解为两个子问题,即对不同的疗法的综合评价和对较优疗法的预测。对于疗法的评价,需要建立一套完
10、整的评价理论体系,我们考虑用累积效应(期望)稳定性(方差)两个指标作为评价标准,最后需要重新定义一个统一的标准即为这两者的合成。在具体建模中,年龄与时间都是要考虑的因素,我们应对目标分层分析,固定一个来考虑另一个的变化。 问题三是对问题二的扩充,在问题二的基础上又增加了费用这一指标。该问题也是两个子问题,首先对于疗法的重新评价,可以转化为双目标评价问题,通过权重的引入将问题简化为单目标。对于最佳治疗终止时间的确定,我们考虑找出花费的金钱与取得的疗效的相关关系,并用简单的函数表示出来进而对疗效预测与优化。模型的建立与求解问题一该问题是一个典型的预测问题,根据需要,我们用以下两种方法建立模型并将得
11、到的结果进行比较:模型一:曲线拟合模型根据对数据的观察,CD4及HIV随时间有比较明显的变化关系,我们首先考虑运用比较简单的曲线拟合模型进行预测。首先对数据进行预处理:我们使题目所给的离散数据有序化,通过对题目所给数据的特点以及355名患者服药后状况的综合考虑,我们对数据进行以下处理:Step1:对不同周次检测的所有患者的CD4和HIV取平均值;将第i周患者CD4和HIV的平均值分别记为和。Step2:记患者服药后药效在第i周前对患者CD4的积累效应为;则: Step3:将时间(周)分为:、和三段。这样我们就得到了患者服药后不同时刻(周)的HIV平均值和药效对CD4的积累效应。如下表所示:表1
12、1 积累效应随时刻变化表t34567891011Ci1.60541.55621.5211.44281.69861.76892.01911.83171.0924HIV2.95663.22642.99183.15453.08032.88562.46562.653.5ni5323161116120864104t212223242526272829Ci1.65461.96122.36652.09421.97571.88941.43321.6081.4612HIV2.28333.06252.72652 .76812.4443.03082.53754.23.94ni616641385926835t3637
13、38394041424344Ci0.64094.55891.82942.48952.28372.03711.0371.51671.5614HIV03.442.76672.51152.42232.56772.28752.93331.6875ni16155294341665根据以上表格中的数据用Matlab可以绘出与的折线图:图一 a b c a b c() ()上面的折线图规律性较差,对此我们采用以下两种方案:方案一:剔除异常点近似曲线拟合 由表11可知:第11,28,29,36,44周的检测者较少,这几周的数据没有普遍性,不能真正反映药物的药效。因此我们剔除这几组数据,然后再用Matlab得到
14、下面的折线图:图二 a b c a b c() ()图二()-a近似的是一条三次曲线;设其方程为:用Matlab求得该方程为: 如图三()-a所示。用同样的方法可以得出其他五组数据对应的方程: (图三()-b) (图三()-c) (图三()-a) (图三()-b) (图三()-c)图三 a b c a b c() ()方案二:变量累加后曲线拟合虽然方案一剔除了异常数据,但有的拟合曲线依然很粗糙。在此我们考虑:将CD4和HIV累加几次后再拟合曲线,得到其累加值与时间的函数关系。累加即令:; (1)其中i表示某区间内的第i周。结果累加一次后得到表12的数据。表12t34567891011Ci1.6
15、0543.16164.68266.12547.8249.592911.61213.44414.536HIV2.95666.1839.174812.32915.4118.29520.76123.41126.911t212223242526272829Ci1.65463.61585.98238.076510.05211.94213.37514.98316.444HIV2.28335.34588.072310.8413.28416.31518.85323.05326.993t363738394041424344Ci0.64095.19987.02929.518711.80213.8414.87716
16、.39317.955HIV03.446.20678.718211.14113.70815.99618.92920.617由表12易得出CD4和HIV的累加值与时间近似的是线性关系。用Matlab拟合出各时间段CD4和HIV的累加值与时间的图像如图四所示:图四 a b c a b c 由图四容易看出CD4和HIV的累加值与时间是较为严格的线性关系。用Matlab得到函数关系为: 由(1)式可知: 因此,由(2)式和(3)可得: (4) (5)(4)式和(5)式反映的是各个区间内CD4和HIV的平均水平。由方案一中图三可以看出CD4和HIV随时间上下波动。也就是说该药的药效不稳定。记CD4的第一个
17、极小值点的函数值为:;HIV的第一个极大值点的函数值为:。由图三可以看出在大部分时间段都有:;。由此可见:患者同时服用zidovudine,lamivudine和indinavir总体上效果是比较明显的。由(1)式可知:按我们的累加规则得到的。故由(4)式可知:服药后总体上患者的CD4在大幅度持续增长。这说明zidovudine,lamivudine和indinavir在提高患者免疫力方面有较好的效果。 由(5)式可知总体上患者的HIV浓度先上升后下降,这与方案一中得到的药效的不稳定性恰好吻合。同时也反映出其药效慢的特点。尽管由(5)式我们可以得知在后一段时间内总体上患者的HIV浓度在下降,但
18、是图三中最后一段HIV浓度呈上升趋势。即患者在服药四十几周HIV浓度会在一段时期内呈上升趋势。这说明该种药物组合在总体上可以控制HIV浓度但在短时间内波动性较强,稳定性较差。综上所述:同时服用zidovudine,lamivudine和indinavir3种药物在总体上可以提高患者免疫力(即提高患者的CD4值)可以控制HIV的浓度。但在短时期内病情会出现反弹现象;而且HIV浓度降低的速度较慢。但总体上有一定的效果,应继续服用该种药物组合。模型二:不等时距的GM(1,1)模型由上面的分析可以知道,这是一个只有部分数据的预测问题。从个体的角度来考虑,每个人的数据是较少的,因而我们考虑运用灰色预测系
19、统进行建模求解。我们对问题的重新讨论及简化处理如下:1、对时间的简化。在该问题中,每个病人的测试时间是不相同的,这给我们对个体进行统计、综合处理带来了麻烦。但经过对数据的研究,我们找出了测试时间的几个集中点,即0、4、8、24、40、56,我们分别记为t1、t2、t3、t4、t5、t6。我们把符合(其中的值根据不同的时间集中点适度取值,一般不超过,i=1,2,3,4,5,6,j=j|0j60,jZ的时间tj对应的数据均视为ti对应的数据。由于我们是运用加权平均的思想和化区间为点的方法,这种简化对预测趋势和精度影响不大。因此我们将问题转化为通过加权平均后的ti对应的数据,可以进行一定区间内的灰度
20、预测。2、对数据的处理。在上面的模型中我们采用了加权平均的方法,在本题中,由于灰度预测适合较少数据的情况,我们仍然采用加权平均,我们设tj时刻第k个人CD4及HIV分别为Cj,k、Hj,k,则有j= j= (n为所有试验的人数) 3、对时距的处理。在GM(1,1)基本模型中,要求做一次累加生成。当k只是一个不等时距的时间序列时,这种累加无法实现。我们考虑对原有的等时距GM(1,1)模型进行拓展。我们不妨假设等时距的原始数据是客观存在的,只是由于某些原因是其中的一些数据缺失,因而出现了不等时距的原始数列。进一步假设,我们已得到了与这些原始数据较为符合的GM(1,1)模型曲线,曲线的离散形式为:式
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 艾滋病 疗法 评价 疗效 预测 数学 建模 竞赛 优秀
链接地址:https://www.31doc.com/p-2175997.html