复杂网络的免疫策略

1、复杂网络的免疫策略第1页，共41页，2022年，5月20日，13点27分，星期二大纲基本的复杂网络免疫策略改变假设条件：局域搜索免疫改变免疫对象：删除边的免疫改变免疫原则：多重图形剖分免疫对于有向网络免疫的思考第2页，共41页，2022年，5月20日，13点27分，星期二基本的免疫策略 目标：通过对部分人接种而有效地控制疾病的传播 基于局域信息 免疫uniform immunization(均匀免疫) acquaintance immunization（熟人免疫） 基于全局信息 targeted immunization（目标免疫） 第3页，共41页，2022年，5月20日，13点27分，星期

2、二均匀免疫均匀免疫，顾名思义完全随机的从网络中选择一部分节点进行免疫。它对于度数大的节点和度数小的节点平等对待 在无标度网络中对应的免疫临界值 第4页，共41页，2022年，5月20日，13点27分，星期二均匀免疫第5页，共41页，2022年，5月20日，13点27分，星期二熟人免疫随机选择比例为p的节点，然后再从这些选择的节点中随机选择一个邻居节点进行免疫 由于度数大的节点也就意味着有更多的节点与之相连，所以熟人免疫比均匀免疫的效率要好得多第6页，共41页，2022年，5月20日，13点27分，星期二熟人免疫第7页，共41页，2022年，5月20日，13点27分，星期二目标免疫根据无标度网络

3、的不均匀特性，可以进行有选择的目标免疫，即选取度数大的节点进行免疫 在BA无标度网络中，目标免疫对应的免疫临界值为 第8页，共41页，2022年，5月20日，13点27分，星期二目标免疫第9页，共41页，2022年，5月20日，13点27分，星期二不同免疫策略的比较在网络规模为106，幂率指数 在2-3.5之间变化的无标度网络中不同策略对应的免疫临界值均匀免疫(空心圆)熟人免疫（空心三角形）目标免疫（空心正方形）图1（参考文献3）第10页，共41页，2022年，5月20日，13点27分，星期二局域搜索免疫熟人免疫假设条件为已知当前节点的度目标免疫假设条件为已知所有节点的度假设已知邻居节点的度信

4、息，怎样进行免疫呢？第11页，共41页，2022年，5月20日，13点27分，星期二1967年，哈佛大学的社会心理学家Stanley Milgram就设计了一个连锁信件实验4。他将一套连锁信件随机发送给居住在内布拉斯加州 奥马哈的160个人，信中放了一个波士顿股票经纪人的名字，信中要求每个收信人将这套信寄给自己认为是比较接近那个股票经纪人的朋友。朋友收信后照此办理。最终大部分信在经过五、六个步骤后都抵达了该股票经纪人。 第12页，共41页，2022年，5月20日，13点27分，星期二Six degrees of separation成功传递信件的前提是 已知朋友中成功传递信件的程度 类似于该实

5、验过程，提出了局域搜索免疫(local search immunization strategy) 第13页，共41页，2022年，5月20日，13点27分，星期二第14页，共41页，2022年，5月20日，13点27分，星期二局域搜索免疫第15页，共41页，2022年，5月20日，13点27分，星期二在模型中实验图2实验采用SIS病毒传播模型，在ER随机网络(a: N为104,=4)，BA无标度网络模型(b：N= 104,m0=8,m=4；c：N= 104,m0=8,m=6)中进行仿真。F为感染节点的密度，q为免疫节点的比例。第16页，共41页，2022年，5月20日，13点27分，星期二在

6、现实网络中实验图3实验采用SIS病毒传播模型在(autonomous system)AS层面的Internet 网络和High Energy Physics-Theory (HEP-Th)网络中测试局域搜索免疫的性能。F为感染节点的密度，q为免疫节点的比例第17页，共41页，2022年，5月20日，13点27分，星期二该免疫与聚类系数之间的关系由于局域搜索免疫是通过搜索邻居节点中度数最大的节点进行免疫，直观来讲该免疫的性能与网络的聚类系数有着某些联系 Assortative wiring 算法5能在保持节点度分布不变的前提下，增加网络的聚类系数。任意选择两条边，对两条边对应的四个顶点重新连接：

