Nonnegative Matrix Factorization with Zellner Penalty