草庐IT

图数据挖掘:基于概率的流行病模型

Orion's Blog 2023-03-28 原文

1 导引

在上一篇博客《图数据挖掘:网络中的级联行为》中介绍了用基于决策的模型来对级联行为进行建模,该模型是基于效用(Utility)的且是是确定性的,主要关注于单个节点如何根据其邻居的情况来做决策,需要大量和数据相关的先验信息。这篇博客就让我们来介绍基于概率的传播模型,这种模型基于对数据的观测来构建,不过不能对因果性进行建模。

2 基于随机树的流行病模型

接下来我们介绍一种基于随机树的传染病模型,它是分支过程(branching processes)的一种变种。在这种模型中,一个病人可能接触\(d\)个其他人,对他们中的每一个都有概率\(q>0\)将其传染,如下图所示:

接下来我们来看当\(d\)\(q\)取何值时,流行病最终会消失(die out),也即满足

\[\lim _{h \rightarrow \infty} p_h=0 \]

这里\(p_h\)为在深度\(h\)处存在感染节点的概率(是关于\(q\)\(d\)的函数)。如果流行病会永远流行下去,则上述极限应该\(>0\)

\(p_h\)满足递归式:

\[p_h=1-\left(1-q \cdot p_{h-1}\right)^d \]

这里\(\left(1-q \cdot p_{h-1}\right)^d\)表示在距离根节点\(h\)深度处没有感染节点的概率。

接下来我们通过对函数

\[f(x)=1-(1-q \cdot x)^d \]

进行迭代来得到\(\lim _{h \rightarrow \infty} p_h\)。我们从根节点\(x=1\)(因为\(p_1=1\))开始,依次迭代得到\(x_1=f(1), x_2=f(x_1),x_3=f(x_2)\)。事实上,该迭代最终会收敛到不动点\(f(x)=x\),如下图所示:

这里\(x\)是在深度\(h-1\)处存在感染节点的概率,\(f(x)\)是在深度为\(h\)处存在感染节点的概率,\(q\)为感染概率,\(d\)为节点的度。
如果我们想要传染病最终消失,那么迭代\(f(x)\)的结果必须要趋向于\(0\),也即不动点需要为0。而这也就意味着\(f(x)\)必须要在\(y=x\)下方,如下所示:

如何控制\(f(x)\)必须要在\(y=x\)下方呢?我们先来分析下\(f(x)\)的图像形状,我们有以下结论:

