Flexible and Hierarchical Prior for Bayesian Nonnegative Matrix Factorization