Assessing Methods for Evaluating the Number of Components in Non-Negative Matrix Factorization

Non-negative matrix factorization is a relatively new method of matrix decomposition which factors an <i>m</i> × <i>n</i> data matrix <i>X</i> into an <i>m</i> × <i>k</i> matrix <i>W</i> and a <i>k</i> × <i>n&l...

Descripción completa

Guardado en:
Detalles Bibliográficos
Autores principales: José M. Maisog, Andrew T. DeMarco, Karthik Devarajan, Stanley Young, Paul Fogel, George Luta
Formato: article
Lenguaje:EN
Publicado: MDPI AG 2021
Materias:
PCA
Acceso en línea:https://doaj.org/article/17df6d63695d43868d77feb0c575755d
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
Descripción
Sumario:Non-negative matrix factorization is a relatively new method of matrix decomposition which factors an <i>m</i> × <i>n</i> data matrix <i>X</i> into an <i>m</i> × <i>k</i> matrix <i>W</i> and a <i>k</i> × <i>n</i> matrix <i>H</i>, so that <i>X</i> ≈ <i>W</i> × <i>H</i>. Importantly, all values in <i>X</i>, <i>W</i>, and <i>H</i> are constrained to be non-negative. NMF can be used for dimensionality reduction, since the <i>k</i> columns of <i>W</i> can be considered components into which <i>X</i> has been decomposed. The question arises: how does one choose <i>k</i>? In this paper, we first assess methods for estimating <i>k</i> in the context of NMF in synthetic data. Second, we examine the effect of normalization on this estimate’s accuracy in empirical data. In synthetic data with orthogonal underlying components, methods based on PCA and Brunet’s Cophenetic Correlation Coefficient achieved the highest accuracy. When evaluated on a well-known real dataset, normalization had an unpredictable effect on the estimate. For any given normalization method, the methods for estimating <i>k</i> gave widely varying results. We conclude that when estimating <i>k</i>, it is best not to apply normalization. If the underlying components are known to be orthogonal, then Velicer’s MAP or Minka’s Laplace-PCA method might be best. However, when the orthogonality of the underlying components is unknown, none of the methods seemed preferable.