Likelihood estimation of sparse topic distributions in topic models and its applications to Wasserstein document distance calculations
This paper studies the estimation of high-dimensional, discrete, possibly sparse, mixture models in topic models. The data consists of observed multinomial counts of words across independent documents. In topic models, the expected word frequency matrix is assumed to be factorized as a word-topic matrix and a topic-document matrix . Since columns of both matrices represent conditional probabilities belonging to probability simplices, columns of are viewed as -dimensional mixture components that are common to all documents while columns of are viewed as the -dimensional mixture weights that are document specific and are allowed to be sparse. The main interest is to provide sharp, finite sample, -norm convergence rates for estimators of the mixture weights when is either known or unknown. For known , we suggest MLE estimation of . Our non-standard analysis of the MLE not only establishes its convergence rate, but reveals a remarkable property: the MLE, with no extra regularization, can be exactly sparse and contain the true zero pattern of . We further show that the MLE is both minimax optimal and adaptive to the unknown sparsity in a large class of sparse topic distributions. When is unknown, we estimate by optimizing the likelihood function corresponding to a plug in, generic, estimator of . For any estimator that satisfies carefully detailed conditions for proximity to , the resulting estimator of is shown to retain the properties established for the MLE. The ambient dimensions and are allowed to grow with the sample sizes. Our application is to the estimation of 1-Wasserstein distances between document generating distributions. We propose, estimate and analyze new 1-Wasserstein distances between two probabilistic document representations.