Log-Normal Matrix Completion for Large Scale Link Prediction

01/28/2016 ∙ by Brian Mohtashemi, et al. ∙ 0

The ubiquitous proliferation of online social networks has led to the widescale emergence of relational graphs expressing unique patterns in link formation and descriptive user node features. Matrix Factorization and Completion have become popular methods for Link Prediction due to the low rank nature of mutual node friendship information, and the availability of parallel computer architectures for rapid matrix processing. Current Link Prediction literature has demonstrated vast performance improvement through the utilization of sparsity in addition to the low rank matrix assumption. However, the majority of research has introduced sparsity through the limited L1 or Frobenius norms, instead of considering the more detailed distributions which led to the graph formation and relationship evolution. In particular, social networks have been found to express either Pareto, or more recently discovered, Log Normal distributions. Employing the convexity-inducing Lovasz Extension, we demonstrate how incorporating specific degree distribution information can lead to large scale improvements in Matrix Completion based Link prediction. We introduce Log-Normal Matrix Completion (LNMC), and solve the complex optimization problem by employing Alternating Direction Method of Multipliers. Using data from three popular social networks, our experiments yield up to 5 AUC increase over top-performing non-structured sparsity based methods.



There are no comments yet.


page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

As a result of widespread research on large scale relational data, the matrix completion problem has emerged as a topic of interest in collaborative filtering, link prediction [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]

, and machine learning communities. Relationships between products, people, and organizations, have been found to generate low rank sparse matrices, with a broad distribution of rank and sparsity patterns. More specifically, the node degrees in these networks exhibit well known Probability Mass Functions (PMFs), whose parameters can be determined via Maximum Likelihood Estimation. In collaborative filtering or link prediction applications, row and column degrees may be characterized by differing PMFs, which may be harnessed to provide improved estimation accuracy. Directed networks have unique in-degree and out-degree distributions, whereas undirected networks are symmetric and thus exhibit the same row-wise and column wise degree distributions. Though originally thought to follow strict Power Law Distributions, modern social networks have been found to exhibit Log Normal degree patterns in link formation

   In this work, we propose Log Normal Matrix Completion (LNMC) as an alternative to typical or Frobenius norm constrained matrix completion for Link Prediction. The incorporation of the degree distribution prior generally leads to a non-convex optimization problem. However, by employing the Lovasz extension on the resulting objective, we reduce the problem to a convex minimization over the Lagrangian, which is subsequently solved with Proximal Descent and Alternating Direction Method of Multipliers (ADMM). Through experimentation on Google Plus, Flickr, and Blog Catalog social networks, we demonstrate the advantage of incorporating structured sparsity information in the resulting optimization problem.

Ii Related Work

Link prediction has been thoroughly researched in the field of social network analysis as an essential element in forecasting future relationships, estimating unknown acquaintances, and deriving shared attributes. In particular, [18]

introduces the concept of the Social Attribute Network, and uses it to predict the formation and dissolution of links. Their method combines features from matrix factorization, Adamic Adar, and Random walk with Restart using logistic regression to give link probabilities. However, the calculations of such inputs may be time-intensive, and shared attributes may be unlikely, leading to non-descriptive feature vectors.

   Matrix Completion for Link Prediction has previously been investigated within the Positive Unlabeled (PU) Learning Framework, where the nuclear norm regularizes a weighted value-specific objective function [19]. Although the weighted objective improves the prediction results, the subsequent optimization is non-convex and thus subject to instability. Binary Matrix completion employing proximal gradient descent is studied in [20], however, sparsity is not considered, and Link Prediction is not included in the experiment section. The structural constraints that must be satisfied for provably exact completion are described in [21]. In this technical report, the required cardinality of uniformly selected elements is bounded based on the rank of the matrix. Unique rank bounds for matrix completion are considered in [22], where the Schatten p-Norm is utilized on the singluar values of the matrix. Matrix Completion for Power Law distributed samples is studied in [23], where various models are compared, including the Random Graph, Chung Lu-Vu, Preferential Attachment, and Forest Fire models. However, link prediction is not considered and the resulting optimization problem is non-convex.
   The concept of simultaneously sparse and low rank matrices was introduced in [24]

, where Incremental Proximal Descent is employed to sequentially minimize the objective, and threshold the singular values and matrix entries. Due to the sequentiality of the optimization, the memory footprint is reduced, however, the objective is non-convex and may result in a local minimum solution. Also, the tested methods employed in simulation are elementary, and more advanced techniques are well known in the link prediction community. Simultaneous row and column-wise sparsity is discussed in

