Statistics for A flexible R package for nonnegative matrix factorization