Signed graphs encode positive (attractive) and negative (repulsive) relations between nodes. We extend spectral clustering to signed graphs via the one-parameter family of Signed Power Mean Laplacians, defined as the matrix power mean of normalized standard and signless Laplacians of positive and negative edges. We provide a thorough analysis of the proposed approach in the setting of a general Stochastic Block Model that includes models such as the Labeled Stochastic Block Model and the Censored Block Model. We show that in expectation the signed power mean Laplacian captures the ground truth clusters under reasonable settings where state-of-the-art approaches fail. Moreover, we prove that the eigenvalues and eigenvector of the signed power mean Laplacian concentrate around their expectation under reasonable conditions in the general Stochastic Block Model. Extensive experiments on random graphs and real world datasets confirm the theoretically predicted behaviour of the signed power mean Laplacian and show that it compares favourably with state-of-the-art methods.
The scalar power mean of two non-negative scalars @math is a one-parameter family of means defined for @math as @math . Particular cases are the arithmetic, geometric and harmonic means, as shown in Table . Moreover, the scalar power mean is monotone in the parameter @math , i.e. @math when @math (see @cite_55 , Ch. 3, Thm. 1), which yields the well known arithmetic-geometric-harmonic mean inequality @math . As matrices do not commute, several matrix extensions of the scalar power mean have been introduced, which typically agree if the matrices commute, see e.g. Chapter 4 in @cite_14 . We consider the following matrix extension of the scalar power mean: where @math is the unique positive definite solution of the matrix equation @math . Please note that this definition can be extended to positive semidefinite matrices @cite_21 for @math , as @math exists, whereas for @math a diagonal shift is necessary to ensure that the matrices @math are positive definite. 5pt
@cite_55 @cite_14 @cite_21
2007036852 1520400698 2082401067 2089135543
Let @math be a positive integer and @math be a positive even integer. Let @math be an @math order @math -dimensional real weakly symmetric tensor and @math be a real weakly symmetric positive definite tensor of the same size. @math is called a @math -eigenvalue of @math if @math for some @math . In this paper, we introduce two unconstrained optimization problems and obtain some variational characterizations for the minimum and maximum @math --eigenvalues of @math . Our results extend Auchmuty's unconstrained variational principles for eigenvalues of real symmetric matrices. This unconstrained optimization approach can be used to find a Z-, H-, or D-eigenvalue of an even order weakly symmetric tensor. We provide some numerical results to illustrate the effectiveness of this approach for finding a Z-eigenvalue and for determining the positive semidefiniteness of an even order symmetric tensor. We discuss schemes for exact and approximate computations of permanents, and compare them with each other. Specifically, we analyze and generalize the Belief Propagation (BP) approach to computing the permanent of a non-negative matrix. Known bounds and conjectures are verified in experiments, and some new theoretical relations, bounds and conjectures are proposed. We introduce a fractional free energy functional parameterized by a scalar parameter @math , where @math corresponds to the BP limit and @math corresponds to the exclusion principle Mean-Field (MF) limit, and show monotonicity and continuity of the functional with @math . We observe that the optimal value of @math , where the @math -parameterized functional is equal to the exact free energy (defined as the minus log of the permanent), lies in the @math range, with the low and high values from the range producing provable low and upper bounds for the permanent. Our experimental analysis suggests that the optimal @math varies for different ensembles considered but it always lies in the @math interval. Besides, for all ensembles considered the behavior of the optimal @math is highly distinctive, thus offering a lot of practical potential for estimating permanents of non-negative matrices via the fractional free energy functional approach. Let B i be deterministic real symmetric m × m matrices, and ξ i be independent random scalars with zero mean and “of order of one” (e.g., @math ). We are interested to know under what conditions “typical norm” of the random matrix @math is of order of 1. An evident necessary condition is @math , which, essentially, translates to @math ; a natural conjecture is that the latter condition is sufficient as well. In the paper, we prove a relaxed version of this conjecture, specifically, that under the above condition the typical norm of S N is @math : @math for all Ω > 0 We outline some applications of this result, primarily in investigating the quality of semidefinite relaxations of a general quadratic optimization problem with orthogonality constraints @math , where F is quadratic in X = (X 1,... ,X k ). We show that when F is convex in every one of X j , a natural semidefinite relaxation of the problem is tight within a factor slowly growing with the size m of the matrices @math . We give efficient algorithms for volume sampling, i.e., for picking @math -subsets of the rows of any given matrix with probabilities proportional to the squared volumes of the simplices defined by them and the origin (or the squared volumes of the parallelepipeds defined by these subsets of rows). In other words, we can efficiently sample @math -subsets of @math with probabilities proportional to the corresponding @math by @math principal minors of any given @math by @math positive semi definite matrix. This solves an open problem from the monograph on spectral algorithms by Kannan and Vempala (see Section @math of KV , also implicit in BDM, DRVW ). Our first algorithm for volume sampling @math -subsets of rows from an @math -by- @math matrix runs in @math arithmetic operations (where @math is the exponent of matrix multiplication) and a second variant of it for @math -approximate volume sampling runs in @math arithmetic operations, which is almost linear in the size of the input (i.e., the number of entries) for small @math . Our efficient volume sampling algorithms imply the following results for low-rank matrix approximation: (1) Given @math , in @math arithmetic operations we can find @math of its rows such that projecting onto their span gives a @math -approximation to the matrix of rank @math closest to @math under the Frobenius norm. This improves the @math -approximation of Boutsidis, Drineas and Mahoney BDM and matches the lower bound shown in DRVW . The method of conditional expectations gives a algorithm with the same complexity. The running time can be improved to @math at the cost of losing an extra @math in the approximation factor. (2) The same rows and projection as in the previous point give a @math -approximation to the matrix of rank @math closest to @math under the spectral norm. In this paper, we show an almost matching lower bound of @math , even for @math .