[23], where a Laplacian based norm is employed on rows and a Dirichlet semi-norm is utilized on columns. A comparison between nuclear and graph based norms is additionally provided. In [25], Kim et. al present a matrix factorization method which utilizes group wise sparsity, to enable specifically targeted regularization. However, the datasets which we utilize do not identify group membership, and thus we will not consider affiliation in our prediction models.
   Structured sparsity was thoroughly investigated in [26], and applied to Graphical Model Learning. However, the paper focuses solely on the Pareto Distribution which characterizes scale-free networks, and does not cover the Log Normal Methods which are presented in this paper. Also, Link Prediction is not considered in the experimental section. Node specific degree priors are introduced in [27], and the Lovasz Extension is additionally employed to learn scale free networks commonly formed by Gaussian Models. However, the stability of the edge rank updating is not proven, and Log Normally distributed networks are not considered.
   The Lovasz Extension and background theory are presented in [28], where Bach provides an overview on submodular functions and minimization.

Iii Proposed Approach

Iii-a Link Prediction

In this paper, we consider social network graphs, since they have been proven to follow Pareto, and more recently discovered, Log Normal, degree distributions. The Social Network Link Prediction problem involves estimating the link status, , between node and node , where is limited to binary outcomes. Together, the set of all nodes, , and links, , form the graph , where is only partially known. Unknown link statuses may exist when either the relationship between and is non-public, or the observation is considered unreliable over several crawls of the social network. Combined, the observations can be expressed in the form of a partial adjacency matrix, , which contains all known values in the set of observed pairs, . Unmeasured states between two nodes are set to 0 in . This matrix can be stored in sparse format for memory conservation, and operation complexity reduction.

Iii-B Structured Sparsity based Matrix Completion for Link Prediction

As demonstrated in [19, 24, 20], Matrix Completion involves solving for unknown entries in matrices by employing the low-rank assumption in addition to other side information regarding matrix formation and evolution. Traditionally, matrix completion problems are expressed as



is the Frobenius norm, and is the nuclear norm (Schatten p-norm with ). The nuclear norm can be defined as


where is the eigenvalue, when arranged in decreasing order, and and are the row count and column count, respectively. In this paper, is assumed equal to . is the estimated complete matrix after convergence is attained. Generally, these problems are solved using proximal gradient descent, which employs singular value thresholding on each iteration [29]. However, this problem generally lacks incorporation of prior sparsity information encoded into the matrix. Thus we augment the problem as


where is defined as follows:


Here, is a sparsity inducing term, where implies that the sparsity is applied on matrix rows, j implies sparsity is applied on matrix columns, is the prior in-degree distribution, and is the out-degree distribution. For the rest of this paper, we will consider the case of symmetric adjacency matrices, and thus set to .

Iii-C Log-Normal Degree Prior

As demonstrated in [17], many social networks, including Google+, tend to exhibit the Log-Normal Degree Distribution


Thus we derive as the Maximum Likelihood Estimate


where is the degree of the row of , which simplifies to the following:


This is equivalent to a summation of scaled Pareto Distributions with shape parameter 1 added to additional square terms. Thus the final optimization problem becomes


Due to the presence of the log term in the optimization, convex methods cannot be directly applied to the minimization, since the problem is not guaranteed to have an absolute minimum. Optimization of this problem is a multi-part minimization, which can be solved using the Alternating Direction Method of Multipliers (ADMM).

Iii-D Optimization

ADMM allows the optimization problem to be split into less complex sub-problems, which can be solved using convex minimization techniques. In order to decouple (8) into smaller subproblems, the additional variable, , is introduced as

Expressing the problem in ADMM update form, the sequential optimization becomes


In practice, step size values, , in the range have been found to work well. Convergence is assumed, and the sequence is terminated once . The initial values, , and are set to zeros matrices. Although ADMM has slow convergence properties, a relatively accurate solution can be attained in a few iterations. Due to the convexity of the initial equation, proximal gradient descent is employed for minimization. The proximal gradient method minimizes problems of the form


using the gradient and proximal operator as


where , and is a multiplier utilized on each gradient descent round. Typically a value of is sufficient for , leading to rapid convergence in rounds, however, a value would result in slower, but more accurate minimization. The optimal value for is determined through experimentation. For Log-Normal Matrix Completion, , and . The proximal operator of becomes a sequential thresholding on the eigenvalues, , of the argument in (13)



is the matrix of eigenvectors. The subproblem reaches convergence when

. The noise of the matrix is reduced through sequential thresholding, leaving only the strongest components of the low rank matrix. This algorithm is advantageous due to rapid convergence properties and automatic rank selection. Known as the Iterative Soft Thresholding Algorithm (ISTA), this method can be parallelized for gradient calculation and recombined for the Eigenvalue decomposition. Although the interim result of each round of minimization is generally not sparse, matrix entries with values below a given threshhold can be forced to 0 to allow sparse matrix Eigenvalue Decomposition (such as eigs in Matlab) to be performed with minimal error.

Iii-E Lovasz Extension

(10) is a non-convex optimization problem due to the log of the set cardinality function. However, the problem can be altered into a convex form using the Lovasz Extension on the submodular set function. As described in [28], the Lovasz Extension takes on the following form:


Here, is a permutation of which ensures components of are ordered in decreasing fashion, , and is a submodular set function. The Lovasz Extension is always convex when is submodular, thus allowing convex optimization techniques to be used on the resulting transformed problem.
In order to transform each individual row of sampled relationship information into a set, , the support function, is utilized. As a result , where is the number of columns present in the matrix X. A submodular set function must obey the relationship


where , and is an additional set element. In this paper, is a log-normal transformation on the degree . The degree, , is modular, and thus follows (16) with strict equality. Thus for to be sub modular, the subsequent transformation of the degree must be submodular as well.
After applying the Lovasz Extension to (7), the result is


Here, is used in order to maintain the positivity required for the Lovasz Extension to remain convex. Further details regarding the optimization of this problem can be obtained in Appendix A.

Iii-F Considerations

In order for (17) to be utilized, (7) must remain a submodular function of the degree. Thus, both the first derivative and the second derivative of the function must remain positive, creating the following constraint:


is introduced to prevent the left side of the inequality from approaching . In practice, a small constant is also subtracted or added from the obtained set function in order to assure that . These small coefficients are determined during the Cross Validation phase, after obtaining the optimal and values which satisfy the given constraints.

Iv Experiment

In order to compare the performance of the LNMC method with other popular Link Prediction methods, an experiment was performed using several data sets from existing literature:

  1. Google + - The Google + dataset [18] contains nodes and links, captured in AUG 2011. The data contains both Graph topology and node attribute information; however, the side-features are removed since our method requires edge status only.

  2. Flickr - Flickr is a social network based on image hosting, where users form communities and friendships based on common interests. The Flickr dataset [30] contains nodes, links, and groups. Group affiliation was discarded due to irrelevance to the LNMC method.

  3. Blog Catalog - Blog Catalog [30] is a blogging site where users can form friendships, and acquire group membership. The utilized dataset contains nodes, links, and groups. Again, for the context of this paper, the group information was removed.

As seen in Fig. 1

, all datasets follow a roughly Log-Normal distribution, with varying amounts of degree sparsity, and variance. Due to the high number of low degree nodes in the Google+ dataset, all points appear constrained to the left of the plot axis; however, as we will illustrate, the Log-Normal Distribution is still superior to the Pareto Distribution for link prediction. During the training phase,

of the data was removed in order to use for future predictions. For the purposes of demonstration, only of the highest degree nodes are maintained for adjacency matrix formation.

(a) Google+
(b) Flickr
(c) Blog Cat.
Fig. 1:

Empirical Node Degree Data and Fitted Log-Normal Probability Distribution Functions

V Results

V-a Baseline Methods and Performance Metrics

In order to understand the advantage of LNMC, the results are compared against the following methods:

  1. Matrix Completion with Pareto Sparsity (MCPS) - MCPS [26] utilizes the same algorithm which we have outlined in the paper with the exception of the prior. MCPS employs the Pareto Distribution .

  2. Matrix Completion with Sparsity (MCLS) - MCLS is used by Richard et al. [24], and represents one of the first attempts at incorporating sparsity with the Low Rank assumption.

  3. Logistic Regression (MF + RwR + AA) - In their paper on Social Attribute Networks, Gong et al. [17] provide a method which combines features from Matrix Factorization, Random Walks with Restart, and Adamic Adar, which effectively solves the link prediction problem with high accuracy. In this paper, the attributes are removed from the network for equal comparison with our method.

In order to provide a fair basis on which to judge the performance, Area Under the Curve (AUC) is employed for comparison. By utilizing the AUC as the performance metric, we avoid the need for data balancing, a process which frequently results in undersampling negative samples. Thus, all methods can benefit from the additional training data.
The results are obtained via Cross Validation, using a random sampling method for hyper-parameter selection. The rounds are averaged to produce the results shown in Table I.

V-B Performance Comparison

As demonstrated in Fig. 2, LNMC outperforms MCPS, MCLS, and LR, on the Google Plus dataset. Due to the highly Log-Normal characteristic [17] of the data set, LNMC’s fine-tuned degree specific prior captures the degree distribution behavior in combination with the low rank features of the data, leading to high AUC values. The high number of true positives compared to the false positive rate leads to jagged graph distribution. In Fig 3

, it is clear that matrix completion with Pareto Sparsity produces low AUC values due to the inaccurate distribution representation. Similarly the LR method fails to capture accurate low rank information because the low rank matrix factorization is done prior to the the gradient descent training for Logistic Regression.

Fig. 2: Receiver Operating Characteristic for Google Plus Data

Due to the Pareto nature of the Flickr dataset, both the LNMC and MCPS methods perform the same. As can be seen in (17), LNMC can adapt to Scale Free Networks when the first term is small compared to the second term. Logistic Regression performs poorly since the features are set, whereas Matrix Completion methods automatically select the number of latent parameters to utilize.

Fig. 3: Receiver Operating Characteristic for Flickr Data

As seen in Fig. 4, LNMC outperforms the Pareto Sparisty based matrix completion, due to the inclusion of the squared log terms. The sparsity used in the MCLS method is insufficiently descriptive for accurate matrix estimation. Thus Logistic Regression, which incorporates more descriptive features outperforms the MCLS method.

Fig. 4: Receiver Operating Characteristic for Blog Catalog Data

For purposes of comparison, AUC values for each method and dataset, are contained in Table I. As highlighted by the AUC Table, LNMC provides optimal results over all datasets.

Google+ .8541 .8439 .8113 .8434
Flickr .9052 .9052 .8504 .8972
Blog Catalog .7918 .7846 .7150 .7727
TABLE I: AUC Performance Comparison

Vi Conclusion