\(f(x)\)是单调的:对\(0 \leq x, q \leq 1, d>1\)\(f'(x)=q \cdot d(1-q x)^{d-1}>0\),故\(f(x)\)是单调的。

\(f'(x)\)是非增的\(f'(x)=q \cdot d(1-q x)^{d-1}\)会着\(x\)减小而减小。

\(f(x)\)低于\(y=x\),则需要满足

\[f'(0)=q\cdot d<1 \]

综上所述,我们有结论:

\[\lim _{h \rightarrow \infty} p_h=0 \text { when } q \cdot d<\mathbf{1} \]

这里\(R_0=q\cdot d\)表示每个被感染的个体在期望意义上所产生的新的病体数,我们将其称为基本再生数(reproductive number),它决定了传染病病是否会流行:

  • \(R_0\geq 1\): 流行病永远不会消失且感染人数会以指数速度上升。
  • \(R_0\leq 1\): 流行病会以指数速度快速消失。

3 SIR与SIS流行病模型

3.1 模型范式

在病毒的传播中,有两个最基本的参数:

  • 出生率\(\beta\) 被已感染邻居攻击的概率
  • 死亡率\(\delta\) 已感染节点治愈的概率

网络中的节点可以在以下四个状态(S+E+I+R)之间做转移:

  • 易感期(susceptible): 节点患病之前,处于容易被邻居传染的时期,也称敏感期。
  • 潜伏期(exposed):节点已被感染,但是还没具备能力去传染别人。
  • 传染期(infectious):节点已被感染,且能够以一定的概率把疾病传染给那些处于易感期的邻居,也称感染期。
  • 移除期(removed):当一个节点经历了完整的传染期,就不再被考虑了,因为它不会再受感染,也不会对其它节点构成威胁,也称隔离期。

状态转移图如下图所示(图中的\(Z\)表示人工免疫):

其中状态转移的概率由我们前面提到的模型参数\(\beta\)\(\delta\)控制。

3.2 SIR模型

在SIR模型中,节点经历S-I-R三个阶段:

事实上,该模型可用于对水痘和鼠疫的建模,也即一旦我治愈了,那我就永远不会再被感染了。

假设模型满足完美混合(即网络是完全图),则模型的动力方程为:

\[\begin{aligned} &\frac{d S}{d t}=-\beta S I \\ &\frac{d R}{d t}=\delta I \\ &\frac{d I}{d t}=\beta S I-\delta I \end{aligned} \]

处于\(S\)\(I\)\(R\)状态的节点数量随着时间变化曲线如下图所示:

3.3 SIS模型

SIS模型中节点只有S-I两个阶段,它假设已经治愈的节点会立即变为易感节点。节点的状态转移图如下:

这里我们把\(s=\frac{\beta}{\delta}\)定义为病毒的“力量”(strength)。

该模型可用于对流感的建模,也即已被感染的节点经过治愈后会重新回到易感状态。

同样我们假设模型满足完美混合(即网络是完全图),则模型的动力方程为:

\[\begin{aligned} &\frac{d S}{d t}=-\beta S I+\delta I \\ &\frac{d I}{d t}=\beta S I-\delta I \end{aligned} \]

处于\(S\)\(I\)状态的节点数量随着时间变化曲线如下图所示:

3.4 传染阈值

接下来我们考虑SIS模型中的传染阈值(epidemic threshold)\(\tau\)。对于\(SIS\)模型而言,流行阈值可以是任意的。

我们设图\(G\)的传染阈值为\(\tau\)。如果病毒的“力量”\(s=\frac{\beta}{\delta} < \tau\)(这里\(\beta\)指病毒的死亡率,\(\delta\)指病毒的出生率),则疾病的流行就不会发生(它最终会消失)。事实上,图\(G\)的传染阈值\(\tau\)可以表示为

\[\tau=\frac{1}{\lambda_{1, A}} \]

这里\(\lambda_{1, A}\)为图\(G\)的邻接矩阵最大的特征值。这个定理看起来非常神奇,因为我们只用\(\lambda_{1,A}\)就捕捉到了整个图的属性!

以下是在AS图上,当\(s\)大于、小于或等于传染阈值\(\tau\)时的感染节点数量随时间变化图:

如果我们再考虑不同的初始感染人数,则会得到以下的感染人数变化图像:

3.5 一个埃博拉的例子

在一个埃博拉的例子[1]中,设置如下的转换状态:

当设置\(R_0=1.5\text{-}2.0\)时,总死亡人数随时间变化如下:

参考

[1] Gomes M F C, y Piontti A P, Rossi L, et al. Assessing the international spreading risk associated with the 2014 West African Ebola outbreak[J]. PLoS currents, 2014, 6.
[2] http://web.stanford.edu/class/cs224w/
[3] Easley D, Kleinberg J. Networks, crowds, and markets: Reasoning about a highly connected world[M]. Cambridge university press, 2010.
[4] Barabási A L. Network science[J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2013, 371(1987): 20120375.

有关图数据挖掘:基于概率的流行病模型的更多相关文章

  1. ruby-on-rails - Rails - 子类化模型的设计模式是什么? - 2

    我有一个模型:classItem项目有一个属性“商店”基于存储的值,我希望Item对象对特定方法具有不同的行为。Rails中是否有针对此的通用设计模式?如果方法中没有大的if-else语句,这是如何干净利落地完成的? 最佳答案 通常通过Single-TableInheritance. 关于ruby-on-rails-Rails-子类化模型的设计模式是什么?,我们在StackOverflow上找到一个类似的问题: https://stackoverflow.co

  2. ruby - 解析 RDFa、微数据等的最佳方式是什么,使用统一的模式/词汇(例如 schema.org)存储和显示信息 - 2

    我主要使用Ruby来执行此操作,但到目前为止我的攻击计划如下:使用gemsrdf、rdf-rdfa和rdf-microdata或mida来解析给定任何URI的数据。我认为最好映射到像schema.org这样的统一模式,例如使用这个yaml文件,它试图描述数据词汇表和opengraph到schema.org之间的转换:#SchemaXtoschema.orgconversion#data-vocabularyDV:name:namestreet-address:streetAddressregion:addressRegionlocality:addressLocalityphoto:i

  3. ruby-on-rails - Rails - 一个 View 中的多个模型 - 2

    我需要从一个View访问多个模型。以前,我的links_controller仅用于提供以不同方式排序的链接资源。现在我想包括一个部分(我假设)显示按分数排序的顶级用户(@users=User.all.sort_by(&:score))我知道我可以将此代码插入每个链接操作并从View访问它,但这似乎不是“ruby方式”,我将需要在不久的将来访问更多模型。这可能会变得很脏,是否有针对这种情况的任何技术?注意事项:我认为我的应用程序正朝着单一格式和动态页面内容的方向发展,本质上是一个典型的网络应用程序。我知道before_filter但考虑到我希望应用程序进入的方向,这似乎很麻烦。最终从任何

  4. ruby-on-rails - 在混合/模块中覆盖模型的属性访问器 - 2

    我有一个包含模块的模型。我想在模块中覆盖模型的访问器方法。例如:classBlah这显然行不通。有什么想法可以实现吗? 最佳答案 您的代码看起来是正确的。我们正在毫无困难地使用这个确切的模式。如果我没记错的话,Rails使用#method_missing作为属性setter,因此您的模块将优先,阻止ActiveRecord的setter。如果您正在使用ActiveSupport::Concern(参见thisblogpost),那么您的实例方法需要进入一个特殊的模块:classBlah

  5. ruby-on-rails - 如何验证非模型(甚至非对象)字段 - 2

    我有一个表单,其中有很多字段取自数组(而不是模型或对象)。我如何验证这些字段的存在?solve_problem_pathdo|f|%>... 最佳答案 创建一个简单的类来包装请求参数并使用ActiveModel::Validations。#definedsomewhere,atthesimplest:require'ostruct'classSolvetrue#youcouldevencheckthesolutionwithavalidatorvalidatedoerrors.add(:base,"WRONG!!!")unlesss

  6. ruby-on-rails - form_for 中不在模型中的自定义字段 - 2

    我想向我的Controller传递一个参数,它是一个简单的复选框,但我不知道如何在模型的form_for中引入它,这是我的观点:{:id=>'go_finance'}do|f|%>Transferirde:para:Entrada:"input",:placeholder=>"Quantofoiganho?"%>Saída:"output",:placeholder=>"Quantofoigasto?"%>Nota:我想做一个额外的复选框,但我该怎么做,模型中没有一个对象,而是一个要检查的对象,以便在Controller中创建一个ifelse,如果没有检查,请帮助我,非常感谢,谢谢

  7. ruby-on-rails - 如何将验证与模型分开 - 2

    我有一些非常大的模型,我必须将它们迁移到最新版本的Rails。这些模型有相当多的验证(User有大约50个验证)。是否可以将所有这些验证移动到另一个文件中?说app/models/validations/user_validations.rb。如果可以,有人可以提供示例吗? 最佳答案 您可以为此使用关注点:#app/models/validations/user_validations.rbrequire'active_support/concern'moduleUserValidationsextendActiveSupport:

  8. ruby-on-rails - Rails 模型——非持久类成员或属性? - 2

    对于Rails模型,是否可以/建议让一个类的成员不持久保存到数据库中?我想将用户最后选择的类型存储在session变量中。由于我无法从我的模型中设置session变量,我想将值存储在一个“虚拟”类成员中,该成员只是将值传递回Controller。你能有这样的类(class)成员吗? 最佳答案 将非持久属性添加到Rails模型就像任何其他Ruby类一样:classUser扩展解释:在Ruby中,所有实例变量都是私有(private)的,不需要在赋值前定义。attr_accessor创建一个setter和getter方法:classUs

  9. ruby-on-rails - Rails - 从另一个模型中创建一个模型的实例 - 2

    我有一个正在构建的应用程序,我需要一个模型来创建另一个模型的实例。我希望每辆车都有4个轮胎。汽车模型classCar轮胎模型classTire但是,在make_tires内部有一个错误,如果我为Tire尝试它,则没有用于创建或新建的activerecord方法。当我检查轮胎时,它没有这些方法。我该如何补救?错误是这样的:未定义的方法'create'forActiveRecord::AttributeMethods::Serialization::Tire::Module我测试了两个环境:测试和开发,它们都因相同的错误而失败。 最佳答案

  10. ruby-on-rails - Ruby 中的内存模型 - 2

    ruby如何管理内存。例如:如果我们在执行过程中采用C程序,则以下是内存模型。类似于这个ruby如何处理内存。C:__________________|||stack|||------------------||||------------------|||||Heap|||||__________________|||data|__________________|text|__________________Ruby:? 最佳答案 Ruby中没有“内存”这样的东西。Class#allocate分配一个对象并返回该对象。这就是程序

随机推荐