2014年第三届数学建模国际赛小美赛
C题 埃博拉病毒:为什么不能遏制传播?
原题再现:
致命的埃博拉出血热最早发现于1976年,自从理查德·普雷斯顿的《热区》出版以来,它就一直困扰着公众20年。然而,在这段时间里,没有一种药物被批准用于治疗这种疾病。据世界卫生组织称,自2013年12月以来,西非持续爆发的埃博拉病毒已在几内亚、塞拉利昂和利比里亚感染至少567人,其中350人死亡。这次爆发似乎是历史上规模最大的一次,超过了2000年乌干达埃博拉爆发时发生的425例病例。埃博拉病毒感染者只接受旨在支持患者的一般治疗。根据疾病控制和预防中心的说法,他们可能需要补充液体(埃博拉患者经常脱水),或者进行旨在维持血压和氧气水平的治疗,如果出现感染,还可以进行治疗。(向两名美国患者提供的实验药物已经用完)埃博拉病毒治疗的缺乏令人不安。但是,考虑到药物开发的资金来源,这也是可以预测的。
当制药公司决定将研发资金投向何处时,他们自然会评估候选药物的潜在市场。这意味着他们有动机针对那些影响更富有的人(首先是发达国家的人)的疾病,这些人有能力支付很多费用。他们有动机去制造许多人会服用的药物。他们有一个动机,使药物,人们将定期服用很长一段时间的药物,如他汀类药物。
这个系统在让西方人得到他们想要的药物(尽管通常价格很高)方面做得很合理。但这也导致了对某些疾病和某些药物种类的巨大投资不足。主要影响穷国穷人的疾病并不是研究的重点,因为这些市场不太可能提供可观的回报。因此,像疟疾和肺结核这样每年夺去200万人生命的疾病,在制药公司受到的关注比高胆固醇要少。世界卫生组织称之为“被忽视的热带疾病”,如查加斯病和登革热;它们影响到10亿多人,每年造成50万人死亡。一项研究发现,在1975年至2004年间上市的1500多种药物中,只有10种是针对这些疾病的。当一种疾病的受害者既穷又不多的时候,这是双重打击。从这两个方面来看,一种治疗埃博拉病毒的药物看起来都是一项糟糕的投资:到目前为止,这种疾病只出现在贫穷国家,影响到的人数相对较少。
然而,这一体系损害的不仅仅是发展中国家。近年来,耐药微生物的增多使得我们使用的抗生素的效果降低,增加了传染病失控的风险。卫生官员同意,西方人需要的是我们可以保留的新药,以防常规抗生素无法控制的疾病爆发。然而,在过去的三十年里,新抗生素的供应已经放缓到涓涓细流。波士顿大学卫生法项目的联合主任凯文·奥特森(Kevin Outterson)是美国疾病预防控制中心(C.D.C.)抗生素耐药性工作组的创始成员,他告诉我:“抗生素耐药性确实有可能使我们生活方式的方方面面发生变化。因此,我们需要推动这条管道。”
问题还是在于商业模式。如果一家制药公司真的发明了一种功能强大的新抗生素,我们就不会希望它被广泛使用,因为它的目标是延缓耐药性。奥特森说:“公共卫生官员会尽可能适当地限制这种药物的销售。”;投资前景不佳。
那么,我们如何才能获得我们所需要的药物,而不神奇地改变开发药物的行业呢?关键是奖励那些创造了大量公共卫生福利的公司。最简单的方法就是为新药提供奖励。Outterson描述了一个场景:“政府将向公司支付一笔或一系列款项,作为交换,公司将放弃销售该产品的权利。”制药公司将获得报酬,并将避免尝试推出新产品的所有费用(无论如何,你不希望使用最后手段的抗生素)。社会将得到一种新药,公共卫生官员将能够控制它的推广和使用方式。
任务#1:选择您认为在非洲埃博拉疫情中最严重的国家。建立一个模型,在没有任何额外药物的情况下,估算2006年至2025年该国埃博拉感染人数的预期变化率。充分解释您的模型和模型的基础假设。此外,估计我们应对危机的成本。
任务#2:建立一个数学模型,确定在此期间该地区是否有任何新疾病需要开发药物。如果是,估计药物开发过程的截止日期和产量。
任务#3:制定制药企业“奖励”政策,支持药品研发工作。您制定的方法应利用公共卫生福利和美国预算。
任务#4:多公司联合开展药物研究,可有效降低成本和风险。如何选择一些企业建立合作关系是一个棘手的问题。您的团队可以帮助政府做出选择吗?
整体求解过程概述(摘要)
针对这一问题,我们从各个方面对埃博拉病毒进行了研究。我们查阅了互联网上有关埃博拉病毒在非洲爆发的数据。比较不同国家的病例、死亡率等影响最严重的因素。找到问题的内部信息和特征对我们来说不是很明显,也就是说,有一部分是未知的,然后找到不同的数学模型来求解,利用收集到的数据绘制出速率曲线,并预测其未来的变化。
围绕任务1,我们从不同方面找到了美国成本的财务状况。预测埃博拉病毒治疗药物的经济发展。并预测美国将花费多少来应对这场危机。收集有关非洲病毒的信息。联系经济、病毒等诸多因素,寻找合适的国家。利用灰色预测模型和GM(1,1)模型在MATLAB中进行编程,绘制各因素的曲线。
任务#2是多项式回归曲线的拟合,我们用Logistic回归模型进行求解,将不同的数据进行拟合。通过这两个模型,我们预测了药物开发的可能性。最后对埃博拉病毒检测结果进行了分析。研究药物研究的影响因素。发现主要影响来源是病毒的发育速度、死亡率、变化的脱色率。将三个因素拟合到变化曲线中。
针对任务#3,建立了“奖励”策略模型,描述了“奖励”策略各方面的客观性。在此基础上,提出了创新的“奖励”政策。激励公司创造更大的社会效益。完善“奖励”政策。
在任务4中,按照各公司的方案处理药品开发风险和不同成本。通过任务分配模型,降低药物开发的成本和风险。优化作业,提高工作效率。
模型假设:
(1) 在没有其他新的病毒感染的情况下,我们只讨论埃博拉病毒的发展。
(2) 摆脱环境和温度的巨大变化。
(3) 确保病毒不会有很多变态。然后,它保持其主要的增长曲线。
(4) 该新药尚未开发,因此不会影响埃博拉病毒的传播速度。
(5) 不同的病毒有不同的发展速度。在我们的模型中,我们关注的是会被感染的病毒。
(6) 我们的模型不考虑非感染病毒。当我们确认研究的合理性时,我们可以在合理的时间内完成。
(7) 我们研究的人口问题是根据人口的平均密度。
问题分析:
问题1分析
我们需要选择一个非洲国家作为埃博拉病毒感染的主要研究目标。根据世界卫生组织的数据,我们关注乌干达。根据2000年、2007年、2011年、2012年和2013年埃博拉病毒感染的现有数据,通过线性插值得到2000年至2014年的相关数据。然后利用灰色预测方法对2014-2025年的感染风险人数进行估计,最后计算出感染人数的变化率。对于美国对这种病毒的财政援助,我们的目标是医药投资的成本。根据相关文献提供的药品不同阶段的机会成本,我们以现金和资本两种形式评估支出。结合美国的贴现率和成本消耗的增长,我们可以计算出2010年每种药物的平均成本。
问题2分析
这个问题要求我们确定2006年至2025年乌干达是否需要新药。我们选择该区域更重要的疾病类型:结核病、艾滋病和疟疾。如果方差检验有显著性差异,则新药的需求是明显的。在确定新药需求后,我们提取每种疾病每年的异常值感染数。药品的总产量是根据一定时期内感染的确切人数来估计的。
问题3分析
这个问题要求我们设计一个方案,在药物开发过程中为制药公司提供资金支持。该方案旨在降低药品价格,压缩药品销售,使非洲较贫穷的国家能够负担费用,使药品不会因广泛使用而降低药效。我们将对制药企业的补贴分为两个阶段:药品研发阶段和销售阶段。政府补贴来自政府自己的预算和公共卫生福利。为了便于计算,我们将公共卫生福利资本化,从而获得最佳的援助方案。
问题4的分析
第四个问题要求政府设计一个选择公司合作药品的方案,以降低成本和风险。为了有效发挥不同企业的优势,发挥其在医药行业的优势,根据公司的相关数据,我们定义了四个指标来评价其药物研发能力,根据各公司的各种参数值,选择各公司的最大优势进行混搭。最后,对于可持续性的考虑,我们考虑了公司数量增加所带来的影响。
模型的建立与求解整体论文缩略图
全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可
部分程序代码:(代码和文档not free)
clc,clear
Number=input('please enter the codes of problem(1 min risks 2lowest cost)Number=:');
if Number==1
F=[5 4 5 1 2 3 2 4 5 5 1 5 1 3 1 2 3 2 4 2 1 3 4 2 3];
elseif Number==2
F=[20 19 20 28 17 18 24 27 20 20 26 16 15 18 15 17 20 24 19 16 15 18 21 17 21];
end
A=[];B=[];
Aeq=[ones(1,5) zeros(1,20);zeros(1,5) ones(1,5) zeros(1,15);zeros(1,10) ones(1,5)
zeros(1,10);...
zeros(1,15) ones(1,5) zeros(1,5);zeros(1,20) ones(1,5);...
1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4);...
zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);...
zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);...
zeros(1,3) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);...
zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
Beq=ones(1,10);
if Number==1
[x,min_fval]=bintprog(F,A,B,Aeq,Beq);
min_fval2=min_fval/5;
fprintf('(1)the match function of company:\n')
fprintf(' (A1) the match function of company A(1 replace match 0 replace not match):
\n')
disp(' A-1 A-2 A-3 A-4 A-5'),disp(x(1:5)')
fprintf(' (B1) the match function of company B(1 replace match 0 replace not match):
\n')
disp(' B-1 B-2 B-3 B-4 B-5'),disp(x(6:10)')
fprintf(' (C1) the match function of company C(1 replace match 0 replace not match):
\n')
disp(' C-1 C-2 C-3 C-4 C-5'),disp(x(11:15)')
fprintf(' (D1) the match function of company D(1 replace match 0 replace not match):
\n')
disp(' D-1 D-2 D-3 D-4 D-5'),disp(x(16:20)')
fprintf(' (E1) the match function of company E(1 replace match 0 replace not match):
\n')
disp(' E-1 E-2 E-3 E-4 E-5'),disp(x(21:25)')
fprintf('min risks:'),disp(min_fval2)
elseif Number==2
[x,min_fval]=bintprog(F,A,B,Aeq,Beq);
min_fval2=min_fval/5;
fprintf('(2)the match function of company:\n')
fprintf(' (A2) the match function of company A(1 replace match 0 replace not match):
\n')
disp(' A-1 A-2 A-3 A-4 A-5'),disp(x(1:5)')
fprintf(' (B2) the match function of company B(1 replace match 0 replace not match):
\n')
disp(' B-1 B-2 B-3 B-4 B-5'),disp(x(6:10)')
fprintf(' (C2) the match function of company C(1 replace match 0 replace not match):
\n')
disp(' C-1 C-2 C-3 C-4 C-5'),disp(x(11:15)')
fprintf(' (D2) the match function of company D(1 replace match 0 replace not match):
\n')
disp(' D-1 D-2 D-3 D-4 D-5'),disp(x(16:20)')
fprintf(' (E2) the match function of company E(1 replace match 0 replace not match):
\n')
disp(' E-1 E-2 E-3 E-4 E-5'),disp(x(21:25)')
fprintf('lowest costs:'),disp(min_fval2)
end
%The result of these as follows:
Please enter the codes of problem(1min risks 2 lowest cost)Number=:1
(1)the match function of company:
(A1) the match function of company A(1 replace match 0 replace not match):
A-1 A-2 A-3 A-4 A-5
0 0 0 1 0
(B1) the match function of company B(1 replace match 0 replace not match):
B-1 B-2 B-3 B-4 B-5
0 1 0 0 0
(C1) the match function of company C(1 replace match 0 replace not match):
C-1 C-2 C-3 C-4 C-5
0 0 0 0 1
(D1) the match function of company D(1 replace match 0 replace not match):
D-1 D-2 D-3 D-4 D-5
0 0 1 0 0
(E1) the match function of company E(1 replace match 0 replace not match):
E-1 E-2 E-3 E-4 E-5
1 0 0 0 0
min risks: 1.4000
Please enter the codes of problem(1min risks 2 lowest cost)Number=:2
(1)the match function of company:
(A2) the match function of company A(1 replace match 0 replace not match):
A-1 A-2 A-3 A-4 A-5
0 1 0 0 0
(B2) the match function of company B(1 replace match 0 replace not match):
B-1 B-2 B-3 B-4 B-5
0 0 0 1 0
(C2) the match function of company C(1 replace match 0 replace not match):
C-1 C-2 C-3 C-4 C-5
0 0 1 0 0
(D2) the match function of company D(1 replace match 0 replace not match):
D-1 D-2 D-3 D-4 D-5
0 0 0 0 1
(E2) the match function of company E(1 replace match 0 replace not match):
E-1 E-2 E-3 E-4 E-5
1 0 0 0 0
lowest costs: 17