CN106789160B - 一种基于时间序列的复杂网络链接预测方法 - Google Patents
一种基于时间序列的复杂网络链接预测方法 Download PDFInfo
- Publication number
- CN106789160B CN106789160B CN201611051368.6A CN201611051368A CN106789160B CN 106789160 B CN106789160 B CN 106789160B CN 201611051368 A CN201611051368 A CN 201611051368A CN 106789160 B CN106789160 B CN 106789160B
- Authority
- CN
- China
- Prior art keywords
- node
- link
- probability
- time series
- prediction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/147—Network analysis or design for predicting network behaviour
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明公开了一种基于时间序列的复杂网络链接预测方法，包括如下步骤：步骤一、使用静态方法对任意两个节点的链接概率进行评分；步骤二、使用时间序列的方法，获取最近不相连的任意两节点的链接出现概率；步骤三、对于使用静态方法计算出的链接出现概率序列p1，p2，…pk，计算其权值：W(i)＝α(pi‑pi‑1)+(1‑α)W(i‑1)i＝2…k；W(1)＝1，得到最终的链接出现的概率pk＝W(k)pk+1；步骤四、对步骤三中计算得到的结果进行排序，输出预测的链接。本发明提供的基于时间序列的复杂网络链接预测方法，通过加权算法，使预测结果更为精确。
Description
技术领域
本发明属于复杂网络链接预测技术领域，特别涉及一种基于时间序列的复杂网络链接预测方法。
背景技术
复杂网络模型是对各种真实存在网络的抽象，比如社交网络、科研合作网络、生物代谢网络等。链接预测，是根据观察到的网络中已有的链接和节点属性信息，来估计两个节点间链接存在的可能性，包括当前存在但是观察不到的链接，以及将来可能出现的链接。作为链接挖掘分析中最重要的问题之一，链接预测有着多种多样的应用：它可以用于推荐系统，帮助人们找到新朋友或潜在的合作者，在网上购物中提供感兴趣的商品；也能用来推断完整的网络结构，更好地理解网络演化；在生物学领域，还能用来发现蛋白质之间的交互情况。
现有的链接预测可以分为两类：静态网络链接预测以及基于时间序列的网络连接预测。
静态方法主要采用网络拓扑信息以及节点属性信息进行链接预测。
时间序列预测方法主要采用不同时间点的网络变化信息对链接进行预测。预测的网络变化需要满足线性变化的条件。
两种方法对任意两个不直接相连节点之间存在链接的概率进行评分，最后通过评分阈值或整体预测比例进行确定存在链接的节点对。然而这种方法对以下两种情况不能很好的预测：第一种是两个节点之间距离较近并且节点较为相似但是网络比较稳定的情况(即：预测出现不应出现的链接)；第二种是两节点之间连接概率随着时间点的推移不断增大但是当前结果并没有达到阈值的情况(即：不能预测应该出现的链接)。
发明内容
本发明的目的是解决现有网络链接预测方法中预测出现不应出现的链接和不能预测应该出现的链接的缺陷，提供了一种基于时间序列的复杂网络链接预测方法。
本发明提供的技术方案为：
一种基于时间序列的复杂网络链接预测方法，包括如下步骤：
步骤一、使用静态方法对任意两个节点的链接概率进行评分；
步骤二、使用时间序列的方法，获取最近不相连的任意两节点的链接出现概率；
步骤三、对于使用静态方法计算出的链接出现概率序列p1，p2，…pk，计算其权值：
W(i)＝α(pi-pi-1)+(1-α)W(i-1)i＝2…k
W(1)＝1
其中α为调整系数，k为网络快照个数；
得到最终的链接出现的概率P＝W(k)pk+1；
步骤四、对步骤三中计算得到的结果进行排序，输出预测的链接。
优选的是，步骤一中，任取节点A和节点B，节点A到节点B的链接概率为
其中为从节点A到节点B的第i条无环路径，为这条路径的长度。
优选的是，步骤一中，节点A到节点B的链接概率评分为：
εt是独立同分布的随机变量序列，满足期望为0的条件，Φi是不同时刻网络快照的权重。
本发明的有益效果是：本发明提供了一种于时间序列的复杂网络链接预测方法，通过加权算法，使预测结果更为精确。
附图说明
图1为本发明所述的基于时间序列的复杂网络链接预测方法流程图。
图2为第一时刻节点连接示意图。
图3为第二时刻节点连接示意图。
图4为第三时刻节点连接示意图。
图5为第四时刻节点连接示意图。
具体实施方式
下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
如图1所示，本发明提供了一种基于时间序列的复杂网络链接预测方法，包括如下步骤：
步骤一：使用静态方法对任意两个节点的链接概率进行评分。在该阶段可以使用网络拓扑，节点属性等信息对。
步骤二：使用时间序列的方法针对所有的静态链接概率评分结果进行计算，得出最近不相连的任意两节点的链接出现概率。
第三步：计算权值。具体方法如下：
对于一个使用静态方法计算出的一个链接出现概率序列p1，p2，…pk，其中k表示时间序列的中时间步的个数，他的权值可以使用以下递归公式计算出来：
W(i)＝α((pi-pi-1)+(1-α)W(i-1)i＝2…k
W(1)＝1
对于最终的链接出现的概率计算加权得到：P＝W(k)pk+1
步骤四：排序所有结果，根据定义好的阈值或比例系数输出预测的链接。
其实现过程如下：
静态方法中使用了拓扑信息，节点A到节点B的链接概率为
其中为从节点A到节点B的第i条无环路径。为这条路径的长度。在时间序列方法中使用了ARMA方法，得到概率评分：
如图2-图5所示，通过静态方法获得的评分如表1-表4所示。
表1：通过静态方法计算出的图2中两个节点链接概率评分
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
| 1 | NA | NA | NA | NA | 2/9 | 5/16 | NA | 1/8 |
| 2 | NA | 5/16 | NA | 5/16 | 2/9 | 1/8 | 34/225 | |
| 3 | NA | NA | 1/2 | NA | 5/16 | 2/9 | ||
| 4 | NA | NA | 1/2 | 2/9 | 5/16 | |||
| 5 | NA | NA | 5/16 | NA | ||||
| 6 | NA | NA | 5/16 | |||||
| 7 | NA | 34/225 | ||||||
| 8 | NA |
表2：通过静态方法计算出的图3中两个节点链接概率评分
表3：通过静态方法计算出的图4中两个节点链接概率评分
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
| 1 | NA | NA | NA | NA | 118/225 | 11/16 | NA | 11/36 |
| 2 | NA | 49/72 | NA | 5/9 | 118/225 | 19/36 | 3439/11025 | |
| 3 | NA | NA | 49/72 | NA | 49/72 | 3782/11025 | ||
| 4 | NA | NA | 11/16 | 118/225 | 57/144 | |||
| 5 | NA | NA | 5/9 | NA | ||||
| 6 | NA | NA | 57/144 | |||||
| 7 | NA | 3439/11025 | ||||||
| 8 | NA |
表4：通过静态方法计算出的图5中两个节点链接概率评分
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
| 1 | NA | NA | NA | NA | 3013/3600 | 131/144 | NA | 411/900 |
| 2 | NA | 1319/1800 | NA | 1061/1200 | 697/900 | 683/900 | 80749/176400 | |
| 3 | NA | NA | 523/600 | NA | 1447/1800 | 81337/176400 | ||
| 4 | NA | NA | 131/144 | 3013/3600 | 571/1200 | |||
| 5 | NA | NA | 56/75 | NA | ||||
| 6 | NA | NA | 571/1200 | |||||
| 7 | NA | 25283/58800 | ||||||
| 8 | NA |
通过时间序列的方法获得的评分如表5所示，其中其中εt取值为0，Φi＝{0.5，0.3，0.16，0.04}。
表5：通过时间序列方法得出的最终评分
根据时间序列的方法获得的链接可能性排列为：<4，6>,<1，6>,<3，5>,<2，5>,<3，7>,<2，3>,<4，7>,<5，7>,<2，6>,<2，7>,<6，8>,<4，8>,<3，8>,<1，8>,<2，8>,<1，5>。
权值为：
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
| 1 | NA | NA | NA | NA | 0.303 | 0.251 | NA | 0.185 |
| 2 | NA | 0.180 | NA | 0.316 | 0.213 | 0.291 | 0.098 | |
| 3 | NA | NA | 0.212 | NA | 0.190 | 0.260 | ||
| 4 | NA | NA | 0.234 | 0.309 | 0.128 | |||
| 5 | NA | NA | 0.227 | NA | ||||
| 6 | NA | NA | 0.128 | |||||
| 7 | NA | 0.098 | ||||||
| 8 | NA |
最终的评分为：
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
| 1 | NA | NA | NA | NA | 0.106656 | 0.191764 | NA | 0.066372 |
| 2 | NA | 0.11447 | NA | 0.213548 | 0.126735 | 0.12179 | 0.034695 | |
| 3 | NA | NA | 0.159 | NA | 0.171402 | 0.10006 | ||
| 4 | NA | NA | 0.180414 | 0.196833 | 0.054144 | |||
| 5 | NA | NA | 0.139151 | NA | ||||
| 6 | NA | NA | 0.05619 | |||||
| 7 | NA | 0.03450 | ||||||
| 8 | NA |
加权之后的链接可能性排列为：<2，5>,<4，7>,<1，6>,<4，6>,<3，7>,<3，5>,<5，7>,<2，6>,<2，7>,<2，3>,<1，5>,<3，8>,<1，8>,<6，8>,<4，8>,<2，8>,<7，8>。
和之前的序列相比，节点2和节点5的链接概率以及节点4和节点7的链接概率被提升到最前方，而节点4和节点6的链接概率以及节点1和节点6的链接概率在序列中被后移。通过网络演化可看出，我们解决了前文提及的两个问题。
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。
Claims (1)
1.一种基于时间序列的复杂网络链接预测方法，其特征在于，包括如下步骤：
步骤一、使用静态方法对任意两个节点的链接概率进行评分；任取节点A和节点B，节点A到节点B的链接概率为
其中为从节点A到节点B的第i条无环路径，为这条路径的长度；
节点A到节点B的链接概率评分为：
εt是独立同分布的随机变量序列，满足期望为0的条件，Φi是不同时刻网络快照的权重；为第i时刻网络快照中的从节点A到节点B的链接概率；
步骤二、使用时间序列的方法，获取最近不相连的任意两节点的链接出现概率；
步骤三、对于使用静态方法计算出的链接出现概率序列p1，p2，…pk，计算其权值：
W(i)＝α(pi-pi-1)+(1-α)W(i-1)；i＝2…k
W(1)＝1
其中α为调整系数，k为网络快照个数；
得到最终的链接出现的概率P＝W(k)pk+1；
步骤四、对步骤三中计算得到的结果进行排序，输出预测的链接。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611051368.6A CN106789160B (zh) | 2016-11-24 | 2016-11-24 | 一种基于时间序列的复杂网络链接预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611051368.6A CN106789160B (zh) | 2016-11-24 | 2016-11-24 | 一种基于时间序列的复杂网络链接预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106789160A CN106789160A (zh) | 2017-05-31 |
| CN106789160B true CN106789160B (zh) | 2019-11-22 |
Family
ID=58912803
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611051368.6A Expired - Fee Related CN106789160B (zh) | 2016-11-24 | 2016-11-24 | 一种基于时间序列的复杂网络链接预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106789160B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107730055A (zh) * | 2017-11-21 | 2018-02-23 | 燕山大学 | 符号网络链接预测方法及终端设备 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104408072A (zh) * | 2014-10-30 | 2015-03-11 | 广东电网有限责任公司电力科学研究院 | 一种基于复杂网络理论的适用于分类的时间序列特征提取方法 |
-
2016
- 2016-11-24 CN CN201611051368.6A patent/CN106789160B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104408072A (zh) * | 2014-10-30 | 2015-03-11 | 广东电网有限责任公司电力科学研究院 | 一种基于复杂网络理论的适用于分类的时间序列特征提取方法 |
Non-Patent Citations (3)
| Title |
|---|
| link prediction using time series of neighborhood based node similarity scores;Ismail Gunes等;《Data Mining and Knowledge Discovery》;20160131;第30卷(第1期);第153-158页 * |
| The time-series link prediciton problem with applications in communication surveillance;Zan Huang 等;《INforms journal on computing》;20090228;第21卷(第2期);第287-293页，图1 * |
| Time sreies based link prediction;da Silva Soares等;《Neural Networks，the 2012 international joint conference on IEEE》;20121231;第4页第1栏 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106789160A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liang et al. | Differential evolution based on fitness Euclidean-distance ratio for multimodal optimization | |
| CN102591915B (zh) | 一种基于标签迁移学习的推荐方法 | |
| CN107391542A (zh) | 一种基于文件知识图谱的开源软件社区专家推荐方法 | |
| CN103399858A (zh) | 基于信任的社会化协同过滤推荐方法 | |
| Zhang et al. | Alleviating new user cold-start in user-based collaborative filtering via bipartite network | |
| Chen et al. | Opinion Dynamics of Social‐Similarity‐Based Hegselmann–Krause Model | |
| Xiong et al. | Dynamic evolution of collective emotions in social networks: a case study of Sina weibo | |
| Zhang et al. | The empirical likelihood for first-order random coefficient integer-valued autoregressive processes | |
| Lu et al. | Real‐Time Corrected Traffic Correlation Model for Traffic Flow Forecasting | |
| CN108573062A (zh) | 一种基于异构社交关系的个性化推荐方法 | |
| Liu et al. | Power grid partition with improved biogeography-based optimization algorithm | |
| Zou et al. | Collaborative community-specific microblog sentiment analysis via multi-task learning | |
| CN106789160B (zh) | 一种基于时间序列的复杂网络链接预测方法 | |
| CN103942298A (zh) | 基于线性回归的推荐方法及系统 | |
| CN107016459A (zh) | 一种基于网络社区信息的点到点最短路径计算方法 | |
| Church et al. | Voltage collapse detection using ant colony optimization for smart grid applications | |
| Huang et al. | A projection neural network with mixed delays for solving linear variational inequality | |
| CN109299849B (zh) | 一种社会网络中群体需求层次计算方法 | |
| Greco et al. | Metaheuristics for combinatorial optimization | |
| Fu et al. | Preferential information dynamics model for online social networks | |
| Sawai | Exploring a new small-world network for real-world applications | |
| Xin et al. | A graph neural network-based algorithm for Point-of-Interest recommendation using social relation and time series | |
| Ho et al. | Discovering communities of users on social networks based on topic model combined with Kohonen network | |
| CN110717085B (zh) | 一种基于虚拟品牌社区的意见领袖识别方法 | |
| Sayyadi et al. | A new view to uncertainty in Electre III method by introducing interval numbers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191122 Termination date: 20201124 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |