网络群体与市场 课件.ppt
传染病的网络传播模型,(第21章),疾病传播模型:基本要素,多大范围(规模)?要不要采取措施?许多传染病会自生自灭,例如季节性的多长时间?会不会反复(“死灰复燃”)对患者影响的严重性(生命、损伤、不适),接触网络,疾病特征,传染概率,关心的问题:,流行病传播与社会行为传播的差别,已经讨论过的传播模型总体模型总体(占比)效应网络模型社会关系网络结构 基于门槛值的模型流行病传播特点传播原因:传染传播途径:接触关系网络结构传播方式:具有随机特性,决策模型,随机模型,分支过程:一个最简单的疾病传播模型,分支过程形成树结构接触网络,最初一个人携带病菌进入人群,以一个独立的概率 p 传染给遇到的每个人。假设疾病感染期间每个人遇到 k 个其他人,这个模型虽然简单,但便于引出刻画疾病传播过程的一个关键概念,p反映了疾病的传染性,传染过程停止的条件?,为什么随着层数往下,被传染的人越来越多?,为什么随着层数往下,被传染的人越来越少?,基本再生数(reproductive number),基本再生数:由单一个体引起的新发病例数期望值,记为R0,R0=p*k 如果R0 1,则疾病持续在每一波以一定的概率至少感染一个人,基本再生数的特点,因为R0=p*k,当R0接近于1时,稍微改变接触人数k或传染概率会对结果产生很大影响。R0接近于1时有一个“刀刃”(knife-edge)特性如果R0略低于1,稍微增加传染概率p;结果可能会使R0 高于1,造成一个突然的疾病爆发略微减少疾病的传染性,可以将R0减小到1以下,消除疾病大范围流行的风险结论:当门槛值R0=1附近,社会应该付出努力减小基本再生数(隔离,戴口罩等就是措施)实际的疾病传播要复杂得多,分支结构没有考虑三角关系,弱连接,SIR疾病传播模型,SIR模型可以描述更一般网络结构的传播情况,接触网络中的节点在流行病传播过程中经历三个状态(susceptible,infectious,removed):易感状态(S):易感者,有可能被一个患病的网络邻居节点感染传染状态(I):患病者,以一定的概率将疾病传染给处于易感状态的邻居移出状态(R):痊愈并免疫者,当一个节点经历了完整的病期康复后,就不再会被传染,也不会再对其他节点造成威胁,相当与从接触网络中移出了,SIR模型描述,SIR模型流行病的传播过程取决于接触网络结构,以及两个量值:p(传染的概率)和 tI(传染期的长度)最初,一些节点处于状态 I,所有其他节点处在状态 S。每个进入状态 I 的节点v在固定的步骤 tI 期间具有传染性。在 tI 的每一步,v 以概率 p 将疾病传染给它的处于易感状态的邻居经过 tI步后,节点 v 进入移出状态(R),成为网络中的惰性节点,不会再被传染,也不会向其他人传播,易感 S,传染 I,移出 R,接触网是一个有向图图中有三种不同的节点白色(易感)粗边界的灰色(传染)细边界的灰色(移出),一个 tI=1的SIR模型示例,接触不一定传染,非树型网络的R0不再有“刀刃”特性,网络由每层两个节点组成向右无限伸延,tI=1,传染概率p为2/3,最左边的两个节点是最初传染者R0是由一个节点造成的新病例数,预期为4/3 1每条边不传播疾病的概率为1/3,所有四条边都没导致疾病传播的概率为(1/3)4,于是每一层成为最后一层(疾病传染终止层)的概率至少是1/81,这种结构相当于形成了一个路障,疾病无法超越它向远方传播,SIR模型(运行)的静态图示,接触网络每条边预先投掷硬币(按照给定传染概率p),投掷硬币成功的边定义为开放边;其余的边定义为 阻塞边利用开放边和阻塞边描述流行病传播过程:最终受感染的节点正是那些能够从最初感染节点沿着在网络中随机选定的 开放边能够达到的节点注意:每次运行结果可能不同,图中,粗黑线为开放边。注意到不是所有开放边都起作用的(例如 wr),SIR模型的形式化描述,易感人数、传染人数、移出人数分别为S(t),I(t),R(t),占比分别为s(t),i(t),r(i):个人单位时间与其他人平均有效接触次数:单位时间治愈人数 SIR模型描述为:,来自最早提出SIR模型的文献:Kermack W O,McKendrick A G.Contributions to the mathematical theory of epidemics,part IJ.Proceedings of the Royal Society of London A,1927,115:700721,一个给定和的SIR时间曲线,设=1,=0.4,s(at start)=0.99,i(at start)=0.01,r(at start)=0疾病什么时候消失呢?取决于:/。如果发病率高于痊愈率,则疾病会持续,否则疾病将消失,利用SIR模型分析SARS传播,利用SIR模型中di/dt方程模拟北京在2003 年4 月18 日到6 月16 日的每天新增病例模拟中使用的参数是:=0.148,=0.85,初始值为:s(0)=1199,i(0)=1(这里总人数应该看作是接触到病源的易感、传染、免疫的人数之和,而不是北京市总人口),SIS流行病传播模型,SIR模型只适用于一次患病终身免疫的疾病。SIS模型描述了更一般的疾病传播。节点结束传染状态后再回到易感状态,节点在 S 和 I 两种状态中交替,因此得名SIS模型最初,一些节点处于状态 I,其余节点处于状态 S每个进入状态 I 的节点 v 在 tI 期间内具有传染性,在 tI的每一步,节点以概率 p 将疾病传染给其处在易感状态的邻居经过 tI 步后,I 节点不再具有传染性,返回状态 S,易感 S,传染 I,SIS模型的特点,SIS模型下的传染病的流行周期会比SIR模型的要长基本再生数具有刀刃特点:可以证明,对于某些类型的接触网络,存在特定传染概率p,在其附近,网络中的SIS传染病在一边表现为“快速消失”,在另一边则会“持续一个很长的时间”。该传染概率临界值以微妙的方式依赖于网络结构。,SIR与SIS的关系,SIS模型可以看成是SIR模型的特例设 tI=1,将患病节点v在每一个时间步视为“不同的个体”,可以将SIS模型表示为SIR模型在每个时间周期 t=0、1、2、3,为每个节点分别创建一个节点副本,对于原网络中v连接到w的边,复制 t 时刻扩展网络中 v 到 t+1时刻w的边,同时被w传染,SIR与SIS的关系,SIRS 流行病传播模型,结合SIR和SIS,疾病有暂时免疫特性:受传染节点治愈后会暂短地经过R状态,然后再次进入S状态SIRS模型开始某些节点在 I 状态,其他在S状态进入 I 状态的节点v在固定时间 tI 内具有传染性,在 tI 的每一步以概率p将疾病传播给易感状态的邻居I 状态节点经过 tI 步后,进入状态R,经过 tR 后再次回到S状态,传染I,易感S,移出R,疾病传播过程受量值 p、tI 以及 tR 的影响,同时受到网络结构的影响,tI,tR,SIRS模型下可能看见传染病流行的震荡效应,T,S,U,W,X,V,p=0.5,形象过程,Watts-Strogatz小世界模型上运行SIRS模型,I=4,R=9Ninf(0)=0.1模拟c=0.01,0.2,0.9三种情况,结论:当弱连接(远程)边占比较高时,所有节点被传染的周期趋于同步。经历被传染,持续相同的时间,进入隔离,之后又同时被传染,受感染人数占比,模拟时间,c:弱连接边占比,为什么会有同步震荡?,当c较小时:00.5,网络呈现较高的有序的聚集特点,疾病周期会在一个局部区域内完成,之后进入免疫期(不活跃期)当c值较大时:0.5,弱连接打破了有序的聚集区域,刚刚痊愈的个体可能又通过一个弱连接被再次感染。因此整个网络成同步状态,模型中可以考虑的一些其他因素,传染期不同时段传染性可以不同不同个体传染性的不同患病者可能在传染期间不同的阶段恢复接触网络的边具有时间特性其他非网络因素(例如蚊子),一个生物学基因遗传的研究,线粒体夏娃理论:1987年,Rebecca Cann、Mark Stoneking和Allan Wilson 在Nature 杂志上发表一篇论文引起轰动线粒体DNA严格按照母系遗传现代全人类的线粒体DNA基本相同于是,每个人都有一条母系祖先的线索,称为母亲的血统。所有这些血统源自10万至 20万年前(14万年)的一个单一女子,很可能是在非洲,她是我们当代所有人的母系祖先。,一点生物遗传知识,线粒体遗传的网络结构动力学解释,Wright-Fisher单亲祖先模型:每一代的人口数量保持在固定的N当前一代N个个体产生新一代的N个个体每个后代由一个单亲产生,从当前一代中独立地均匀随机选择,这样,当前一代的一个个体可能会有多个孩子 这样,就可以想到一个网络,节点是个体,代际分层,边代表单亲关系,每个个体都向上连接到自己的母亲;时间自上向下展开,最底层表示当前有N个个体从底部开始,任何一个个体都可以沿着向上的边追溯他她的单亲血统能否说明这些边最终将汇聚到一点?,简单的概率分析,可以发现当前给定数量的个体在向上追溯单亲祖先的过程中会不断发生汇聚(发现同一个母亲),于是独立的线索越来越少,最终汇聚到一点这个距当代最近的汇聚点就是线粒体夏娃,基因遗传学的解释,当时肯定同时生活着许多女人(以及男人),只不过她们的线粒体基因没有遗传到现在。不仅是线粒体基因存在这种情形,人的任何一种基因,只要把它们的突变过程追溯得足够远,总能找到一个共同祖先,这些祖先生活的时间可能各不相同,甚至未必是人,有的可能要追溯到人进化出来之前的某个生物,Wright-Fisher单亲祖先模型的意义,可以分析单亲物种的遗传特点可以分析双亲物种中那些只由单亲遗传的现象,如线粒体DNA分析社会继承行为(包括学术继承等),我,父母,祖父母,曾祖父母,曾曾祖父母,上溯到第十代,有多少不同的人是与你有血缘关系的祖辈?,小结,疾病传染的简单网络模型:分支结构、SIR、SIS、SIRS不同模型下的不同现象与问题生物遗传过程也可以自然看成是一个网络,从中可能发现出人意料(或者令人困惑)的现象,作业:练习21.2,设想你了解一群人的接触关系,但不知道其中接触发生的时间段。现观察到某种疾病在某些接触关系中总不会传染,而在另一些接触关系中会传染,称前者为阴性关系,后者为阳性关系。很自然要问是否能求得接触网络中每条边存在的时间段,对应上述阳性关系和阴性关系的特征。我们用简单的五人接触网络描述这样的问题,如图21.23所示。你能否求得每条边存在的时间段,使得疾病有可能从每个节点流向所有其他节点,只有一个例外,就是不可能从节点a传到节点e?如果你认为这是可能的话,描述这样一组时间段;如果你认为不可能,解释为什么没有这样的时间段存在。你能否求得每条边存在的时间段,使得疾病可能从节点a传到节点d,以及从节点b传到节点e,但不能从节点a传到c?同样,如果你认为它是可能的,描述这样一组时间段,如果你认为这是不可能的,解释为什么没有这样时间段存在。,