Nonnegative Matrix Factorization (NMF) with Heteroscedastic Uncertainties and Missing data

Abstract

Dimensionality reduction and matrix factorization techniques are important and useful machine-learning techniques in many fields. Nonnegative matrix factorization (NMF) is particularly useful for spectral analysis and image processing in astronomy. I present the vectorized update rules and an independent proof of their convergence for NMF with heteroscedastic measurements and missing data. I release a Python implementation of the rules and use an optical spectroscopic dataset of extragalactic sources as an example for demonstration. A future paper will present results of applying the technique to image processing of planetary disks.

Vectorized update rules for NMF with heteroscedastic measurements and the proof. The code NonnegMFPy is available at https://github.com/guangtunbenzhu/NonnegMFPy and can be installed through PyPI. Comments are most welcome

References in corpus (1)

• A New View of Classification in Astronomy with the Archetype Technique: An Astronomical Case of the NP-complete Set Cover Problem

Cited by in corpus (9)

• The Spectroscopic Data Processing Pipeline for the Dark Energy Spectroscopic Instrument
• The circumgalactic medium of eBOSS emission line galaxies: signatures of galactic outflows in gas distribution and kinematics
• SCARLET: Source separation in multi-band images by Constrained Matrix Factorization
• Using Data Imputation for Signal Separation in High Contrast Imaging
• A statistical study of the optical spectral variability in gamma-ray blazars
• Sensitivity of the Roman Coronagraph Instrument to Exozodiacal Dust
• Reconstructing redshift distributions with photometric galaxy clustering
• A Deep Search for a Strong Diffuse Interstellar Band in the Circumgalactic Medium
• NMF-based GPU accelerated coronagraphy pipeline