{SiteName}
医学概念
医学核心
医学条件
医学类别
医学创始
国内机构

硬核循证R0951传10R0,不能

                            

当你看到有人使用R0来估算奥密克戎BA.2大流行期间上海的感染人数时,你要警惕了,因为你很有可能正在被R0的错误使用所误导。

为什么这么说?因为

社交媒体上所用的R0其实指的多是Re/Rt。

R0并不是病毒本身的属性,它会随时间和空间波动,很可能使用的R0值已经过时了或者并不适用于上海的地方语境;

R0不是通过测量得到的,而是通过计算得到的;

R0带给你的是一种平均值假象,它是统计学指标,代表着某种数学期望值;

R0只是一个生硬的监控指标……

什么是R0?

R0(Rnaught,基本传染指数)指的是,在完全易感人群中,由典型病例引起的二次感染人数的平均数量。世界卫生组织在《冠状病毒病(COVID-19)基本用品预测工具》[1]给出了R0的计算方式,即:

R0=传染期*每人每天的接触者数*每个接触者的感染概率

这三个变量可以更具体地被描述为:传染性持续时间、个体之间接触的类型和频率、以及传染性个体和易感个体之间的传播概率[2]。

R0源于人口统计学,直到年才被用于流行病学中[3]。今天,当我们提到传染病的感染数时,R0是一个非常重要的衡量参数。以至于在新冠疫情期间,我们经常能看R0出现在社交媒体的公开讨论中。实际上,不止在新冠疫情期间,历次传染病爆发时,都有研究者对其R0进行计算,维基百科整理了历次广为人知的传染病的R0值如下(仅作参考):

R0,能告诉我们什么?

仅仅通过R0这个数值,我们可以获得关于传染病的什么信息?

首先,R0可以帮助我们衡量病毒的传染性。我们经常能够看到类似的这么一句话:

“如果R0<1,则意味着疾病得到控制或传播速度不太快;

如果R0=1,则平均1人可以传播到另1人;

如果R0>1,那么这种疾病可以从1个人传播到更广泛的人群中”。

此外,R0还可以帮助我们预测逆转传染病流行的免疫者比例P。一般来说群体免疫值,即阻断人群传播所需的免疫个体的临界比例,可以作为免疫接种计划的目标,以阻止传染病的传播。这也是强调疫苗接种率的原因之一。

P=1-1/R0

不仅如此,R0还可以帮助预测流行病的规模、平均感染年龄等[4]。

因此R0是一个有价值的流行病学概念。自新冠疫情突发以来,在评估无论是早期在武汉流行的毒株,还是之后的德尔塔,抑或奥密克戎及其变异毒株的传染性时,研究人员都分别计算出了它们不同的R0值,这些不同的R0值也成为国家或组织的应对疫情不可或缺的参考。

①中国疾病预防控制中心在《新型冠状病毒疫情进展和风险评估》(年1月28日)中指出新冠病毒原型株的R0在2-3之间[5];周涛,刘权辉等人以《人民日报》和丁香园发布的新型冠状病毒感染肺炎疫情实时动态数据为基准,估计新冠病毒原型株的基本再生数在2.8-3.3之间[6]。

②德尔塔毒株最早于年10月在印度发现,高传染性使得德尔塔毒株迅速在世界范围内传播。我国学者李文艳,杜志成等人根据中国疾病预防控制信息系统传染病报告信息管理系统和广州市荔湾区CDC发布的感染者(确诊病例和无症状感染者)相关信息得到德尔塔毒株的R0=5.1[7];根据此前科技日报报道,德尔塔+变异株的R0值约为5-8[8]。

③年11月9日,南非首次从病例样本中检测到一种新冠病毒变异株,即奥密克戎。根据此前科技日报,奥密克戎毒株的R0值是7-8[9]。由香港中文大学(中大)医学院(中大医学院)JC公共卫生及基层医疗学院和赛马会兽医及生命科学学院领导的研究团队给出奥密克戎的R0值为8[10]。不过随着奧密克戎本身的持续演变,其亚型株BA.2迅速取代了最初流行的BA.1,也被认为是此次香港、上海疫情主要流行毒株。3月25日,医院感染科主任张文宏在疫情防控发布会上表示,此次在上海流行的奥密克戎BA.2毒株的R0值是9.5[11]。医院院长吴明贤给出的国际上估算的R0则为10[12]。

