Abstract
Link prediction in networks has applications in computer science, graph theory, biology, economics, etc. Link prediction is a very well studied problem. Out of all the different versions, link prediction for unipartite graphs has attracted most attention. In this work we focus on link prediction for bipartite graphs that is based on two very important concepts—potential energy and mutual information. In the three step approach; first the bipartite graph is converted into a unipartite graph with the help of a weighted projection, next the potential energy and mutual information between each node pair in the projected graph is computed. Finally, we present Potential EnergyMutual Information based similarity metric which helps in prediction of potential links. To evaluate the performance of the proposed algorithm four similarity metrics, namely AUC, Precision, Predictionpower and Precision@K were calculated and compared with eleven baseline algorithms. The Experimental results show that the proposed method outperforms the baseline algorithms.
Introduction
The applications of link prediction in realworld networks has been attracting the attention of researchers from various domains. The real world networks can be represented with the help of graphs, where nodes represent entities that are connected by edges. The edges represent associations or interactions between nodes. These networks are dynamic in nature, where new nodes and edges can be observed at future timestamps. Link prediction task aims to identify the most probable links that may appear in the near future. Some of the popular realworld applications of link prediction are: Prediction of hidden relationships between terrorists, ecommerce recommendation system, prediction of drug side effects, proteinprotein interaction prediction, finding missing reactions in metabolic networks, predicting coauthorship. Link prediction problem takes as input a snapshot of the network and computes the likelihood of a connection between two unconnected vertices^{1,2,3}. The algorithms not only predicts the addition of new edges, it also predicts the connections that may disappear from the network in future. Most of the link prediction algorithms have been designed for unipartite networks. The structure and properties of bipartite networks are different from unipartite ones. For this reason the algorithms that perform well for unipartite networks may not work well for bipartite networks. Link prediction on Bipartite networks have applications in Recommendation systems. Li et al.^{4} used applied link prediction to build recommender systems based on the concept of graph kernels. They generated random walks starting from a focal useritem pair in the graph kernel. Due to the fact that defining features in the complex graph is challenging, authors used a kernelbased machine learning framework that works on kernel function and this function defines the similarity between data instances. The authors performed testing on three realworld datasets and the it was observed that the proposed method outperformed the baseline algorithms. Kurt et al.^{5} also proposed a graphbased recommender system (SIMLP) that computed similarity between two nodes in a bipartite network. They assigned complex numbers as the weight of an edge connecting a node pair. SIMLP algorithm was tested on two datasets and the results show that it algorithm performed better than the complex numberbased baseline algorithms.
We have used the terms graphs and networks interchangeably throughout the paper. The Link Prediction problem can formally be defined as follows.
Problem statement Given a bipartite graph G(U, V, E), where U and V are two disjoint and independent sets of vertices and E is set of edges at timestamp \( T_1 \). The link prediction algorithm aims to predict new links among nodes that can be observed at timestamp \( T_2 \), such that \( T_2 > T_1 \). These links may be due to one of the two reasons; either due to the formation of a new connection, or due to revival of an existing missing connection.
In this work, the authors make the following contributions in order to propose new algorithm for the above mentioned problem:

1.
Two parameters have been introduced with respect to link prediction: potential energy and mutual information.

2.
Potential Energy and Mutual Information based similarity metric (PMIS) proposed to compute the weight of the patterns.

3.
The proposed algorithm works on each disconnected node pair in the network instead of the only candidate node pair.

