cstring只获取到第一个数_传染病中有效再生数的计算

最近在忙毕业论文的事情，里边需要评估新冠肺炎疫情里各国疫情防控的情况，通过查阅文献最终选择有效再生数作为评估的依据，但是国内关于计算有效再生数的文献并不多，因此头秃了很久，好在最终还算有了个好的结果，现在把这个问题总结一下。

一、什么是有效再生数

有效再生数Rt（effective reproduction number）代表某一时刻t确诊的某一位感染者，在其感染期内平均会传染多少二代确诊病例。

一般地，可以用Rt 来衡量疫情的实时可传播性，当Rt>1时，流行病疫情将会在人群中迅速传播，而随着Rt值逐渐降低，疫情的发展速度也会逐渐下降，最终当Rt<1后，感染人数会逐渐降低，疫情逐渐解除。

二、有效再生数的计算方法

对Rt的计算与R0类似，通常可以有三种方法：一是直接从病例调查报告入手，通过查看传染病的传染链路弄清谁将疾病传给了谁，而后有效再生数的计算只需要简单地统计每一个染病个体感染的下一代个体数量，最后取一个平均值即可，这种方法是一线公共卫生人员常用的方法；二是用通过指数函数对疫情爆发初期的病例数据进行拟合，得到确诊病例增长率r，尔后有效再生数等于代间隔的概率分布的矩量生成函数在-r取值的倒数；三是在了解传染病传播的特点后，写出其传播的动力学模型，通过拟合动力学模型及每日新增病例数据拟合有效再生数。

But，说实话，这三种方法我都没搞明白怎么操作，第一种不说了，主要难点在于一般人难以获取有效可靠的数据；第二种除了增长率r的计算以外，其他我一概没看懂；第三种的话，难点也在于数据的获取，以新冠肺炎疫情来看，比较适用于SEIR模型，但以SEIR模型来拟合的话，有效再生数的计算如下：