The Mahonian probability distribution on words is asymptotically normal According to the corrigendum, notions and results from the authors’ work are well-known (see, e.g., {[it P. Diaconis}, Group representations in probability and statistics, IMS Lecture Notes-Monograph Series, 11. Hayward, CA: Institute of Mathematical Statistics. vi, 198 p. (1998; Zbl 0695.60012)], p. 128-129).par Summary: The Mahonian statistic is the number of inversions in a permutation of a multiset with $a_i$ elements of type $i, 1 leqslant i leqslant m$. The counting function for this statistic is the $q$ analog of the multinomial coefficient $inom {a_1+cdots +a_m}{a_1,cdots ,a_m}$, and the probability generating function is the normalization of the latter. We give two proofs that the distribution is asymptotically normal. The first is computer-assisted, based on the method of moments.par The Maple package MahonianStat, available from the webpage of this article, can be used by the reader to perform experiments and calculations. Our second proof uses characteristic functions. We then take up the study of a local limit theorem to accompany our central limit theorem. Here our result is less general, and we must be content with a conjecture about further work. Our local limit theorem permits us to conclude that the coefficients of the $q$-multinomial are log-concave, provided one stays near the center (where the largest coefficients reside).