7、用一条边连接两个度数比较大的节点，另一条边连接两个度数比较小的节点。第18页，共41页，2022年，5月20日，13点27分，星期二图4 在BA无标度网络中，聚类系数与局域搜索免疫性能之间的关系。F为算法的免疫临界值，c为网络的聚类系数对BA无标度网络（N=104,m0=8,m=4）使用assortative wiring 算法对网络增加聚类系数第19页，共41页，2022年，5月20日，13点27分，星期二对于局域搜索免疫的改进局域免疫算法是随机选择一个节点，然后按照一定要求搜索。如果一个网络是由几个小的不连通的网络组成，那么这种策略就有可能一直在一个小的网络中进行循环搜索。解决方案：n种局

8、域搜索免疫同时进行 第20页，共41页，2022年，5月20日，13点27分，星期二改进的局域搜索免疫问题:n=?第21页，共41页，2022年，5月20日，13点27分，星期二删除边的免疫无论是熟人免疫还是目标免疫，基本思想都是找到度数大的节点进行免疫，也就相当于对度数大节点的所有的边进行删除，但是并不是所有的边都有必要删除的。比如节点i的度数很大，而节点j的度数很小，因为度数小的节点在疾病传播过程中起的作用很小，所以边E(i,j)也就没有必要删除。如果是通过物理的方式对网络进行免疫，那么对节点进行免疫，就极大的破坏了网络的连通度。第22页，共41页，2022年，5月20日，13点27分，星

9、期二连通度指的是两个随机选择的个体之间存在路径相连接的概率，其决定了网络的活跃性，可以通过宽度优先搜索算法6来计算。宽度优先搜索算法是一种图形搜索策略，从一个源节点开始搜索其邻居节点，然后搜索与邻居节点最近的节点，直到满足条件为止。第23页，共41页，2022年，5月20日，13点27分，星期二为了有效地降低感染节点的密度，并且提高网络的连通度，我们提出了删除边的免疫策略(Edges Cut Immunization Strategy, EC免疫策略)。首先是按照节点的度数进行排序，从高到低选择一定数目的节点，删除节点与节点直接相连的边。为了降低病毒在度数大节点之间的传播，也要删除边E(i,j

10、)，如果其余节点i具有多于一条边连接到给定数目节点j。 第24页，共41页，2022年，5月20日，13点27分，星期二第25页，共41页，2022年，5月20日，13点27分，星期二删除边的免疫第26页，共41页，2022年，5月20日，13点27分，星期二在模型中测试免疫策略性能图5实验采用SIS病毒传播模型，在ER随机网络(图a: N为104,=4)，BA无标度网络 (图b：N= 104,m0=8,m=4；图c：N= 104,m0=8,m=6)中进行仿真。F为感染节点的密度，q为免疫边的比例。第27页，共41页，2022年，5月20日，13点27分，星期二在模型中测试连通度图6 研究目标

11、免疫和EC免疫策略对网络模型连通度C的影响。其中q为免疫边的比例第28页，共41页，2022年，5月20日，13点27分，星期二在现实网络中测试免疫的性能图7基于SIS病毒传播模型，分别采用目标免疫和EC策略对于(a)AS网络,(b)HEP-Th网络和(c)PGP网络，进行免疫，根据删除边的比例q的变化研究感染节点概率F的变化。第29页，共41页，2022年，5月20日，13点27分，星期二在现实网络中测试连通度图8研究目标免疫和EC免疫策略对现实网络连通度C的影响。其中q为免疫边的比例第30页，共41页，2022年，5月20日，13点27分，星期二对于EC免疫策略的思考 EC免疫是从全局角度

