In this paper we develop a likelihood ratio test for independence between multivariate discrete random vectors whose components may be dependent or may have different marginal distributions. A special case that tests independence between two discrete random variables is also discussed. The model under the null hypothesis is constructed based on the marginal distributions of the random variables, while the model under the alternative is built by adding a common random effect which is not restricted to follow the same distribution as those of the other variables. The test statistic is shown to admit an asymptotic distribution of a weighted sum of chi-square random variables, each with one degree of freedom. In order to avoid the high computational cost of finding the parameters of the asymptotic distribution, a permutation-based procedure is introduced. It is shown in simulations that the permutation-based realization yields empirical level close to the nominal level, while achieving power comparable to or higher than that of existing competitors in the literature. The application of the test is illustrated with a real data set.