This book is devoted to the study of multivariate discrete q-distributions, which is greatly facilitated by existing multivariate q-sequences and q-functions. Classical multivariate discrete distributions are defined on a sequence of independent and identically distributed Bernoulli trials, with either being a success of a certain rank (level) or a failure. The author relaxes the assumption that the probability of success of a trial is constant by assuming that it varies geometrically with the number of trials and/or the number of successes. The latter is advantageous in the sense that it permits incorporating the experience gained from the previous trials and/or successes, which leads to multivariate discrete q-distributions. Furthermore, q-multinomial and negative q-multinomial formulae are obtained. Next, the book addresses q-multinomial and negative q-multinomial distributions of the first and second kind. The author also examines multiple q-Polya urn model, multivariate q-Polya and inverse q-Polya distributions.