12、来对边进行免疫，也同样可以从局部信息的角度来处理。关于边的免疫，一直感觉不是很切实际，毕竟在现实生活中，都是对整个节点进行免疫，比如某人患有H1N1，就把他完全隔离，并没有要求这个人只能见某些人或不能见某些人，所以对于EC免疫策略的实用性方面一直存在疑惑。第31页，共41页，2022年，5月20日，13点27分，星期二多重图形剖分免疫以往的免疫策略的免疫原则为：根据度数或者介数，对重要的节点进行免疫。Yiping Chen 通过对目标免疫分析发现：目标免疫策略把网络分成好几种小的网络。小的网络在病毒传播过程中起的作用很小，所以把网络分成好几个小的网络实际上浪费了代价。Yiping 通过嵌入分割

13、算法（nested dissection algorithm）8把网络分成几个近似大小的网络，然后对分割集团进行免疫，提出了EGP策（equal graph partitioning immunization strategy）。 第32页，共41页，2022年，5月20日，13点27分，星期二EGP免疫策略可以比目标免疫少用5%-50%的免疫剂量，达到相同的感染密度。原则是 免疫一组节点(separator group)，节点把网络分成几个相似大小的集团。 图9来自文献7第33页，共41页，2022年，5月20日，13点27分，星期二类似于EGP算法的策略，可以同样采用嵌入式分割算法，用边对

14、网络进行划分，提出了多重图形剖分算法。 图10 Nested dissection 的执行过程（取自文献8）第34页，共41页，2022年，5月20日，13点27分，星期二在linux环境下通过metis软件中的kmetis和pmetis程序来对网络划分，结果是把每个顶点对应的集团编号存放在文本中，然后对于不同集团之间的边进行免疫。实验如图11：图11第35页，共41页，2022年，5月20日，13点27分，星期二对于有向网络免疫的思考M. E. J. Newman的 Email networks and the spread of computer viruses9文章，对email有向网络

15、进行分析免疫，首先是把Email网络进行分析图12 Email的分析来自文献9第36页，共41页，2022年，5月20日，13点27分，星期二然后根据出度对于Email网络进行目标免疫图13 对Email网络进行免疫来自文献9第37页，共41页，2022年，5月20日，13点27分，星期二 针对有向网络的免疫，我思考的是使用类似于pagerank 算法来求解。 Pagerank的思想是对网页进行打分，原理：网页A指向网页B，则B=A的分值/A的出度+。针对SIS病毒传播模型，比如节点B,C,D三个节点指向节点A，那么节点A感染病毒的概率为至少有一个邻居节点为感染节点A=1-(1-B)(1-C)

16、(1-D)所以针对有向无权网络使用SIS病毒传播模型： 如果n与i有边连接, E(n,i)=1，否则为0。value(i)为节点i感染疾病的可能性问题：大型稀疏矩阵的求解第38页，共41页，2022年，5月20日，13点27分，星期二参考文献1 Reuven Cohen, Shlomo Havlin, Daniel ben.Avraham, Phys Rev Lett 91 (2003) 2779012 Pastor-Satorras R, Vespignani A, Phys. Rev. E. 65 (2002) 0361043 Madar, N.; Kalisky, T.; Cohen,

17、R.; Ben-Avraham, D，etc.Eur.Phys.J.B, 38(2004):269-2764 Stanley Milgram, The Small World Problem, Psychology Today, 1967, Vol. 2, 60-67 5 Shi Zhou, Raul J. Mondragon, New Journal of Physics 9(2007) 1736 Andy Yoo, Edmond Chow, etc,Proceedings of the 2005 ACM/IEEE conference on Supercomputing, 2005.7 Chen Y, Paul G, etc, Phys Rev Lett. 101(2008) 058701.8 Bruce Hendrickson, Robert Leland， A multilevel algorithm for partitioning graphs，Supercomputing，Tech. report SAND93-1301, Sandia National