数学建模竞赛-艾滋病疗法评价及疗效预测模型_精品文档.pdf
- 文档编号:3216407
- 上传时间:2022-11-20
- 格式:PDF
- 页数:23
- 大小:688KB
数学建模竞赛-艾滋病疗法评价及疗效预测模型_精品文档.pdf
《数学建模竞赛-艾滋病疗法评价及疗效预测模型_精品文档.pdf》由会员分享,可在线阅读,更多相关《数学建模竞赛-艾滋病疗法评价及疗效预测模型_精品文档.pdf(23页珍藏版)》请在冰豆网上搜索。
12005年四川省数学建模竞赛艾滋病疗法评价及疗效预测模型张慧勇四川大学摘要:
本文研究了艾滋病的疗法评价及疗效预测的问题。
艾滋病(AIDS)是当前人类社会最严重的瘟疫之一,评价治疗艾滋病方法的优劣,预测继续治疗的效果对于治疗艾滋病有非常重要的意义。
问题1,通过分析、统计数据,做出CD4比值与时间、HIV数量与时间的散点图,根据散点图用三次多项式拟合CD4比值的曲线,用负指数函数拟合HIV的曲线;由HIV拟合曲线得出治疗10周,HIV就能被控制在稳定值2.84;已知患者开始时CD4数量Q,通过CD4比值曲线求出t时刻的比值f(t),就能预测出t时刻患者的CD4数量f(t)*Q。
问题2,本文建立了两个评价疗效模型:
多元回归模型、灰色关联度模型,两个模型都做出疗效4最优的评价。
本文还建立了三个预测模型:
多元回归模型、灰色预测模型、BP神经网络模型,并都对疗法4进行了预测,其中灰色预测模型能画出疗效与时间的关系曲线图,看出疗效4能确保患者体内CD4数量不减少,患者应继续治疗。
问题3,考虑4种疗法的费用,本文继续沿用问题2中的灰色关联度模型,但对其中的评价值矩阵H进行修改,对此引入了效益函数f=f(CD4,money,),表示患者的贫富程度,且f函数是的增函数;用此修改后的灰色关联的模型来评价4种疗法,大致结果是:
比较大时,疗法4最优;比较小时疗法1最优;本文在建立模型求解问题时,十分注重数据信息的挖掘,有的模型用CD4的比值,有的模型用CD4的差值,能很好地反映疗法对患者的疗效;问题三中的灰色关联度模型引入了自定义的效益函数,十分好地解决了评价带有费用的疗效的问题。
关键词:
曲线拟合、回归模型、灰色关联度、灰色预测、效益函数2一、问题的重述艾滋病(AIDS)是当前人类社会最严重的瘟疫之一,它是由艾滋病病毒(HIV)引起的。
这种病毒破坏人的免疫系统,使人体丧失抵抗各种疾病的能力,从而严重危害人的生命。
人类免疫系统的CD4细胞,作为一种识别因子,能帮助辅助T细胞识别其它细胞的感染信号,对于艾滋病病毒(HIV)感染来说,它是病毒进入和破坏细胞功能的重要的识别因子,当CD4被HIV感染而裂解时,其数量会急剧减少,HIV将迅速增加,导致AIDS发作,因此,在AIDS的治疗疗法评价及疗效预测中,CD4是一个很重要的指标。
AIDS治疗的目的,是尽量减少人体内HIV的数量,同时产生更多的CD4,至少要有效地降低CD4减少的速度,以提高人体免疫能力。
而迄今为止人类还没有找到能根治AIDS的疗法,目前的一些AIDS疗法不仅对人体有副作用,而且成本也很高,因此,对AIDS的治疗疗法评价及疗效预测是一个需要的问题。
解决的基本问题:
1、利用附件1的数据,预测继续治疗的效果,或者确定最佳治疗终止时间。
2、利用附件2的数据,评价4种疗法的优劣(仅以CD4为标准),并对较优的疗法预测继续治疗的效果,或者确定最佳治疗终止时间。
3、已知艾滋病药品的主要供给商对不发达国家提供的药品价格,如果病人需要考虑4种疗法的费用,对2中的评价和预测(或者提前终止)有什么改变。
二、基本假设1、题中所给出的数据真实,能反映现实情况;2、最佳治疗终止时间是指:
CD4的浓度(个/l)达到某值的初始时间;3、治疗有效果是指:
HIV的数量降低,同时产生更多的CD4细胞,至少能有效地降低CD4减少的速度;4、问题2中疗法的优劣评价,仅以CD4为标准;5、在疗法1中,以日用量600mgzidovudine或400mgdidanosine按每月交替服用,1个月按30天计算,且在第1月中,每日服用600mgzidovudine;6、在测试过程中,忽略病人年龄的变化;7、处于同一年龄段的病人,生理状况一致。
三、问题1的分析、建模与求解1、问题的分析通过分析统计数据,做散点图,拟合出疗效与时间的曲线,就能解决此问题,预测继续治疗的效果,或者确定最佳治疗终止时间。
2、建立数学模型31)对个别测试的CD4(和HIV)的时刻调整由于有个别测试CD4(和HIV)的开始时刻并不是第0周,故将其开始时刻调整为第0周。
即:
设某个体第一次测试的时刻为0t,若00t,则将其各次测试时刻均减去0t。
2)调整CD4的值为相对第0周测试时的比值我们认为服用药物后的CD4值的改变可以反映出疗效的好坏,由于所给出的CD4浓度值差异很大(最小值为0,最大值为9330.2/l个),故对其处理如下:
规定第0周CD4的调整值(4(0)CD)为1,第i周的调整值(4()CDi)为相应原始值(4iCD)与第0周原始值(04CD)的比值,即044()4iCDCDiCD。
3)统计统计第i周调整值4()CDi总和sumCD4(i)和测试人数numCD4(i),用三项移动平均法的思想,令4
(1)4()4
(1)4()4
(1)4()4
(1)sumCDisumCDisumCDimeanCDinumCDinumCDinumCDi用4()meanCDi的值代替4()sumCDi。
统计第i周HIV的测得值总和sumHIV(i)与测试人数numHIV(i),用三项移动平均法的思想,令
(1)()
(1)()
(1)()
(1)sumHIVisumHIVisumHIVimeanHIVinumHIVinumHIVinumHIVi用()meanHIVi的值代替()sumHIVi。
4)作图、拟合由于numCD4(i)=20时对应meanCD4(i)的散点图,根据图形做3次多项式拟合。
用MATLAB中的polyfit函数得到:
拟合多项式系数:
0.00030.02010.44140.7203a拟合曲线:
320.00030.02010.44140.7203yxxx相关系数的平方:
rsquare=0.9839残差平方和:
S=0.1540由此认为拟合比较好散点图、拟合曲线及置信区间如下图:
4图1:
CD4比值的曲线拟合图同样由于numHIV(i)=20时对应meanHIV(i)的散点图,根据散点图确定用负指数函数来拟合。
拟合曲线方程213kxykekgg用Matlab求得:
k1=3.4668,k2=0.4673,k3=2.83650.46733.46682.8365xye残差平方和:
S=0.0966由残差平方和与拟合曲线图看出拟合结果比较好图2:
HIV的曲线拟合图53、预测治疗效果由图1可知在0到45周治疗效果比较好,并能推测出继续治疗10周左右疗效也应比较好,但拟合曲线模型不能预测更长时间治疗效果。
由图2可知,治疗10周后HIV数量趋于稳定值k3;若某艾滋病患者测试开始时CD4数量为0Q(*0.2/l)第i周后CD4数量为:
32i0Q0.2(0.00030.02010.44140.7203)Qiii(/l)我们认为在HIV降低到稳定值k3的前提下,CD4浓度(个/l)上升到某值(根据文献2,CD4数量达到=350/l,就可以停药观察)的时间为最佳终止时间,因此解方程()iQi,求出i就得到最佳终止时间,即在第i周终止治疗。
4、模型评价及改进该模型建立思路比较简单,曲线拟合效果较好,在对数据处理中用了三项移动平均法,可以剔除数据中的偶然因素;该模型不能够预测较长时间的疗效;可以考虑通过某函数,将CD4浓度的改变量与HIV值结合在一起,进行综合考虑。
四、问题2的分析、建模与求解回归模型1、疗法评价1)对附件2的数据处理为了反映某种疗法的疗效,我们将附件2所给出的CD4值(为log(countCD4+1))进行处理。
规定第0周CD4变化值(04CD)为0,第i周的变化值(4iCD)为相应测得值(4iCD)与第0周测得值的差值,即0444iiCDCDCD。
2)评价4种疗法的优劣及较优疗法预测的多元回归模型令10jx采用第j种疗法否则,j=1,2,3,45x表示治疗时间(单位:
周),6x表示年龄(单位:
岁),y表示相应的疗效(即:
4iCD),建立多元回归方程:
2201122334455667586yxxxxxxxx通过求出1234,中的最大值,则对应的疗法最优。
3)疗法评价结果利用MATLAB函数regress(程序见附录)求得:
60.138240.0659900.129860.315160.002410.001410.000360.00001统计量F=46.02571(8,5027)1.94F(=0.05),知回归模型成立。
比较(1,2,3,4)ixi前系数i得,4321,由此可知对于所有年龄段的人,第4种疗法最优,第3种疗法次之,第1种疗法最差。
2、预测疗效用多元回归模型对最优疗法4进行预测,则方程为22055667586yxxxx(56,xx含义与以上一致)将采用疗法4病人的数据代入上述方程,用MATLAB函数regress求得:
1.16280.03050.05730.00090.0007且的置信区间都不包括0,知方程中y(疗效)与治疗时间及年龄的一次项、二次项都有线性关系。
统计量F=10.84651(4,1292)2.37F,知回归模型成立。
将病人的治疗时间及年龄代入上述多元回归方程,就可预测其采用疗法4的治疗效果。
3、模型评价及改进该回归模型能够对疗法优劣做出评价及疗效预测,但相关系数较低,可以通过对每一个点的残差置信区间的分析来剔除一些异常点,以此提高其回归的显著性。
灰色模型1、基于灰色关联度的方法评价4种疗法的优劣1)考虑年龄与不同疗法疗效的关系和各个年龄段统计量的多少,将年龄分段为:
12345age=2525age=3535age=4545age=5555age表1:
年龄段的划分2)数据处理为了反映疗效,规定第0周CD4的测得值(04CD)为0,第i周的变化量(4iCD)为相应测得值(4iCD)与第0周测得值的差值,即0444iiCDCDCD。
已知1300多名病人大约每隔8周测试的CD4浓度,为了数据处理方便,将第8(i-1)(i=1,2,3,4,5,6)周附近进行测试治疗的数据都归入第8(i-1)周。
3)对某一个具体的年龄段,选定母系列和比较序列7定义评价值矩阵64H,其中的元素ijh(i=1,2,3,4,5,6;j=1,2,3,4)表示该年龄段的采用第j种疗法病人在第8(i-1)周的8
(1)4iCD的平均值。
令014maxiijjhh,则母序列0ih(i=1,2,6)是构造出来的最理想评判值序列。
由1234,iiiihhhh分别构成4个子序列(比较序列)。
4)计算关联度子序列中每个数值与母序列中相应数值的关联系数为:
()()jjabviib,其中0()|jijiihh,1416minmin()jjiai,1416maxmax()jjibi,=0.5。
把关联系数取均值得到关联度:
611(),(1,2,3,4),6jjiuvij关联度表示疗法评判值序列和理想评判值序列相关联的程度,数值越大,表明疗法越好(对于该年龄段的人来说)。
5)模型的求解利用MATLAB编程计算得到各个年龄段的各关联度:
1234age=250.55770.66970.59570.958325age=350.50330.51470.60130.998835age=450.47450.58270.87000.850545age
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 竞赛 艾滋病 疗法 评价 疗效 预测 模型 精品 文档