来源：李迅雷金融与投资（ID： lixunlei0722）文章已获授权前言截止到2020年2月2日24:00，我国确诊的新型冠状病毒患者已超过1.4万人，短短一个来月时间确诊病例数已经远超2003年“非典”疫情的全部确诊数，目前每日新增的确诊数仍在攀升，疫情的传播速度超过了大多数人的预期。科学客观地评估新冠肺炎的传染性强弱以及预测患病人数规模和峰值时间，对决策者实施必要的防控措施、评估对经济的影响以及投资者如何应对都具有重要的现实意义。作为数学建模的爱好者，而非流行病学的专业人员，作者通过搜集和学习海内外文献，对经典的流行病传播模型（SEIR）有了较准确的理解，通过适当的简化，较直观地向读者介绍模型的原理和各参数的意义。对于经典且成熟的模型来说，输出结果的可靠性完全由输入参数的准确性决定，因此本文从新冠病毒的具体情况入手，着重从多个角度来评估、检验和校正模型的输入值，提高模型输出结果的可靠性。主要内容&结论
其中S(t)、E(t)、I(t)和R(t)分别表示t时刻社群中处于易感染态、潜伏期、传染期和隔离态个体数目，N表示社群中个体总数目，且N=S(t)+E(t)+I(t)+R(t)。据参考文献[3],基本传染数可表示为：
来源：百度地图，各地卫健委简化微分方程组，建立数值模拟模型 逐一来理解四个微分方程：
b.在1月14日至21日期间潜伏期人群的平均人数在5257左右。两个待定参数，两个约束条件，理论上是可以求解的。4）模拟优化过程：参数R0、TI各自模拟区间内随机取值，通过3个迭代方程可以计算出E(t)、I(t)、R(t)的时间序列（以2019年12月8日为起点出现首个确诊患者），通过时间序列可以计算λ以及1月14日至21日期间潜伏期人数，判断约束条件是否接近。大量模拟后取约束条件最接近的参数值，作为近似解。 
其中k是一个有传染能力的患者平均每天与易感人群的接触次数，b是每次接触传染成功的概率，D是可以传播的时间，可以认为等同于TI。 新冠病毒R0明显高于SARS主要是因为平均每个病人传播的时间更长 比较新型冠状病毒和SARS，考虑到新冠病毒传播初期刚好遇到春运，人员流动性高，假设单位时间接触次数比SARS时期高10%。图表5：新型冠状病毒与SARS对比
可见，新冠病毒明显高于SARS主要是因为平均每个病人的传染期时间更长，而单次接触的传染概率可能略低于SARS。新冠病毒传染期格外长可能主要是因为其有大量的轻症患者，不会及时就医或采取措施，而会较长时间处于活跃的传染状态。国家卫健委的专家：“说它是新型冠状病毒肺炎，不如说是新型冠状病毒感染，因为检测阳性的患者不见得一定出现肺炎。就像流感、重感冒，有的轻症患者不需要住院便能治愈。病例数虽多，但大多数都是轻症。”而恰恰是轻症患者较多增加了不小的防控难度。对防控措施应该有信心采取防控措施就是要降低k、b、D这几个参数值，从而降低传染数。虽然新增确诊数还在不断攀升，但从1月23日之后采取的一系列措施来看，我们应该对控制疫情有信心。在流行病传播模型里，只要将基本传染数R0降到1以下，疫情就不会失控。 假定R0就如前文测算的5.38，远高于SARS。首先，交通管制、取消聚会等措施使得人员流动大幅下降，单位时间的接触次数k下降70%甚至更多都有可能。其次，通过宣传和提醒使得民众对疫情具有极高的警惕性，有轻微症状的患者大多也会尽快就医，从而减少传染期的长度，假设传染期的长度由7天降为4天，则D下降43%。第三，通过佩戴口罩、消毒等防护措施，单次接触传染的概率假如能下降30%。则传染数就下降为5.38*（1-70%）（1-43%）（1-30%）=0.64，明显小于1，疫情就能较快控制住了。事实上，WHO专家曾表示，中国在2003年非典期间后期，严格的防疫措施使得SARS的基本传染数下降到0.4。防控措施的效果体尚未体现到新增确诊数 为什么采取了这么严格的防控措施，新增确诊数还在不断攀升呢？首先，可能由于确诊检测试剂盒上市较晚，一些地方的基层医院缺乏确诊的条件或者确诊效率很低，造成确诊数据失真。可能存在前些天发病的病例过了几天才有条件确诊，这就可能造成近日新增确诊数大幅增加。其次，1月24日采取了进出武汉的交通管制，之前大量的人口从武汉外流到各地，按照10天左右的潜伏期，现在是发病的高峰期，随着他们的确诊，跟他们密切接触者大多会进入隔离观察，隔离观察者中还是会有比较多的潜伏期感染者，因此接下来新增确诊数可能还会保持高位，直到这批隔离观察者的潜伏期过去，被确诊，之后的新增确诊数就会大幅下降。从时间上看，可能要到2月10日之后。 以上海为例，其近日公布的新增确诊病例中，大部分还是曾湖北居住或旅游，说明目前新增的确诊病例主要还是原来潜伏期的感染者，严格的管控措施产生的效果还没有反映到新增确诊数据上来。图表6：上海每日新增确诊病例中有湖北居住或旅行史的占比
数据来源：Wind，WHO图表11：SARS期间美股、港股、A股表现
主要内容：
1）基于SEIR流行病传播模型测算基本传染数R0（也叫基本再生数）来衡量传染性的强弱普遍存在输入参数不易确定、输出结果对参数非常敏感的问题。国内外不同文献对同一疫情（比如SARS）的测算结果可能相差甚远。简单取一组参数得到的输出结果并不可靠。
2）本文对SEIR模型做适当简化，减少参数，再针对新冠病毒已明确的具体信息设定参数区间范围，再从多个维度选取校验条件，利用计算机数值模拟的方法求取参数值。通俗的来讲就是让计算机程序在参数区间范围内做大量尝试，试探出较合理的参数值，使得各个校验条件都能较好的相恰。
3）通过1月14日到22日武汉迁出人口排名前8的省市随后一周的确诊病例数与武汉迁出人数的关系，来反推武汉潜伏期感染者的比例和数目，获得求解模型的重要约束条件。通过数值求解得到参数值，代入模型对新冠疫情进行评估和预判。
主要结论：
1）新冠病毒基本传染数明显大于SARS，但单次接触的传染概率可能略低，民众无需恐慌。根据数值模拟确定的参数值，新冠病毒的基本传染数（R0）远高于SARS，即不采取防疫措施的情况下平均每个病人传染的人数更多。但进一步分析表明，新冠病毒R0更高的主要原因可能是潜伏期末期具有传染性和轻症患者比较多，潜伏期末期和轻症患者由于症状不明显，不会及时就医住院或限制自己活动，使得暴露在易感人群中的时间（传染期）更长。而单次接触的传染概率则可能略低于SARS。可见，潜伏期长和轻症患者多增加了疫情防控的难度，对病例的接触者采取更早、更全、更严格的隔离可能是最有效的防控措施。
2）当前新增确诊病例数大幅攀升（日新增病例超2500）可能与检测确诊效率提高有关，预计将很快回落（模型中性假设下预测值是1700左右）。1月25日之前，公布新增确诊数大幅低于模型预测值，26日之后则大幅高于模型预测值，可能是因为确诊检测试剂盒供应充分了，将之前已经发病而没有确诊的病例确诊了。
3）假定隔离和防控措施能达到或接近SARS后期的效果，则预计未来总患病人数将达到88500人（中性假设下）或58000（乐观假设下）左右，新增确诊人数将在2月17日（中性假设下）或2月11日（乐观假设下）左右出现拐点。
4）新冠病毒的重症及致死率低于SARS、禽流感、H1N1等，由于潜伏期长且轻症患者较多，控制疫情传播需要更加严格的隔离防控措施，对短期经济会产生明显影响。
5）从2009年美国H1N1大爆发以及2003年“非典”期间美股、港股和A股表现来看，在疫情快速传播期，股市会明显调整，但幅度不会很大，且往往在疫情形势好转之前率先见底。
风险提示：1）计算结果对参数较为敏感；2）政策实施对疫情有较大影响。
如何定量评估流行病的传染性强弱？
SEIR模型
SEIR是经典的流行病传播模型，这里我们只介绍其原理和框架，然后引用一些参考文献的结论计算基本传染数，而不对微分方程的推导过程展开讨论。SEIR模型将整个社群人口分为四类：
易感人群S类：对该病毒没有抗体，接触传染期人群可能被传染；
潜伏期人群E类：被感染后处于潜伏期的人群，还没有传染性，平均时长记为TE；
传染期人群I类：代表潜伏期之后已具有传染性的人群，传染期的平均时长为TI；
隔离态人群R类：治愈并获得免疫或被有效隔离或死亡，既不能传染他人也不能被传染。特别说明：
1）SEIR模型中潜伏期（E类）是特指已感染但还没有传染性的时期，设这一时期的平均时长为TE。但通常对潜伏期的理解是已感染但还没有症状的这段时间（卫健委专家表示新冠病毒的平均潜伏期在10天左右），本次新冠病毒已发生过在没有症状的时候传染给下一个人的案例。国家医疗专家组成员、北京地坛医院感染性疾病诊疗与研究中心首席专家李兴旺表示，“从一般的呼吸道传染病规律来讲，潜伏期末期，病人将要发病的时候可能具有一定的传染率，早期一般不应该有传染性”。因此，本文假定平均为10天的无症状潜伏期中，没有传染性的潜伏期平均为6.5天（参考了SARS的相关数据），最后的3.5天具有传染性，计入传染期TI。2） 传染期的平均时长TI是指病人具备传染性后到不会接触不特定的易感人群为止，比如住院或居家卧床等，而不是从病人发病到痊愈的时间。查询了大量文献，发现在对各类疫情建立SEIR模型时，设置的传染期时长通常都在3到5天，甚至更短，这显然不是病人发病到痊愈的时间，比如SARS很多患者都需要住院治疗很长时间。这样理解下，症状越重的病毒其平均传染期TI就越短，因为患者会很快住院或散失活动能力，而症状越轻的病毒则可能传染期较长。 假设一个传染期的人（I类）与易感人群（S类）接触后易感者被感染进入潜伏期（E类）的概率为β，一个处于潜伏期（E类）的个体单位时间内将以的概率转变为传染期（I类）个体，一个传染期（I类）个体单位时间内将以的概率转变为隔离态（R类）个体。图表1：SEIR模型其中S(t)、E(t)、I(t)和R(t)分别表示t时刻社群中处于易感染态、潜伏期、传染期和隔离态个体数目，N表示社群中个体总数目，且N=S(t)+E(t)+I(t)+R(t)。据参考文献[3],基本传染数可表示为：来源：百度地图，各地卫健委简化微分方程组，建立数值模拟模型 逐一来理解四个微分方程：4)第四个微分方程是单位时间隔离态的人数等于新增的从传染期进入隔离态的人数。离散形式：
由此，我们得到3个离散的方程式，反映的是三个状态量的迭代关系，即如果知道所有参数值，就可以从t-1时刻推导出t时刻的值，确定一个起始点就可以得到整个序列。 数值模拟求解三个迭代方程涉及到3个参数：β、γ1、γ2，根据它们的定义可以转化为另外三个参数：R0、TE、TI，其中:
3）约束条件：a.
b.在1月14日至21日期间潜伏期人群的平均人数在5257左右。两个待定参数，两个约束条件，理论上是可以求解的。4）模拟优化过程：参数R0、TI各自模拟区间内随机取值，通过3个迭代方程可以计算出E(t)、I(t)、R(t)的时间序列（以2019年12月8日为起点出现首个确诊患者），通过时间序列可以计算λ以及1月14日至21日期间潜伏期人数，判断约束条件是否接近。大量模拟后取约束条件最接近的参数值，作为近似解。 参数优化结果：
基本传染数R0比之前大多估计都高一些，传染期在不加干预的情况下长达7天，可能与有大量的轻症患者以及发病前夕无症状时就具有传染性有关。因为没有症状或者症状很轻会使得病人不会停止活动、更不会住院，平均传染期就变长了。
高传染性下，防控措施能达效果吗？自然传染数是指在前期没有额外的防疫干预措施的情况下，平均每个病人能传染的人数。根据传播动力学[1]：
其中k是一个有传染能力的患者平均每天与易感人群的接触次数，b是每次接触传染成功的概率，D是可以传播的时间，可以认为等同于TI。 新冠病毒R0明显高于SARS主要是因为平均每个病人传播的时间更长 比较新型冠状病毒和SARS，考虑到新冠病毒传播初期刚好遇到春运，人员流动性高，假设单位时间接触次数比SARS时期高10%。图表5：新型冠状病毒与SARS对比可见，新冠病毒明显高于SARS主要是因为平均每个病人的传染期时间更长，而单次接触的传染概率可能略低于SARS。新冠病毒传染期格外长可能主要是因为其有大量的轻症患者，不会及时就医或采取措施，而会较长时间处于活跃的传染状态。国家卫健委的专家：“说它是新型冠状病毒肺炎，不如说是新型冠状病毒感染，因为检测阳性的患者不见得一定出现肺炎。就像流感、重感冒，有的轻症患者不需要住院便能治愈。病例数虽多，但大多数都是轻症。”而恰恰是轻症患者较多增加了不小的防控难度。对防控措施应该有信心采取防控措施就是要降低k、b、D这几个参数值，从而降低传染数。虽然新增确诊数还在不断攀升，但从1月23日之后采取的一系列措施来看，我们应该对控制疫情有信心。在流行病传播模型里，只要将基本传染数R0降到1以下，疫情就不会失控。 假定R0就如前文测算的5.38，远高于SARS。首先，交通管制、取消聚会等措施使得人员流动大幅下降，单位时间的接触次数k下降70%甚至更多都有可能。其次，通过宣传和提醒使得民众对疫情具有极高的警惕性，有轻微症状的患者大多也会尽快就医，从而减少传染期的长度，假设传染期的长度由7天降为4天，则D下降43%。第三，通过佩戴口罩、消毒等防护措施，单次接触传染的概率假如能下降30%。则传染数就下降为5.38*（1-70%）（1-43%）（1-30%）=0.64，明显小于1，疫情就能较快控制住了。事实上，WHO专家曾表示，中国在2003年非典期间后期，严格的防疫措施使得SARS的基本传染数下降到0.4。防控措施的效果体尚未体现到新增确诊数 为什么采取了这么严格的防控措施，新增确诊数还在不断攀升呢？首先，可能由于确诊检测试剂盒上市较晚，一些地方的基层医院缺乏确诊的条件或者确诊效率很低，造成确诊数据失真。可能存在前些天发病的病例过了几天才有条件确诊，这就可能造成近日新增确诊数大幅增加。其次，1月24日采取了进出武汉的交通管制，之前大量的人口从武汉外流到各地，按照10天左右的潜伏期，现在是发病的高峰期，随着他们的确诊，跟他们密切接触者大多会进入隔离观察，隔离观察者中还是会有比较多的潜伏期感染者，因此接下来新增确诊数可能还会保持高位，直到这批隔离观察者的潜伏期过去，被确诊，之后的新增确诊数就会大幅下降。从时间上看，可能要到2月10日之后。 以上海为例，其近日公布的新增确诊病例中，大部分还是曾湖北居住或旅游，说明目前新增的确诊病例主要还是原来潜伏期的感染者，严格的管控措施产生的效果还没有反映到新增确诊数据上来。图表6：上海每日新增确诊病例中有湖北居住或旅行史的占比基于数学模型的全国新冠病毒患病人数预测：
我们能打赢这场防疫战
利用前文建立的迭代方程组和估算的参数，可以对未来疫情发展进行预测。但在采取严格防控措施后，基本传染数下降到什么水平目前还没有数据可以估算，因为新冠病毒的潜伏期比较长，当前的新增确诊数据还更多地受之前潜伏期感染人数的影响。 疫情发展预测：中性与乐观假设中性假设：
1）假设自1月23日后，民众减少出行和聚会，单位时间接触次数下降50%，由于提高了对疫情的警惕，有轻微症状的患者会及时就医，使得传染期的时长从7天降为4.5天，由于佩戴口罩、消毒等，单次接触传染的概率下降15%。则基本传染数将从5.38降为1.47。2）假设自2月15日起，由于隔离观察制度的严格实施，大部分新增病例都在隔离观察名单里，基本失去了再次传染的可能，基本传染数下降到0.5，稍高于非典后期（WHO研究认为非典后期在采取严控措施后R0降到了0.4）。乐观假设：1）假设自1月23日后单位时间平均接触次数下降55%，其余与中性假设1）不变。2）假设自2月10日起，隔离观察制度开始起作用，基本传染数下降到0.4。在这两种假设下，可以对未来累计病例数、新增病例数等作出预测，结果如下：1）当前（2月1日）估算的累计病例人数为22732人，与公布的累计确诊病例相差较大，可能主要是因为前期缺乏高效的确诊技术手段，叠加上流感季，武汉医疗资源严重不够，加上新冠病毒的轻症患者较多，症状跟感冒类似，可能未确诊前就已经痊愈了。2）预计未来新冠病毒总患病人数最高将达到88500人（中性假设下）或58000（乐观假设下）人左右，新增确诊人数将在2月17日（中性假设下）或2月11日（乐观假设下）左右出现拐点。3）当前新增确诊病例数大幅攀升（新增病例超2500）可能与检测确诊效率提高有关，预计将很快回落（模型中性假设下预测值是1700左右）。1月25日之前，公布新增确诊数大幅低于模型预测值，26日之后则大幅高于模型预测值，可能是因为确诊检测试剂盒供应充分了，将之前已经发病而没有确诊的病例确诊了。如果看1月24到2月1日期间总的新增确诊量，则模型预估和官方公布的数据很接近，分别为13469人和13582人。图表7：中性假设下模型预测的新冠病毒确诊数和实际公布的新增确诊数说明：2009年H1N1美国的数据是CDC事后建模预估的，当时由于疫情发展失控，没有统计和确认病例数了，据事后的研究，当时H1N1在墨西哥境内死亡率为2%，在其他地区的死亡率为0.1%左右。
股市短期可能反应过度，若超跌是买入机会 2009年美国爆发H1N1流感疫情，最后完全失控。期间美国股市仅在2009年5月、6月有小幅回调，幅度-6%左右。美国股市对疫情冲击反应较小。 2003年SARS疫情期间，港股和A股做出了明显反应，下跌幅度-14%左右，但随后都反弹明显，且反弹时间领先于疫情拐点。可见，股市可能在情绪影响下短期反应过度，随后会修复。 在1月23日发布的《关注新型肺炎“概念股”的估值偏差》文中提到，“事件性的、一次性的利好，即使对当期业绩影响很大（盈利翻倍），对股价的理论影响也很小”，其实反过来也一样，事件性、一次性的利空影响，对公司价值的影响也是较小的，股价不应该反应过度。当然利空与利好不同的是：短期利空（如收入剧减、资金回笼困难等）可能引起公司短期周转困难从而引发更大风险。因此也建议政府相关部门出台一些措施，对疫情中遭遇短期困难的企业给予支持和帮助。图表10：2009年H1N1流感爆发期间美股表现数据来源：Wind，WHO图表11：SARS期间美股、港股、A股表现郑重声明：本文版权归原作者所有，转载文章仅为传播更多信息之目的，如作者信息标记有误，请第一时间联系我们修改或删除，多谢。