4.
The proposed algorithm evaluates on ten real world datasets and demonstrates superior performance of the proposed algorithm compared to baseline link prediction techniques.
Link prediction algorithms^{6,7,8,9,10} usually follow one of the following four approaches: Node Based link prediction, Neighborbased link prediction, Pathbased and Social Theorybased link prediction. LibenNowell et al.^{1} analysed various proximity measures and suggested proximity measures for best prediction results among node neighborhoods approach, pathbased approach and Higherlevel (meta) approaches. Common neighbors, Jaccard coefficient, Adamic/Adar coefficient are common neighbor based approaches. Triangleclosing model states that nodepairs with a high number of common neighbors try to form a triangle in the graph. Such nodepairs have a high probability of forming connections in the future. This concept does not apply to bipartite networks because of their distinct topological structure. Thus common neighbor based approaches can not apply directly on bipartite networks. Hasan et al.^{11,12} further extended the work of LibenNowell et al.^{1}. Link prediction for bipartite networks has been addressed by many researchers^{4,13,14,15,16,17,18,19}. Cannistraci et al.^{20} proposed an algorithm to target not only common neighbors and neighbor’s common neighbors but also to their connections structure. This was the first attempt of bipartite formulations of the Common Neighbor index. Further, the authors in^{21} used the concept of local community paradigm (LCPtheory) for link prediction and states that the cohort of CNs and their cross interactions form a localcommunity edge. The cross interactions between common neighbors are called local community edges or links. The LCPbased method presented in^{21} improves topological prediction in bipartite complex networks^{13} used the concept of projection and supervised learning for link prediction. They introduced three link prediction metrics for bipartite graphs and implemented these metrics on DBLP dataset.Baltakiene et al. in^{22} presented the concept of entropy and used the Bipartite Configuration Model (BiCM)^{23} as a score function for predicting links. As used in Statistical Mechanics, probability per graph can be derived by maximizing the Shannon entropy under the constraint of the degree sequence. So by using the concept of maximal entropy they significantly improved the performance of link prediction. Gao et al. used the concept of projection and CNP(candidate node pairs) for link prediction in the bipartite networks^{24}. The authors first converted the bipartite graph into unipartite graph and then computed the CCNP (Connectivity of candidate node pairs) with the help of pattern weights for link prediction. They evaluated the performance of the algorithm on three datasets and experienced better results than baseline predictor(CN, Katz, ILP^{18}) on AUC. Shakibian et al. proposed another similarity measure based on mutual information and metapath in heterogeneous networks^{25}. They presented a framework in which link entropy is characterized as a semantic measure for link prediction. To measure the effectiveness of the algorithm, authors compared it with different classes of link prediction algorithms namely, mainstream metapath based link predictors, effective pathbased homogeneous link predictor, and LCPbased link indicators. The authors analysed the performance of algorithm on DBLP network. In 2019 Serpil et al. used strengthened projection technique for link prediction in evolving bipartite graph^{26}. They tried to predict link in large scale bipartite networks. Authors did link prediction in mainly two steps. In the first steps, they extracted potential link set and in the second step, they computed the prediction score of each potential link. For this purpose authors proposed a time aware proximity measure based on network evolution. For the result analysis, they used AUC and precision metrics. In the experiment, the authors compared his method with four baseline algorithms (AA, CN, JA, PA). And they found that his algorithm outperforms the baseline algorithm. It is often observed that the complete information about realworld complex networks is not available. Some examples of such networks are financial networks, social networks and biological networks. Cimini et al.^{27} presented an exhaustive review on statistical physics based approaches to predict statistically significant patterns in complex networks. They also addressed the reconstruction of network structure in the absence of complete information. Boguna et al.^{28} presented a review on network based approaches that very effectively identify both physical properties and mathematical properties that are fundamental to networks. They discussed three approaches and proposed interesting future directions. In the review they presented that in case of heterogeneous networks, the models based on hyperbolic space could be better than Euclidean space. In addition the hyperbolic space could be used for link prediction.
Preliminaries
Projection of bipartite graph
Projection of bipartite graphs can be used to convert it into unipartite graphs. For a given bipartite graph \( G =( U,V,E) \), its Uprojected graph can be represented as a unipartite graph \( G_u =(U,E_u) \) in which \((A, B) \in E_u \) if A and B have at least one neighbor common in G. This means, \( N(A) \cap N(B)\ne \phi \) and \( E_u \) can be described as follows:
Here two projections of the graph can be taken, one for U and another for V. Vprojection can be defined for the graph similar to its Uprojection. For a graph G, it’s Vprojected graph will be \( G_v=(V, E_v) \). For example Fig. 1a denotes the graph G where values of U and V are \(\{A,B,C,D,E\}\) and \(\{h,i,j,k,l\}\) respectively, Fig. 1b denotes their Uprojection and Fig. 1c denotes the Vprojection.
Pattern and pattern covered by a Node pair
Pattern Suppose A and B are two nodes in a bipartite graph G(U, V, E) and \(\{A, B\} \in U\). Then (A, B) forms a pattern if they have at least one common neighbor in V^{24}. Therefore a link (A, B) will be available in the projected graph. We can also say that each link \((A,B) \in E_u \) in the projected graph represents a pattern in bipartite graph G.
Pattern covered by a Node pair Suppose (A, i) be a node pair in bipartite graph G and \(G_u\) be the projected graph of G. For each node \( C \in N_u (A) \cap N (i) \), we call {A, C} a pattern covered by node pair(A, i)^{24}. A node pair may cover one or more patterns in the projected graph. Pattern covered by a node pair simply says that a similar edge has already existed in bipartite graph G. The more patterns a node pair covers, more are the chances that this node pair will be connected in the future. In this manner, the number of patterns secured by a node pair can be utilized to measure the likelihood of its edge presence.
Mutual information for link prediction
Our proposed algorithm uses the concept of mutual information and potential energy for link prediction in the bipartite graph.
Selfinformation: “Let X be a random variable and x be an outcome of X with probability p(x). Then, the selfinformation of x quantifies the uncertainty of the outcome x and is defined as follows”^{29}:
Mutual information: “Let X and Y be two random variables and x and y be their outcomes, respectively. The mutual information of X and Y measures the amount of reduction in uncertainty of the outcome x when the outcome y is known, or vice versa, and is defined as follows”^{30}:
Let x, y represent the two nodes in a graph and \(\Gamma (x), \Gamma (y)\) represent their set of neighbors. Also, common neighbors of x, y is represented by \(O_{xy}\). So we can say that \(O_{xy}=\Gamma (x)\cap \Gamma (y)\). Now for given node pair(x, y) and common neighbors \(O_{xy}\); The likelihood score of node pair(x, y) can be computed by the following equation^{31}.
Here \(I(L^{1}_{xy}\mid O_{xy})\) is the conditional selfinformation of the existence of an edge between node pair (x, y) when they have common neighbors \(O_{xy}\). On analysis of selfinformation of the nodepair, it is found that the smaller \(I(L^{1}_{xy}\mid O_{xy})\) is, the higher the probability of the existence of an edge. According to Eq. (2), we can derive the value of \(I(L^{1}_{xy}\mid O_{xy})\) as follows:
where \(I(L^{1}_{xy})\) represents the selfinformation of that node pair (x,y) that are already connected. \(I(L^{1}_{xy};O_{xy})\) represents the mutual information between node pair (x,y) that has one link between them and the node pair’s common neighbors are known. Now let’s consider that the elements of \(O_{xy}\) are independent of each other, then we can find the value of \(I(L^{1}_{xy};O_{xy})\) as follows:
Now \(I(L^{1}_{xy};z)\) can be calculated by \(I(L^1;z)\). \(I(L^1;z)\) is defined as the mean mutual information over all node pairs connected to node z.
Now we can find the value of \(I(L^{1}_{mn};z)\) with the help of Eq. (2).
Here \(I(L^{1}_{mn})\) indicates the selfinformation of node pair (m, n) is connected. \(I(L^{1}_{mn}\mid z)\) is simply the conditional information of that connected node pair (m, n) when node z is one of their common neighbors.
Now \(I(L^{1}_{mn}\mid z)\) can be calculated by the clustering coefficient of node z and the clustering coefficient of z can be calculated as follows:
where \(C_z\) represents the clustering coefficient of z. \(t_z\) represents the number of triangles passing through node z and \(d_z\) represents the degree of node z. Since we have value of \(p(L^{1}_{mn}\mid z)\) we can easily find the value of \(I(L^{1}_{mn}\mid z)\).
\(p(L^{1}_{mn})\) can be computed with the help of \(p(L^{0}_{mn})\). \(L^{0}_{mn}\) represents the event that there is no edge that exists between node m and node n. It is considered here that no degree correlation exists. The value of \(p(L^{0}_{mn})\) can be calculated with the help of path entropy.
Here \(d_m\) and \(d_n\) are the degree of node m and n, respectively. \(T_l\) is the total no of edges in the graph. This formula is symmetric. Thus
So now \(p(L^{1}_{mn})\) and \(p(L^{1}_{nm})\) can be calculated as follows.
With the help of Eq. (1) we can find the value of \(I(L^{1}_{mn})\) and \(I(L^{1}_{nm})\). Collecting these results, we can get the following things.
So with the help of the previous derivation, we have
Here chances of the existence of an edge between node pair(x, y) are directly proportional to \(S^{MI}_{xy}\). It simply means higher the \(S^{MI}_{xy}\), the more likely the nodes will be connected.
Methods
In this section, we are presenting a novel algorithm for link prediction in the bipartite networks based on potential energy and mutual information. The proposed algorithm majorly works on the four concepts projection, potential energy mutual information and PMIS. Figure 2 briefly describes the process for PMIS score calculation.
In this work we define Potential Energy in the context of graphs. We assumed that a pair of nodes act as an object and that the product of the degree of nodes of the graph can represent mass. The gravitational acceleration g can be represented by the sum of the clustering coefficient of common neighbor between two nodes. For a given pair of nodes, sum of the clustering coefficient of common neighbor will be constant, but for a different pair of nodes the value of the sum of the clustering coefficient of common neighbors will be different. Also, we replaced distance (h) by inverse of the shortest distance (sd) between the pair of nodes.
Potential energy PE(A, B) represents the potential energy between nodes A and B. In the context of a graph, this can be defined as the product of three terms; the product of the degree of nodes, the sum of the clustering coefficient of common neighbors and the shortest distance between nodes.
where \(d_A\), \(d_B\) are degree of nodes A and B. \(cl_z\) is clustering coefficient of z \(\forall z\ne \phi \). sd(A, B) represents the shortest distance between node A and B. When \(z = \phi \); that is no common neighbor between node pair(A, B), and is such cases the value of \(cl_z\) will be .1(constant).
To simplify the understanding of PE, we illustrate it with an example. In Fig. 3, the values of PE between node pairs(B, E), (C, E) and (E, G) are .25, .166, and .04 respectively. We used Eq. (16) for PE calculation and Eq. (8) for calculation of the clustering coefficient of the node. Here PE of node pair(B, E) is greater than (C, E). So it shows that Node pair(B, E) is more likely to be connected than node pair(C, E). PE can distinguish node pairs even if they have no common neighbour. For example, PE of node pair(E, F) and (E, G) are .1 and .04 respectively. So this shows that node pair(E, F) is more likely to be connected than node pair(E, G). The value of the clustering coefficient, product of the degree of nodes and the shortest distance of illustrative examples are given in Table 1.
The potential energy of the network incorporates three different properties in itself. One part of potential energy talks about the product of the degree of vertices. Higher the values of product, higher the PE. If we think from the social networks point of view, vertex having a higher degree has always high likelihood to connect with another vertex. For example, a new person joining Twitter has a higher likelihood to follow a celebrity than not so popular people. Because the degree of celebrities is usually higher. The second part of PE is the clustering coefficient. The significance of the clustering coefficient in social networks is that a person tends to have friends who are also friends with each other. It is very closely related to triadic closure. Triadic closure plays a very important role in link prediction. So the clustering coefficient has the inherent properties of the link prediction. And the third part of PE is the shortest distance between vertices. This is another important feature of networks. Kleinberg found that most of the nodes in the social networks are connected with a very short distance^{32}. Distance between nodes has inverse effect on link prediction. This is also related to the small world phenomenon. In real social networks if the distance between two people is smaller that means they have a higher chance of becoming friends in the future. Illustrative example and results also show the effectiveness of PE.
Algorithm framework
The proposed algorithm initially takes a bipartite graph G(U, V, E) as input and using weighted projection transforms it into a unipartite graph \(G_u(U,E_u)\). In the unipartite graph \(G_u(U, E_u)\), the proposed algorithm computes PE and MI for each node pair using Eqs. (15) and (16) respectively. Then algorithm calculates weight of edge(pattern) of unipartite graph using Eq. (17). After the calculation of the weight of the pattern, the proposed algorithm uses Eq. (18) and compute the PMIS score for each node pair of bipartite graph G(U, V, E).
The proposed algorithm takes weighted projection instead of simple projection. The main reason of taking weighted projection is that in simple projection we lose the topological information of the original bipartite graph. To keep such information we use weighted projection. For example, two bipartite graphs in Fig. 4a,b. are different but their simple projection is the same. Figure 4c shows the simple projection of bipartite graph Fig. 4a,b. But if we take weighted projection then we get two different projected graph Fig. 4d,e. In Fig. 4d edge weight is 1 because only one common neighbor is present in the bipartite graph. But in Fig. 4e edge weight is 2 because there are 2 common neighbors in the bipartite graph. WP(A, B) represents the weight of the edge of the projected graph and we calculate it as follows.
Since we have a projected graph of the original bipartite graph, so by using Eq. (15) we can estimate the \(S^{MI}_{AB}\) for each node pair of the projected graph. Or we can say in another way that for every pattern we have \(S^{MI}_{AB}\).
The weight of the edge in the projected graph plays a very important role in link prediction. This weight is nothing but the weight of pattern, which we have already defined. So now we can say that for each pattern we have three values PE(A, B), WP(A, B) and \(S^{MI}_{AB}\). These values declare the importance of pattern.
Definition 0.1
Total Weight of pattern, \(W_t(A, B)\) represents the sum of all types of weight of the pattern (edge A, B in projected graph) and is defined as follows:
Definition 0.2
Potential EnergyMutual Information based similarity metrics(PMIS): PMIS is the sum of the weight of all pattern covered by a node pair.
If (A, i) is a node pair in bipartite graph \(G=(U, V, E)\) then PMIS can also be defined as
Here W(A, B) is the weight of the pattern \(\{A, B\}\). And \(\Gamma (A, i)\) is the set of patterns covered by node pair (A, i). We can find the value of \(\Gamma (A, i)\) as follows:
So it is clear that every node pair in the bipartite graph cover some pattern. And each pattern has some weight. Therefore PMIS computes the sum of the weight of all pattern covered by a node pair. Here we are using PMIS value as the final link prediction score. Higher PMIS value indicates higher likelihood of the existence of an edge.
Complexity analysis The time complexity of proposed algorithm for a given bipartite graph, G(U, V, E) with two vertex sets, U and V is \(O(U^3)\), where \(V<U\). The algorithm has three steps. Bipartite graph is converted into unipartite graph in the first step with time complexity \(O(Vd^2)\). Here, d denotes the maximum vertex degree of given bipartite graph. Since d is constant, the complexity of first step can be rewritten as O(V). The second step of the algorithm computes PE and MI and the complexity of this step is \(O(V^3)\). PMIS is computed in the last step with \(O(Vd^3)\) time complexity. Therefore, total time complexity of the algorithm is \(O(V^3)\). Moreover, the complexity improves to O(V) when algorithm focuses on candidate node pair instead of focusing on each node pair. This complexity is better than \(O((U+V)^3)\) time complexity of Katz algorithm.
Results and discussion
All experiments were conducted on a Linux server with an Intel XeonE52630 v3 2.40 GHz CPU and 64GB memory running CentOS 7.41708. We implemented PMIL and all other algorithms in Python 3.7.0. In all experiments majorly used networkx, pandas, sklearn, numpy, and matplotlib library.
Evaluation metrics
Two standard metrics are generally used to quantify the accuracy of any prediction algorithms, one is area under the receiver operating characteristic curve (AUC)^{33} and another is Precision^{34}. We have performed an extensive experiment and used four metrics to test the performance of the proposed algorithm. Following are the names of four metrics used for the performance evaluation.
AUC AUC value can be defined as the probability that a randomly chosen missing link (i.e., a link in \(E^p\) ) is given a higher score than a randomly chosen nonexistent link (i.e., a link in \(UE\)). Here if suppose among \(N^t\) independent comparisons, \( N^h, \) times the existing edge having a higher score and \( N^s \) times they have the same score,then AUC score can be calculated by following equation.
In general, a larger AUC value indicates high performance. The AUC value of the ideal result is 1.0.
Precision Precision is defined as the ratio of relevant items selected to the number of items selected. After sorting the scores, if there are \(L_r\) links belonging to the test set among topL candidate links, then Precision is obtained by the following equation.
PredictionPower (PP) This metric is used to check the deviation from the mean randompredictor performance^{20}. PP is computed as follows:
where \(Precision_{Random}\) is the result of randompredictor. And we can compute it by \(\frac{L}{UV(EL)}\).
Precision@K It is the fraction of correct predictions in top k predictions^{34}. In our paper, we computed Precision@10, Precision@20 and Precision@50. Precision@10 means precision at the top 10 position in the ranking result. The higher the values of the metrics are, the better the algorithm is.
To evaluate the performance of our model, we used the Kfold CrossValidation. KFold CV is a technique in which a given data set is split into a K number of sections/folds. Each time one subset is chosen as a probe set and the rest \(K1\) used as training set. Here we have taken the value of K is 10.
Datasets
Datasets We used ten realworld datasets to test the performance of the proposed algorithm. These ten datasets are the following: (1) MovieLens (ML)^{35} dataset contains 100,000 ratings from 943 users on 1682 movies. (2) Enzyme (EN)^{36} is a biological network of drugs and enzymes proteins. It contain 445 drugs nodes, 664 proteins nodes and 2926 drug–target interactions. (3) Southern Women network dataset (SWN)^{37} represents 18 women who participated in 14 social events. (4) Corporate Leadership bipartite graph dataset (CL)^{38} contains the person’s name and company name. (5) Club membership dataset (CM)^{39} contains participation data of corporate officials in social associations. (6) Ionchannels (IC)^{36} is biological network of drug and ionchanel proteins. (7) Countryorganization (C2O) is global network of country and various organization (8) Drug target (Drug)^{40} is a chemical network of drug target interaction. (9) Gprotein coupled receptors (GPC)^{36} is biological network. (10) Malaria(mal)^{41} is a genetic network. Table 2 shows the topological features of all the datasets.
Results To test the strength of the PMIL algorithm, we performed extensive experiments on ten different realworld datasets and compared it with eleven baseline link prediction techniques. Since our proposed algorithm comes under the similaritybased technique, so for the comparison purpose we considered mainly similaritybased algorithm. The baseline link prediction techniques include Common Neighbors (CN), Jaccard Coefficient (JC), Preferential attachment (PA), Cannistraci–Alanis–Ravasi (CAR), Cannistraci–Jaccard (CJC), Cannistraci–Adamic–Adar (CAA), Cannistraci resource allocation (CRA), Nonnegative Matrix Factorization (NMF), Cosine (CS), Potential Link prediction (PLP) and Bipartite projection via Randomwalk (BPR)^{10,21,42}. Out of these 11 baseline algorithms, three works on node neighbourhood mechanism, four works on LCP mechanism, three works on projection mechanism and one works on latent feature mechanism.
The AUC and Precision values of the proposed algorithm and other baseline algorithms are listed in Tables 3 and 4 respectively. In these tables, each row represents the method used in the experiment and each column represents the datasets. The largest value in each column is represented in bold text. In each of these 10 datasets, the test set contains 10% edges and training set contains 90% edges. Table 3 shows the proposed PMIL algorithm outperforms the ten baseline link prediction algorithms on seven datasets for AUC values. But on the CL and IC dataset winner is CAA and PLP respectively. Interestingly, the value of AUC for CAA and CRA are same on C2O dataset. Since the AUC value of the PMIL algorithm is better so if we draw ROC curve by plotting truepositive rates (TPR) versus falsepositive rates (FPR) for varying L values then the total area under the ROCcurve (AUC) will be more. Thus it indicates the better prediction result quality, where L is the list of top links as predicted links. The results in Table 4 demonstrate that PMIL algorithm gives best precision values on six datasets (ML, EN, SWN, CM, Mal, GPC). However, CAA and PLP are winners for CL and IC datasets respectively; and performance of BPR algorithm is best on both C2O and Drug dataset based on the Precision value. Figure 5 shows the effects of the size of the training set on AUC for Drug dataset. We experimented by changing the size of the training set from 40% to 90%. It can be observed from Fig. 5 that on increasing the size of the training set to 0.9, all baseline algorithms as well as PMIL gives better AUC scores.
Figures 6, 7 and 8 describes the precision@K for different values of K. To keep the image clear, we presented the result of twelve algorithms on four datasets. In Fig. 6 we presented the value of precision@10 for the Mal, IC, GPC and Drug datasets. For the IC dataset, CJC and PLP have better performance than PMIL. But our proposed algorithm shows a better result than all eleven baseline algorithms on the Mal and GPC dataset. This improvement is very useful in the recommender system. Especially in Ecommerce, where we are interested to show only the top 10 or top 20 or top K results among the best results to the customer.
Figures 7 and 8 show the precision@20 and precision@50 values respectively. The proposed algorithm shows the highest value of Precision@20 for Mal and GPC dataset. The proposed algorithm gives the best result for precision@50 on Mal, IC, GPC and Drug dataset. Interestingly, all the algorithms based on LCPtheory (CAR, CJC, CAA, CRA) and our proposed algorithm show almost similar results for drug dataset.
Figure 9 describes the predictionpower of all twelve algorithms over the four datasets. For the IC and GPC datasets, PMIL secured the first position. And for Mal datasets, JC and CRA both are at first position, whereas in case of the Drug dataset winner is BPR. PP metric is very useful when we are interested to find which algorithms have a minimum or maximum deviation from the mean randompredictor. So basically it characterises the deviation of the algorithm from the randomness.
Conclusion
In this paper, we introduced a novel approach for link prediction in bipartite networks which is based on the concepts of potential energy and mutual information. The performance of the proposed algorithm is evaluated on ten datasets under different classes and compared with eleven baseline predictors on the basis of AUC, Precision, PredictionPower and Precision@K. We used AUC for evaluation of PMIL, in which we assume that all the links in the networks are independent of each other. However in the realworld networks, links may or may not be independent of each other. PMIL algorithm showed best performance on seven out of the ten datasets used based on AUC score and was reasonably close for the three remaining datasets. There is a possibility that the structural properties of nonsocial network graphs effect the performance; thus it be interesting to study the structural properties, such as the community structure, of different datasets and compare any significant differences between social networks and other datasets. In social networks a connection implies that there was an interaction in the past, but the complete information about the connections may be missing. The other datasets may be noisy and an example could be biological datasets, where the connections are identified with the help of experiments that are not always accurate. Thus, expanding PMIL such that it is capable of link prediction in weighted bipartite graphs is another future direction of research, where weight could be the probability of the occurrence of edge. There are many applications of link prediction, such as recommender systems, community detection, finding hidden relationships, etc. Thus, it would be interesting to explore the future directions and design more efficient algorithms for link prediction in bipartite networks.
References
 1.
