数学与应用数学级班张云乾.docVIP

题 目 基于概率和预测模型的流感疫苗投放问题 关 键 词 概率模型 贝叶斯预测模型 评估体系 投放 摘 要： 问题1：我们设想将三联装疫苗改成双联装疫苗，并不是在所有的区域内实行，根据题目给与的数据，根据概率模型找出适合投放双联装疫苗的区域，而其他区域依旧投放三联装疫苗，从而达到降低成本的目的。 现行的投放方案是三联装疫苗的投放形式，我们通过基于因子分析定权评价模型来判评双联装疫苗的优越性。 问题2：通过问题1中的数据统计，依照贝叶斯预测模型，对不同区域检测到的流感记录进行整合处理，筛选出次年爆发可能性最高的病毒类型，从而缩小病毒专家的研究范围，便于疫苗成分的确定和生产。  目 录 一、问题重述 2 二、模型假设 2 三、符号说明 2 四、问题的分析 2 五、模型的建立与求解 3 1．划分小区域并建立概率模型 3 2．关于投放疫苗的评估标准和效果 7 3．利用评估体系，分析三联与两联的对比评价 8 4．从检测的流感记录中取一部分样本 10 5．建立贝叶斯预测模型 10 6. 筛选下一年的新病毒疫苗 11 六、模型的评价与改进 12 参考文献 13 附录： 14 问题重述 流感是一种广泛流行于世界范围内的疾病，每次流感大流行都会造成多人死亡和巨大损失。世界卫生组织大力推荐将疫苗作为一种有效的预防措施来抗击这种潜在的致命性疾病。依其两种主要抗原 (HA，NA) 的不同，区分为不同的亚型。B 型病毒的变异比较缓慢，C 型病毒甚少对人类造成威胁。也就是每份疫苗中有三种经过灭活或裂解处理3的毒株，分别为两个 A 型和一个 B 型。将北半球和南半球分别划分为稍小的区域，并使用不同的疫苗针对不同区域来进行投放表示发病的概率； DLM{1，1，V，δ} 每一时刻t常均值折模型 δ 发病因子 MAPE对预测结果的准确度采用平均绝对百分误差分析 四、 问题分析 问题1中只考虑病毒亚型，不考虑具体的毒株选择，应考虑以下四个问题： 1.建立相应的预测模型； 2.需要划分南北半球为小区域的问题； 3.关于二联用药后的评估标准和效果； 4.三联与两联的对比，利用因子分析定权的评估体系； 问题2中考虑选择具体的毒株，应考虑以下三个问题： 1.从检测的流感流行记录中取一部分样本数据； 2.由这些数据建立一个类似于贝叶斯预测模型； 3.利用模型筛选下一年威胁较大的病毒可能在哪里出现； 五、模型的建立与求解 1．取划分后南北半球中几个小区域并进行建立概率模型 根据题目要求，我们可以从题目提供的数据中分离出流感病毒A型的两个亚型H1N1和H3N2以及B型毒株的分布区域，以及注射带有A型亚型H1N1和H3N2及B型灭活或裂解毒株疫苗对不同人群的效果情况，此效果是有注射疫苗后在一定时间后人体内抗原的减少情况来设定的。 由题意知，我们可以只考虑病毒亚型，而不考虑具体毒株的选择。 疫苗中可能含有的毒株成分： 三联装疫苗由Ａ型的H1N1和H3N2亚型，以及一个Ｂ型灭活或裂解的毒株组成。 双联装疫苗有三种组成： 一是，Ａ型的H1N1亚型和Ｂ型灭活或裂解的毒株组成。　　　　　　　　 二是，Ａ型的H3N2亚型和Ｂ型灭活或裂解的毒株组成。 　　　　　　　　　　　　三是，由Ａ型的H1N1和H3N2亚型的两个灭活或裂解的毒株组成。 对2007年九月到2009年一月数据（见附录1）分大洲选取2-3个国家按照发病次数以及发病范围分析得到下表： H1 H3 B 中国 15 18 18 菲律宾 10 9 14 蒙古 2 3 3 智利 15 1 10 埃及 3 3 9 马达加斯加 18 3 6 保加利亚 3 1 0 罗马尼亚 5 3 5 澳大利亚 10 17 19 新西兰 3 19 20 说明：“*”代表小范围有感染，用1表示；“**”表示区域爆发，用2表示； “***”代表较大规模爆发，用3表示；“****”大规模爆发，用4表示；表中数据为爆发强度加和。 根据数据得出柱状图如图所示： 根据图表中数据若这些地区使用双联装疫苗对发病率具体分析如下； 中国： 图表数据容易看出最合适的双联装为Ｈ3+B类的疫苗，则其中会出现一部分投放不合适的情况，此部分占 = 此部分人群一定会发病，还有 = 的人群可能会发病。则中国地区若投放双联型疫苗总的发病率为 = +=43.530% 菲律宾： 同样从图片中可以看出对菲律宾若使用双联装疫苗则H1+B类的比较合适， 则投放不合适的情况 =27.273% 投放合适的人群为 =69.697% 总发病率为 =+=96.970% 蒙古： 从图表中可以看出最好的双联装为H3+B类疫苗 则投放不合适的情况为 =25.000% 投放合适的人