一种基于动态网络标志物的区域性传染病疫情预警方法技术

本发明专利技术公开了一种基于动态网络标志物的区域性传染病疫情预警方法，该方法使用正规化处理之后的区域平均就诊数据作为计算数据，根据各个区域的地理分布信息和区域间的邻接关系构建城市网络；然后提出景观动态网络标志物方法来计算每个区域对应的L‑DNM得分；最后根据不同时刻的L‑DNM得分，使用合理阈值判断出传染病疫情的预警信号。相对于传统方法中使用的就诊人数数据直接判断，本发明专利技术所使用的景观动态网络标志物方法将传染病的传播与爆发分为三个阶段，即正常阶段、预爆发阶段和爆发阶段。与临床就诊量陡然增加的爆发阶段不同，预爆发阶段与正常阶段没有显著差异。另外，本发明专利技术提供多种数据可视化形式，直观展现疫情中的预警信号。

【技术实现步骤摘要】
一种基于动态网络标志物的区域性传染病疫情预警方法
本专利技术涉及数据挖掘
，具体涉及一种基于动态网络标志物(DynamicNetworkMarker，DNM)的区域性传染病疫情预警方法。
技术介绍
随着经济全球化的发展和人员流动的增多，传染性疾病的传播和蔓延变得更快更广，且难以控制。世界各国面临着新老传染病的双重威胁，少数传染病将被消灭，已经基本控制的传染病又卷土重来，新传染病不断被陆续发现；近年有多起重大传染性疾病的爆发流行，如SARS的爆发流行、甲型肝炎流行和血吸虫病的再次增多，特别是最近发生的新型冠状病毒引发的肺炎疫情。可见，不管对于常规季节性传染病还是突发性新型传染病，建立一个基于城市网络的监测预警系统非常重要。有效的传染病疫情检测方法通过持续而系统性地收集、分析和解释健康数据，有可能尽早探测到爆发，有助于疾病的预防和控制。
技术实现思路
本专利技术的目的是为了解决现有技术中的上述缺陷，提供一种基于动态网络标志物的区域性传染病疫情预警方法，并且最终识别出传染病传播过程中的分级预警信号。本专利技术的目的可以通过采取如下技术方案达到：一种基于动态网络标志物的区域性传染病疫情预警方法，给定区域邻接表数据和区域对应的疫情确诊人数数据，构建区域性城市网络来识别疫情爆发信号，该预警方法包括如下步骤：S1、构建区域性城市网络。具体地，根据各个区域的地理分布信息和区域间的邻接关系构建有权无向图。在图中，每一个节点代表一个地理区域(如县、市、省等)，图中的边代表两个节点的邻接关系。两个节点通过一条边相连，则表示两个节点可以直接通达。S2、将传染病监测数据映射到城市网络图中。具体地，每一个区的传染病确诊人数即为对应节点的值。随着时间的推移，构建的城市网络将成为一个动态变化的网络图，其中，节点大小和边的权重都在时刻变化。S3、计算景观动态网络标志物得分，简称DNMScores。具体地，首先将城市网络切分为若干个小的局部网络。每一个局部网络包含一个中心节点和它的所有一阶邻居。对于有n个节点的局部网络，即一个中心节点和n-1个一阶邻居节点，在t时刻局部网络的得分It计算如下：It＝|ΔSDt(in)|[|ΔPCCt(in)|+|ΔPCCt(out)|]其中，|ΔSDt(in)|是局部网络中所有节点标准差的绝对值，详细的计算公式如下：其中，n表示局部网络内的节点数量，SDt(i)表示局部网络中的节点i在t时刻的标准差；|ΔPCCt(in)|是局部网络内节点i和节点j的皮尔逊相关系数，详细的计算公式如下：其中，n表示局部网络一共有n个节点，PCCt(i,j)表示局部网络内节点i和节点j在t时刻的皮尔逊相关系数；|ΔPCCt(out)|是局部网络内的一个节点x与局部网络外一个节点y的皮尔逊相关系数。详细的计算公式如下：其中，n表示局部网络一共有n个节点，PCCt(x,y)表示局部网络内的一个节点x与局部网络外一个节点y在t时刻的皮尔逊相关系数。值得注意的是，在以上计算过程中，以滑动窗口的方式来处理正规化后的数据，即基于每五周滑动窗口大小的传染病确诊数据来计算标准差和皮尔逊相关系数。S4、基于DNM得分判断出多级预警信号。详细地，设置一个合适的阈值Mthreshold作为预警信号的判断依据，当达到相应的阈值时发出相应级别的预警。对于不同的疾病，阈值的设置不尽相同，即传染性强的传染病阈值设置要低。阈值Mthreshold的计算公式如下：其中，It(i)代表在时间t，n个节点的城市网络中节点i的得分It。本专利技术相对于现有技术具有如下的优点及效果：1)本专利技术使用正规化的数据(该地区确诊总人数/该地区定点医院数量)作为计算原始数据，因为每个地区的人口数与诊所的数量大致成正比。2)本专利技术方法将传染病的传播与爆发分为三个阶段，即正常阶段、预爆发阶段和爆发阶段。与临床就诊量陡然增加的爆发阶段不同，预爆发阶段与正常阶段没有显著差异。同时因为复杂的生物系统和社会系统，一般方法预测出预爆发阶段是非常有难度的。DNM方法则可以很好地克服以上难点，识别出预爆发阶段，给出传染病疫情的预警信号。3)本专利技术根据区域的地理位置信息和区域间的邻接关系构建出宏观的城市网络。传统方法中只是根据单一地区的就诊人数来判断，忽略了传染病的传播能力和地区间的相关信息。使用城市网络带来的增益信息可以克服样本量小的问题，并通过局部网络的重叠性修正测不准的局部网络。4)相比于传统的预警信号给出标准，本专利技术采用全新的阈值计算方式，这种方法包含了当前时刻及其之前的所有时刻的数据，大大提高了预警信号的置信度。附图说明图1是本专利技术公开的一种基于景观动态网络标志物的区域性传染病疫情预警方法的流程图；图2是本专利技术实施例中日本东京都23个区的城市网络图；图3是本专利技术实施例中设计滑动窗口处理的计算方法示意图；图4是本专利技术实施例中2012年和2017年的预警信号结果图；图5是本专利技术实施例中公开的方法与机器学习方法的效果对比图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。实施例本实施例是以收集到的日本东京都23个区2009-2018年的手足口病就诊数据集为原始数据，利用基于动态网络标志物的区域性传染病疫情预警方法来识别每一年份的预警信号。本实施例不仅指导说明本专利技术的方法步骤，同时佐证方法的有效性。原始数据收集自日本东京都传染病监测中心官方网站。图1是本专利技术的基于动态网络标志物的区域性传染病疫情预警方法的流程图，具体包括以下步骤：S1、如图1的S1.1中，输入区域邻接表，并根据输入的区域邻接表构建区域性城市网络。邻接表是存储无向图的常用结构。邻接表是一个二维容器，第一维描述某个节点，第二维描述该点所对应的邻接节点。本实施例中，根据输入的区域邻接表就可以得到每个子区域与其他子区域的邻接关系，从而构建出一个无向图，即区域性城市网络。无向图中的每个节点与某个子区域对应，无向图中相邻的两个节点，表示实际地理位置相邻的两个子区域。详细地，如图2所示，举例说明，minato区对应编号1，相关邻接表对应的信息为{…,1:[2，3，4，5，6],…},即编号1的区域与编号2、3、4、5、6五个区域是相邻的。以此类推，最终构建出一个23个节点的有权无向图。详细的23个节点信息见图2。如图1的S1.2中输入各区域疫情数据，并对疫情数据进行正规化处理。由于每个地区的人口数与诊所的数量大致成正比，所以用该地区就诊总人数除以该地区定点医院数量，得到的该地区平均就本文档来自技高网...