LibenNowell, D. & Kleinberg, J. The linkprediction problem for social networks. J. Am. Soc. Inf. Sci. Technol. 58, 1019–1031 (2007).
 2.
Yu, H. et al. Highquality binary protein interaction map of the yeast interactome network. Science 322, 104–110 (2008).
 3.
Mutinda, F. W., Nakashima, A., Takeuchi, K., Sasaki, Y. & Onizuka, M. Time series link prediction using nmf. In 2019 IEEE International Conference on Big Data and Smart Computing (BigComp) 1–8 (IEEE, 2019).
 4.
Li, X. & Chen, H. Recommendation as link prediction in bipartite graphs: a graph kernelbased machine learning approach. Decis. Supp. Syst. 54, 880–890 (2013).
 5.
Kurt, Z., Ozkan, K., Bilge, A. & Gerek, O. N. A similarityinclusive link prediction based recommender system approach. Elektron. IR Elektrotech. 25, 62–69 (2019).
 6.
Lü, L. & Zhou, T. Link prediction in complex networks: a survey. Physica A Stat. Mech. Appl. 390, 1150–1170 (2011).
 7.
Bilgic, M., Namata, G. M. & Getoor, L. Combining collective classification and link prediction. In Seventh IEEE International Conference on Data Mining Workshops (ICDMW 2007) 381–386 (IEEE, 2007).
 8.