As demonstrated both theoretically, and experimentally, LNMC is able to sufficiently encapsulate the advantages of Pareto Sparsity in addition to Log Normal Sparsity. Previously described by Gong et al. in [17], many modern social networks with undirected graph topologies exhibit Log Normal degree distributions. Thus by incorporating the degree-specific prior the optimization encourages convergence to a Log-Normal degree distribution. Due to the non-convexity of solving the joint low-rank and structured sparsity inducing prior, the Lovasz Extension is introduced to solve the complex problem efficiently. Through analysis on three datasets, and using 3 top performing methods, we provide results which exceed the current optimum. These results reveal the fundamental value of prior degree information in Link Prediction, and can provide insight into understanding the complex dynamics which cause links to form in a similar way across different networks.
In future research we plan to investigate the incorporation of side information into the objective. Node attributes introduce additional challenges, including missing features, and additional training complexities.


  • [1] H. R. Sa and R. B. Prudencio, “Supervised Learning for Link Prediction in Weighted Networks,” Center of Informatics, Federal University of Pernambuco, Tech. Rep.
  • [2] Y. Sun, R. Barber, M. Gupta, C. C. Aggarwal, and J. Han, “Co-Author Relationship Prediction in Heterogeneous Bibliographic Networks,” University of Illinois at Urbana-Champaign, Tech. Rep.
  • [3] G.-J. Qi, C. C. Aggarwal, and T. Huang, “Link Prediction across Networks by Biased Cross-Network Sampling,” University of Illinois at Urbana-Champaign, Tech. Rep.
  • [4] P. Sarkar, D. Chakrabarti, and M. I. Jordan, “Nonparametric Link Prediction in Dynamic Networks,” University of California Berkeley, Tech. Rep.
  • [5] Z. Lu, B. Savas, W. Tang, and I. Dhillon, “ Supervised Link Prediction Using Multiple Sources,” in IEEE 10th International Conference on Data Mining, 2010, pp. 923–928.
  • [6] J. Zhu, “ Max-Margin Nonparametric Latent Feature Models for Link Prediction,” in Proceedings of the 29th International Conference on Machine Learning, 2012.
  • [7] K. T. Miller, T. L. Griffiths, and M. I. Jordan, “ Nonparametric Latent Feature Models for Link Prediction,” in Proceedings on Neural Information Processing Systems, 2009.
  • [8] L. Lu and T. Zhou, “ Link Prediction in Complex Networks: A Survey,” University of Fribourg, Chermin du Musee, Fribourg, Switzerland, Tech. Rep.
  • [9] J. Leskovec, D. Huttenlocher, and J. Kleinberg, “ Predicting Positive and Negative Links in Online Social Networks,” in International World Wide Web Conference, 2010.
  • [10] Y. Dong, J. Tang, S. Wu, and J. Tian, “ Link Prediction and Recommendation across Heterogeneous Social Networks,” in IEEE 12th International Conference on Data Mining, 2012.
  • [11] P. S. Yu, J. Han, and C. Faloutsos, Link Mining: Models, Algorithms, and Applications.   New York, NY: Springer, 2010.
  • [12] D. Li, Z. Xu, S. Li, and X. Sun, “Link Prediction in Social Networks Based on Hypergraph,” in International World Wide Web Conference, 2013.
  • [13] H.-H. Chen, L. Gou, X. Zhang, and C. L. Giles, “Capturing Missing Edges in Social Networks Using Vertex Similarity,” in K-CAP-11, 2011.
  • [14] E. Perez-Cervantes, J. M. Chalco, M. Oliveira, and R. Cesar, “ Using Link Prediction to Estimate the Collaborative Iinfluence of Researchers,” in IEEE 9th International Conference on eScience, 2013, pp. 293–300.
  • [15] Z. Yin, M. Gupta, T. Weninger, and J. Han, “ A Unified Framework for Link Prediction Using Random Walks,” in International Conference on Advances in Social Networks Analysis and Mining, 2010, pp. 152–159.
  • [16] P. Symeondis, E. Tiakas, and Y. Manolopoulos, “ Transitive Node Similarity for Link Prediction in Social Networks with Positive and Negative Links,” in RecSys2010, 2010.
  • [17] N. Z. Gong, W. Xu, and L. Huang, “Evolution of Social-Attribute Networks: Measurements, Modeling, and Implications using Google+,” in IMC, 2012, pp. 1–14.
  • [18] N. Z. Gong, A. Talwalkar, and L. Mackey, “Joint link prediction and attribute inference using a social-attribute network,” ACM Tranactions on Intelligent Systems and Technology, vol. 5, pp. 1–14, 2014.
  • [19] C.-J. Hsieh, N. Natarajan, and I. Dhillon, “PU Learning for Matrix Completion,” in International Conference on Machine Learning 32, 2015, 2015, pp. 1–10.
  • [20] M. A. Davenport, Y. Plan, E. van den Berg, and M. Wootters, “1-Bit Matrix Completion,” Georgia Institute of Technology, Tech. Rep., 2014.
  • [21] Y. Chen, S. Bhojanapalli, S. Sanghavi, and R. Ward, “Completing Any Low Rank Matrix Provably,” University of California Berkeley, Tech. Rep., 2014.
  • [22] F. Nie, H. Wang, X. Cai, H. Huang, and C. Ding, “Robust Matrix Completion via Joint Schatten p-Norm and lp-Norm Minimization,” in IEEE International Conference on Data Mining, 2012, pp. 1–9.
  • [23] R. Meka, P. Jain, and I. S. Dhillon, “Matrix Completion from Power-Law Distributed Samples,” in Neural Information Processing Systems, 2015, pp. 1–9.
  • [24] E. Richard, P.-A. Savalle, and N. Vayatis, “Estimation of Simultaneously Sparse and Low Rank Matrices,” in Proceedings of the 29th International Conference on Machine Learning, 2012, pp. 1–8.
  • [25] J. Kim, R. Monteiro, and H. Park, “Group Sparsity in Nonnegative Matrix Completion,” in Proceedings of the SIAM International Conference on Data, 2012, pp. 1–12.
  • [26] A. Defazio and T. S. Caetano, “A Convex Formulation for Learning Scale-Free Networks via Submodular Relaxation,” in Advances in Neural Information Processing Systems 25, 2015, pp. 1–9.
  • [27] Q. Tang, S. Sun, C. Yang, and J. Xu, “Learning Scale Free Network by Node Specific Degree Prior,” Toyota Technical Institute, Tech. Rep., 2015.
  • [28] F. Bach, “Learning with Submodular Functions: A Convex Optimization Perspective,” Ecole Normale Superieure, Tech. Rep., 2013.
  • [29] S. Gu, L. Zhang, W. Zuo, and X. Feng, “Weighted Nuclear Norm Minimization with Application to Image Denoising,” in Computer Vision and Pattern Recognition, 2014.
  • [30] L. Tang and H. Liu, “Scalable Learning of Collective Behavior based on Sparse Social Dimensions,” in Proceedings of the 18th ACM Conference on Information and Knowledge Management, 2009, pp. 1–10.