Identifying Low-Dimensional Structures in Markov Chains: A Nonnegative Matrix Factorization Approach