Doppa, J. R., Yu, J., Tadepalli, P. & Getoor, L. Chanceconstrained programs for link prediction. In NIPS Workshop on Analyzing Networks and learning with graphs (2009).
 9.
Wang, C., Satuluri, V. & Parthasarathy, S. Local probabilistic models for link prediction. In Seventh IEEE International Conference on Data Mining (ICDM 2007) 322–331 (IEEE, 2007).
 10.
Newman, M. E. Clustering and preferential attachment in growing networks. Phys. Rev. E 64, 025102 (2001).
 11.
Al Hasan, M., Chaoji, V., Salem, S. & Zaki, M. Link prediction using supervised learning. In SDM06: Workshop on Link Analysis, Counterterrorism and Security (2006).
 12.
Sa, H. R. & Prudencio, R. B. Supervised learning for link prediction in weighted networks. In III International Workshop on Web and Text Intelligence (2010).
 13.
Benchettara, N., Kanawati, R. & Rouveirol, C. Supervised machine learning applied to link prediction in bipartite social networks. In 2010 International Conference on Advances in Social Networks Analysis and Mining 326–330 (IEEE, 2010).
 14.
Chang, Y.J. & Kao, H.Y. Link prediction in a bipartite network using Wikipedia revision information. In 2012 Conference on Technologies and Applications of Artificial Intelligence 50–55 (IEEE, 2012).
 15.
