Abstract
Given any nonnegative matrix \(A \in \mathbb{R}^{m \times n}\), it is always possible to express A as the sum of a series of nonnegative rank-one matrices. Among the many possible representations of A, the number of terms that contributes the shortest nonnegative rank-one series representation is called the nonnegative rank of A. Computing the exact nonnegative rank and the corresponding factorization are known to be NP-hard. Even if the nonnegative rank is known a priori, no simple procedure exists presently that is able to perform the nonnegative factorization. Based on the Wedderburn rank reduction formula, this paper proposes a heuristic approach to tackle this difficult problem numerically. Starting with A, the idea is to recurrently extrat, whenever possible, a rank-one nonnegative portion from the previous matrix while keeping the residual nonnegative and lowering its rank by one. With a slight modification for symmetry, the method can equally be applied to another important class of completely positive matrices. No convergence can be guaranteed, but repeated restart might help alleviate the difficulty. Extensive numerical testing seems to suggest that the proposed algorithm might serve as a first-step numerical means for exploring the intriguing problem of nonnegative rank factorization.
Similar content being viewed by others
References
Bárány, I.: Sylvester’s question: the probability that n points are in convex position. Ann. Probab. 27(4), 2020–2034 (1999)
Barioli, F., Berman, A.: The maximal cp-rank of rank k completely positive matrices. Linear Algebra Appl. 363, 17–33 (2003). Special issue on nonnegative matrices, M-matrices and their generalizations (Oberwolfach 2000)
Beasley, L.B., Laffey, T.J.: Real rank versus nonnegative rank. Linear Algebra Appl. 431(12), 2330–2335 (2009)
Berman, A., Plemmons, R.J.: Nonnegative Matrices in the Mathematical Sciences, Classics in Applied Mathematics, vol. 9. Society for Industrial and Applied Mathematics (SIAM), Philadelphia (1994). Revised Reprint of the 1979 Original
Berman, A., Rothblum, U.G.: A note on the computation of the CP-rank. Linear Algebra Appl. 419(1), 1–7 (2006)
Berman, A., Shaked-Monderer, N.: Completely Positive Matrices. World Scientific Publishing Co. Inc., River Edge (2003)
Bertsekas, D.P.: Nonlinear Programming, 2nd edn. Athena Scientific, Belmont (1999)
Bocci, C., Carlini, E., Rapallo, F.: Perturbation of matrices and nonnegative rank with a view toward statistical models. SIAM J. Matrix Anal. Appl. 32(4), 1500–1512 (2011)
Bro, R., de Jong, S.: A fast non-negativity-constrained least squares algorithm. J. Chemom. 11(5), 393–401 (1997)
Campbell, S.L., Poole, G.D.: Computing nonnegative rank factorizations. Linear Algebra Appl. 35, 175–182 (1981)
Chu, M.T., Diele, F., Plemmons, R.J., Ragni, S.: Optimality, computation and interpretation of nonnegative matrix factorizations. Available online at http://www4.ncsu.edu/mtchu/Research/Papers/nnmf.ps (2005)
Chu, M.T., Funderlic, R.E., Golub, G.H.: A rank-one reduction formula and its applications to matrix factorizations. SIAM Rev. 37(4), 512–530 (1995)
Chu, M.T., Funderlic, R.E., Golub, G.H.: Rank modifications of semidefinite matrices associated with a secant update formula. SIAM J. Matrix Anal. Appl. 20, 428–436 (1999, electronic)
Chu, M.T., Lin, M.M.: Low-dimensional polytope approximation and its applications to nonnegative matrix factorization. SIAM J. Sci. Comput. 30(3), 1131–1155 (2008)
Cline, R.E., Funderlic, R.E.: The rank of a difference of matrices and associated generalized inverses. Linear Algebra Appl. 24, 185–215 (1979)
Cohen, J.E., Rothblum, U.G.: Nonnegative ranks, decompositions, and factorizations of nonnegative matrices. Linear Algebra Appl. 190, 149–168 (1993)
Donoho, D., Stodden, V.: When does nonnegative matrix factorization give a correct decomposition into parts? In: Proc. 17th Ann. Conf. Neural Information Processing Systems. NIPS, Stanford University, Stanford, CA (2003)
Elsner, L., Nabben, R., Neumann, M.: Orthogonal bases that lead to symmetric nonnegative matrices. Linear Algebra Appl. 271, 323–343 (1998)
Gillis, N.: Nonnegative matrix factorization: complexity, algorithms and applications. Ph.D. thesis, Université catholique de Louvain (2011)
Gillis, N., Glineur, F.: On the geometric interpretation of the nonnegative rank. Linear Algebra Appl. 437(11), 2685–2712 (2012)
Goemans, M.X.: Smallest compact formulation for the permutahedron (2009, preprint)
Hannah, J., Laffey, T.J.: Nonnegative factorization of completely positive matrices. Linear Algebra Appl. 55, 1–9 (1983)
Hopke, P.K.: Receptor Modeling for Air Quality Management. Elsevier, Amsterdam (1991)
Householder, A.S.: The Theory of Matrices in Numerical Analysis. Dover Publications Inc., New York (1975). Reprint of 1964 edition
Hoyer, P.O.: Nonnegative sparse coding. In: Proc. IEEE Workshop Neural Networks for Signal Processing. Martigny (2002)
Hoyer, P.O.: Non-negative matrix factorization with sparseness constraints. J. Mach. Learn. Res. 5, 1457–1469 (2004)
Hubert, L., Meulman, J., Heiser, W.: Two purposes for matrix factorization: a historical appraisal. SIAM Rev. 42(1), 68–82 (2000, electronic)
Jain, S.K., Tynan, J.: Nonnegative rank factorization of a nonnegative matrix A with \(A\sp \dagger A\geq 0\). Linear Multilinear Algebra 51(1), 83–95 (2003)
Jeter, M.W., Pye, W.C.: A note on nonnegative rank factorizations. Linear Algebra Appl. 38, 171–173 (1981)
Jeter, M.W., Pye, W.C.: Some nonnegative matrices without nonnegative rank factorizations. Ind. Math. 32(1), 37–41 (1982)
Kawamoto, T., Hotta, K., Mishima, T., Fujiki, J., Tanaka, M., Kurita, T.: Estimation of single tones from chord sounds using non-negative matrix factorization. Neural Netw World 3, 429–436 (2000)
Kim, H., Park, H.: Nonnegative matrix factorization based on alternating nonnegativity constrained least squares and active set method. SIAM J. Matrix Anal. Appl. 30(2), 713–730 (2008)
Klain, D.A., Rota, G.C.: Introduction to Geometric Probability. Lezioni Lincee. [Lincei Lectures]. Cambridge University Press, Cambridge (1997)
Lawson, C.L., Hanson, R.J.: Solving Least Squares Problems, Classics in Applied Mathematics, vol. 15. Society for Industrial and Applied Mathematics (SIAM), Philadelphia (1995). Revised Reprint of the 1974 Original
Lee, D.D., Seung, H.S.: Learning the parts of objects by non-negative matrix factorization. Nature 401, 788–791 (1999)
Lee, D.D., Seung, H.S.: Algorithms for non-negative matrix factorization. In: Advances in Neural Information Processing Systems, vol. 13, pp. 556–562 (2001)
Levin, B.: On calculating maximum rank one underapproximations for positive arrays. Tech. rep., Division of Biostatistics, Columbia University, New York (1985)
Lin, M.M., Chu, M.T.: On the nonnegative rank of Euclidean distance matrices. Linear Algebra Appl. 433(3), 681–689 (2010)
Paatero, P.: The multilinear engine—a table-driven, least squares program for solving multilinear problems, including the n-way parallel factor analysis model. J. Comput. Graph. Stat. 8(4), 854–888 (1999)
Paatero, P., Tapper, U.: Positive matrix factorization: a non-negative factor model with optimal utilization of error. Environmetrics 5, 111–126 (1994)
Shaked-Monderer, N.: Minimal cp-matrices. ELA 8, 140–157 (2001)
Shaked-Monderer, N.: A note on the cp-rank of matrices generated by a soules matrix. ELA 12, 2–5 (2004)
Sra, S., Dhillon, I.S.: Nonnegative matrix approximation: algorithms and applications. Tech. rep., Department of Computer Sciences, University of Texas at Austin (2006)
Sylvester, J.J.: On a special class of questions on the theory of probabilities. Birmingham British Assoc. Rept., pp. 8–9 (1865)
Thomas, L.B.: Solution to problem 73-14: rank factorization of nonnegative matrices by A. Berman and R. J. Plemmons. SIAM Rev. 16(3), 393–394 (1974)
Vavasis, S.A.: On the complexity of nonnegative matrix factorization. SIAM. J. Optim. 20(3), 1364–1377 (2009)
Wedderburn, J.H.M.: Lectures on Matrices. Dover Publications Inc., New York (1964)
Author information
Authors and Affiliations
Corresponding author
Additional information
Bo Dong’s research was supported in part by the National Natural Science Foundation of China (Grant No. 11101067 and 11171051), TianYuan Special Funds of the National Natural Science Foundation of China (Grant No. 11026164) and the Fundamental Research Funds for the Central Universities.
Matthew M. Lin’s research was supported in part by the National Science Council of Taiwan under grant 101-2115-M-194-007-MY3.
Moody T. Chu’s research was supported in part by the National Science Foundation under grant DMS-1014666.
Rights and permissions
About this article
Cite this article
Dong, B., Lin, M.M. & Chu, M.T. Nonnegative rank factorization—a heuristic approach via rank reduction. Numer Algor 65, 251–274 (2014). https://doi.org/10.1007/s11075-013-9704-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-013-9704-0
Keywords
- Nonnegative matrix
- Nonnegative rank
- Nonnegative matrix factorization
- Nonnegative rank factorization
- Wedderburn rank reduction formula
- Completely positive matrix
- cp-rank