④年1月19日英国首先检测到了BA.1和BA.2变异株重组后的新变体:XE。一些国家也已经出现了XE的感染病例,但截至目前我国暂未检出。世卫组织指出,XE重组毒株的传播性相比BA.2毒株要高出10%,但由于数据样本太小,还需要进一步观察[13]。

R0,不能告诉我们什么?

是不是新冠变异毒株的R0值越高?新冠疫情就会越来越严重?正在上海流行的奥密克戎BA.2,甚至已经开始在英国出现的XE的R0值越来越高,甚至能够达到一传十?

当我们以R0作为奥密克戎BA.2以及新近变异毒株XE的判断标准时,需要极为警惕。

R0容易解释,但是难以计算。R0值具有“外柔内刚”的特性,所谓“外柔”是因为它的定义及解释使其看起来是一个十分简单、容易理解且容易利用其进行计算的指标,而“内刚”则是其专业性与复杂性的内核,这恰恰被其表征所掩盖。

1、R0是平均值假象

R0代表的是一个感染病例平均产生的二代病例数。而平均值是统计学中十分重要的一个概念,这种随机数的平均值可以被称为数学期望,这种期望值与我们现实生活中存在很大的差距。

以新冠毒株德尔塔为例,据澎湃新闻此前报道,德尔塔+变异株的R0值约为5-8[14]。当R0=5时,那就意味着每个感染者都会将病毒传染给另外5个人。但是,这也可能意味着,在10个人中,倘若有感染者是感染人的“超级传播者”,那么其余9名感染者就不会感染任何人。当然,这是十分理想化且极端的假设。不过在现实生活中,也出现了类似的情况。据央视新闻,通过对陕西(西安)近一周通报的确诊病例详情梳理发现,年1月4日发布的本土确诊第60号病例已感染至少43名密切接触者[15]。而此前根据俄罗斯塔斯社年8月7日报道,俄罗斯曾出现一名传染了人的新冠病毒超级传播者[16]。因而,我们无法确定人群中“超级传播者”的数量,更无法排除“超级传播者”的存在。

2、R0是在完全易感人群前提假设下计算得出

什么是「完全易感人群」?完全易感人群即是指没有任何预防手段介入并且所有人对此病原体没有免疫力的情况下。因此R0值可能准确的测量时间是新出现的传染性病原体爆发期间,因为这时候没有任何阻断措施,病毒是在是完全易感人群中传播[17]。而随着各种阻断措施如封控、接种疫苗的出现,R0值的计算已经与现实世界中的情况大不相同。

因此,当你看到有人在使用R0预测未来感染规模时,它并不是完全准确,因为预测结论成立的前提要么是被计算的人群是完全易感人群,要么是我们能够提前计算出人们对于这种病毒预先存在的免疫力[18]。

3、R0是回顾性的估计值,会随时间波动

R0是回顾性(retrospectively)估计。R0值通常有两种计算方式,一种是根据连续流行病学数据,即使用在流行病开始时获得的个人层面的接触者追踪数据来计算R0;另一种则是使用数学模型(构建常微分方程)进行回顾性估计,即是根据人口层面的数据[19],随着总体患病率的上升,个体接触追踪变得困难或不可能,因此研究者会使用人口水平的累积发病率数据进行计算,这也是更常用的一种方法。个体层面和人口层面的方法可能会产生不同的数字,这也是我们常常看到不同研究人员给出的R0值各不相同的原因之一。

R0值几乎总是从数学模型中估计出来的,它取决于建模过程中做出的许多决定。对于新冠病毒R0值的计算,不少研究者采用了SIR传播模型即易感-感染-恢复模型,或SEIR传播模型即易感-暴露-感染-恢复模型等等,不同传播模型估算出的R0值并不相同。此外,还要考虑人口混合和接触模式等等[20]。Pandit还指出,不论是个人层面的建模还是人口层面建模,R0的估计都取决于感染者每天接触的平均人数、每次暴露被感染的概率、感染持续时间(包括无症状感染期)、人口规模(尤其是易感人群的比例)以及康复或死亡率[21]。而这些变量,其真实值通常是未知的、难以或不可能直接测量得到。这也是我们常常看到的R0值通常是在一个范围区间内的原因之一。

而另一种情况,正如Delamater在论文中指出的那样:经常报告的麻疹R0值为12-18,这是基于-年在美国(R0为12.5)和-年在英格兰和威尔士(R0为13.7-18.0)获得的数据,但是随着时间,人类在社会和地理中的组织方式已然发生的重大变化,进而使得这些历史数值极不可能与当今的流行病学现实相匹配[22]。

所以,毫无疑问你所使用的R0值本身就是使用过往数据计算出来的估计值,并且它很可能已经过时。

4、R0不是毒株本身的属性,是人们创造出来的数值,还会随空间波动

R0是一个无量纲量,所谓无量纲量是相对于有量纲量而言,量纲(dimension)是指物理量的基本属性。这看起来很专业,那么通俗地说,“有量纲量”就是能够依赖测量量纲得到的“常数”,而“无量纲量”实际上就是我们自己“创造”的[23]。因此,R0并不是毒株本身的属性,是人们创造出来的指标。

根据R0的计算公式,可以看到R0值取决于传染性持续时间、个体之间接触的类型和频率以及传染性个体和易感个体之间的传播概率。即使公式中毒株传染性持续时间以及传播概率是生物学常数,当个体之间接触的频率或者类型发生改变,那么R0也会随着时间与空间发生波动[24]。

我们具体来看一个例子。Benjamin等人在其论文中指出,世界不同地区的R0值存在显著的异质性。他的举例子是这样的:在印度的各个邦中,年H1N1病毒的R0值在1.03到1.之间;秘鲁的则是从1.2到2.2不等。甚至地理位置很近地区也有完全不同的R0估值:中国估计R0均值为1.68,而日本最初估计R0均值为2.3,后来降至1.21至1.35,等等。他还举例指出,并不是所有R0的后续估计都是向下偏倚的,比如弗雷泽等人是墨西哥最早估算R0的人之一,他们提出R0为1.4到1.6,但是几个月后,另一个团队却估计R0在2.3到2.9之间。进而他指出,几乎没有证据表明,一个地理区域的R0值适用于另一个地理区域[25]。

所以,很可能你所使用的R0值并不是根据你想要计算地区的人口结构、接触方式等计算得出。

5、更大的R0值并不一定意味着更高的传染数量

流行性感冒的R0值约在2-3,据中国疾控中心每年发布的全年法定传染病疫情概况,年流感的发病数约为.5万;于年在中国广东发生的SARS,其R0值约在2-5,根据世界卫生组织,全球累计报告非典型肺炎临床诊断病例例[26]。由于R0值是一个波动范围,且不同的持续时间、不同的地区以及人群都会不同,因而这样类比似乎并不科学和专业。

但是,当我们仅从数学意义上使用R0这个无量纲量来预估的话,相同的R0值应该会有相同的感染人数,但是根据以上数据我们可以看到在实际社会中,由于各种阻断措施的介入,实际的传染量并不会如同R0所预估的那样发展。

不仅如此,根据上述统计的数据来看,目前流行的奥密克戎BA.2以及最新的变异毒株XE的R0值并不是历次传染病中的最大值,大大低于麻疹的R0值(12-16),因此不必因为R0值的变化而感到过度恐慌。

6、大众/社交媒体误用的R0值,其实是Re值

前面提到,R0值计算的前提假设是完全易感人群,所以即便通过疫苗手段使群体获得免疫,R0值并不会降低。当我们去衡量疫苗效果或者群体中已经有部分成员获得免疫效果时,更合适的度量标准是有效繁殖数Re(有时也被称作Rt),它与R0相似,并经常与R0混淆,但它不假设人群完全易感,因此可以用具有免疫成员的人群来估计。

Re指的是一个群体中在任何特定时间可以被一个人感染的人数。它会受到感染者人数以及易感人群人数的影响,人们的行为(如社交距离)也会影响Re[27]。当Re1时,每个受感染的个体将疾病传播给不止一个人,疾病可能在人群中传播;如果Re1,并不是每个病例都将导致另一个人发生新的感染,传播将停止(然而,可能发生小的传播链)[28]。

前面提到,群体免疫值(P=1-1/R0),可以作为免疫接种计划的目标,以阻止疾病的传播,而在测试疫苗接种的有效性时,我们应该使用可用于具有免疫成员的人群的有效繁殖数Re[29]。

也就是说我们可以粗暴的将R0理解为在完全易感人群中的预估值,而现实社会中的实际值我们需要随时

转载请注明:http://www.lingshanl.com/yxcs/86564.html

  • 上一篇文章:
  • 下一篇文章: 没有了