Kunegis, J., De Luca, E. W. & Albayrak, S. The link prediction problem in bipartite networks. In International Conference on Information Processing and Management of Uncertainty in Knowledgebased Systems 380–389 (Springer, 2010).
 16.
Xia, S., Dai, B., Lim, E.P., Zhang, Y. & Xing, C. Link prediction for bipartite social networks: the role of structural holes. In Proceedings of the 2012 International Conference on Advances in Social Networks Analysis and Mining (ASONAM 2012) 153–157 (IEEE Computer Society, 2012).
 17.
Shams, B. & Haratizadeh, S. Sibrank: Signed bipartite network analysis for neighborbased collaborative ranking. Physica A Stat. Mech. Appl. 458, 364–377 (2016).
 18.
Allali, O., Magnien, C. & Latapy, M. Link prediction in bipartite graphs using internal links and weighted projection. In 2011 IEEE Conference on Computer Communications Workshops (INFOCOM WKSHPS) 936–941 (IEEE, 2011).
 19.
Li, R.H., Yu, J. X. & Liu, J. Link prediction: the power of maximal entropy random walk. In Proceedings of the 20th ACM International Conference on Information and Knowledge Management 1147–1156 (ACM, 2011).
 20.
Cannistraci, C. V., AlanisLobato, G. & Ravasi, T. From linkprediction in brain connectomes and protein interactomes to the localcommunityparadigm in complex networks. Sci. Rep. 3, 1613 (2013).
 21.