【技术保护点】
1.一种基于动态网络标志物的区域性传染病疫情预警方法，其特征在于，所述的区域性传染病疫情预警方法包括下列步骤：/nS1、输入区域邻接表，并根据输入的区域邻接表构建区域性城市网络，同时，输入各区域疫情数据，并对疫情数据进行正规化处理；/nS2、将经过正规化的各区域疫情数据映射到所述的区域性城市网络中；/nS3、对于映射后的区域性城市网络，采用景观动态网络标志物方法计算景观动态网络标志物得分，简称DNM Scores；/nS4、根据计算得到的DNM Scores，判断当前DNM Scores是否达到阈值标准，若达到则发出预警信号；反之则继续到下一时刻，重复执行步骤S3和步骤S4继续判断并识别预警信号。/n

【技术特征摘要】
1.一种基于动态网络标志物的区域性传染病疫情预警方法，其特征在于，所述的区域性传染病疫情预警方法包括下列步骤：
S1、输入区域邻接表，并根据输入的区域邻接表构建区域性城市网络，同时，输入各区域疫情数据，并对疫情数据进行正规化处理；
S2、将经过正规化的各区域疫情数据映射到所述的区域性城市网络中；
S3、对于映射后的区域性城市网络，采用景观动态网络标志物方法计算景观动态网络标志物得分，简称DNMScores；
S4、根据计算得到的DNMScores，判断当前DNMScores是否达到阈值标准，若达到则发出预警信号；反之则继续到下一时刻，重复执行步骤S3和步骤S4继续判断并识别预警信号。


2.根据权利要求1所述的一种基于动态网络标志物的区域性传染病疫情预警方法，其特征在于，所述的步骤S1中根据输入的区域邻接表构建区域性城市网络的过程如下：
根据输入的区域邻接表得到每个子区域与其他子区域的邻接关系，从而构建出一个无向图，即区域性城市网络，该无向图中每个节点与一个子区域对应，无向图中相邻的两个节点，表示实际地理位置相邻的两个子区域。


3.根据权利要求1所述的一种基于动态网络标志物的区域性传染病疫情预警方法，其特征在于，所述的步骤S1中对疫情数据进行正规化处理的过程如下：
将该地区就诊总人数除以该地区定点医院数量得到该地区平均就诊人数，该地区平均就诊人数即为正规化处理后的数据。


4.根据权利要求1所述的一种基于动态网络标志物的区域性传染病疫情预警方法，其特征在于，所述的步骤S2中将正规...

【专利技术属性】
技术研发人员：张旭航，刘锐，陈培，钟佳元，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：广东;44