Bi-level Optimization for hyperparameters in Nonnegative Matrix Factorizations