Daminelli, S., Thomas, J. M., Durán, C. & Cannistraci, C. V. Common neighbours and the localcommunityparadigm for topological link prediction in bipartite networks. N. J. Phys. 17, 113037 (2015).
 22.
Baltakiene, M. et al. Maximum Entropy Approach to Link Prediction in Bipartite Networks. arXiv preprint arXiv:1805.04307 (2018).
 23.
Saracco, F., Di Clemente, R., Gabrielli, A. & Squartini, T. Randomizing bipartite networks: the case of the world trade web. Sci. Rep. 5, 10595 (2015).
 24.
Gao, M. et al. Projectionbased link prediction in a bipartite network. Inf. Sci. 376, 158–171 (2017).
 25.
Shakibian, H. & Charkari, N. M. Mutual information model for link prediction in heterogeneous complex networks. Sci. Rep. 7, 44981 (2017).
 26.
Aslan, S., Kaya, B. & Kaya, M. Predicting potential links by using strengthened projections in evolving bipartite networks. Physica A Stat. Mech. Appl. 525, 998–1011 (2019).
 27.
Cimini, G. et al. The statistical physics of realworld networks. Nat. Rev. Phys. 1, 58–71 (2019).
 28.
Boguna, M. et al. Network Geometry. arXiv preprint arXiv:2001.03241 (2020).
 29.
