Characterization of Stability of Non-Negative Matrix Factorization Models: An Application to Single-Cell Data