2016
DOI: 10.1109/tsp.2016.2591510
|View full text |Cite
|
Sign up to set email alerts
|

Efficient and Non-Convex Coordinate Descent for Symmetric Nonnegative Matrix Factorization

Abstract: Given a symmetric nonnegative matrix A, symmetric nonnegative matrix factorization (sym-NMF) is the problem of finding a nonnegative matrix H, usually with much fewer columns than A, such that A ≈ HH T . SymNMF can be used for data analysis and in particular for various clustering tasks. In this paper, we propose simple and very efficient coordinate descent schemes to solve this problem, and that can handle large and sparse input matrices. The effectiveness of our methods is illustrated on synthetic and real-w… Show more

Help me understand this report
View preprint versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
1
1
1

Citation Types

0
61
0

Year Published

2016
2016
2025
2025

Publication Types

Select...
6
2

Relationship

1
7

Authors

Journals

citations
Cited by 43 publications
(61 citation statements)
references
References 31 publications
(56 reference statements)
0
61
0
Order By: Relevance
“…where e 0 is the initial error (see Section 4.1), and e(t) is the error X −X F achieved by an algorithm for a given initialization within t seconds. Since our algorithms are nonincreasing, we have E(t) ∈ [0, 1] for all t, with m n k slack matrix of the 12-gon 12 12 5 slack matrix of the 16-gon 16 16 5 slack matrix of the 20-gon 20 20 6 slack matrix of the 24-gon 24 24 6 slack matrix of the 28-gon 28 28 E(0) = 1 and E(t) → t→∞ 0 if the corresponding algorithm converges towards an exact factorization. In order to illustrate the efficiency of a given algorithm, (11) has the advantage that it makes sense to take the average of E(t) for several initializations and data sets and display a single curve.…”
Section: Comparisons For Different Values Of the Parametersmentioning
confidence: 99%
“…where e 0 is the initial error (see Section 4.1), and e(t) is the error X −X F achieved by an algorithm for a given initialization within t seconds. Since our algorithms are nonincreasing, we have E(t) ∈ [0, 1] for all t, with m n k slack matrix of the 12-gon 12 12 5 slack matrix of the 16-gon 16 16 5 slack matrix of the 20-gon 20 20 6 slack matrix of the 24-gon 24 24 6 slack matrix of the 28-gon 28 28 E(0) = 1 and E(t) → t→∞ 0 if the corresponding algorithm converges towards an exact factorization. In order to illustrate the efficiency of a given algorithm, (11) has the advantage that it makes sense to take the average of E(t) for several initializations and data sets and display a single curve.…”
Section: Comparisons For Different Values Of the Parametersmentioning
confidence: 99%
“…However, CD may not converge to the set of KKT points of SymNMF. Instead, there is an additional NS-SymNMF PGD [22] PNewton [22] ANLS [11] SNMF [10] CD [23] (a) N = 500, K = 60.…”
Section: Numerical Resultsmentioning
confidence: 99%
“…condition given in [23] for checking whether the generated sequence converges to a unique limit point. A heuristic method for checking the condition is additionally provided, which requires, e.g., plotting the norm between the different iterates.…”
Section: Numerical Resultsmentioning
confidence: 99%
“…We remark that the previous work [23] has made the observation that solving SymNMF with the additional constraints X i 2 ≤ 2 Z F , ∀i will not result in any loss of the global optimality. Lemma 2 provides a stronger result, that all KKT points of SymNMF are preserved within a smaller…”
Section: The Proposed Algorithmmentioning
confidence: 96%