Shannon, C. E. A mathematical theory of communication ACM sigmobile mob. Comput. Commun. Rev. 5, 3–55 (2001).
 30.
Cover, T. M. & Thomas, J. A. Elements of Information Theory (Wiley, New York, 2012).
 31.
Tan, F., Xia, Y. & Zhu, B. Link prediction in complex networks: a mutual information perspective. PLoS ONE 9, e107056 (2014).
 32.
Kleinberg, J. M. Navigation in a small world. Nature 406, 845–845 (2000).
 33.
Hanley, J. A. & McNeil, B. J. The meaning and use of the area under a receiver operating characteristic (roc) curve. Radiology 143, 29–36 (1982).
 34.
Herlocker, J. L., Konstan, J. A., Terveen, L. G. & Riedl, J. T. Evaluating collaborative filtering recommender systems. ACM Trans. Inf. Syst. (TOIS) 22, 5–53 (2004).
 35.
Harper FM, Konstan JA (2016) The movielens datasets: history and context. ACM Trans. Interact. Intell. Syst. (tiis) 5 19.
 36.
Yamanishi, Y., Araki, M., Gutteridge, A., Honda, W. & Kanehisa, M. Prediction of drugtarget interaction networks from the integration of chemical and genomic spaces. Bioinformatics 24, i232–i240 (2008).
 37.
Davis, A., Gardner, B. B. & Gardner, M. R. Deep South: A Social Anthropological Study of Caste and Class (The University of Chicago Press, Chicago, 1941).
 38.
Barnes, R. & Burkett, T. Structural redundancy and multiplicity in corporate networks. Int. Netw. Soc. Netw. Anal. 30, 420 (2010).
 39.
Faust, K. Centrality in affiliation networks. Soc. Netw. 19, 157–191 (1997).
 40.
Yamanishi, Y. et al. Dinies: drugtarget interaction network inference engine based on supervised analysis. Nucl. Acids Res. 42 W39–W45 (2014).
 41.
Rask, T. S., Hansen, D. A., Theander, T. G., Pedersen, A. G. & Lavstsen, T. Plasmodium falciparum erythrocyte membrane protein 1 diversity in seven genomesdivide and conquer. PLoS Comput. Biol. 6, e1000933 (2010).
 42.
Lee, D. D. & Seung, H. S. Algorithms for nonnegative matrix factorization. Adv. Neural Inf. Process. Syst. 13 556–562 (2001).
Acknowledgements
The authors would like to thank Prof. Karmeshu for insightful discussions and constant support.
Author information
Affiliations
Contributions
D.S. intial idea and conception(s), P.K. conducted the experiment(s), P.K. and D.S. analysed the results.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kumar, P., Sharma, D. A potential energy and mutual information based link prediction approach for bipartite networks. Sci Rep 10, 20659 (2020). https://doi.org/10.1038/s41598020773649
Received:
Accepted:
Published:
Further reading

Multimorbidity prediction using link prediction
Scientific Reports (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.