The impact of survey design on standard multinomial-based methods for a multiway contingency table is studied, under nested loglinear models. The asymptotic null distribution of the Pearson chi-squared test statistic, $X^2$ (or the likelihood ratio test statistic, $G^2$) is obtained as a weighted sum of independent $\chi^2_1$ random variables, and the weights are then related to the familiar design effects (deffs) used by survey samplers. A simple correction to $X^2$ (or $G^2$) is also obtained which requires the knowledge of only the cell deffs and the deffs for collapsed tables (marginals), whenever the model admits a direct solution of likelihood equations under multinomial sampling. Finally, an example on the relative performance of $X^2$ and some corrected $X^2$ statistics in a three-way table is given, using some data from the Canada Health Survey, 1978-1979.
"On Chi-Squared Tests for Multiway Contingency Tables with Cell Proportions Estimated from Survey Data." Ann. Statist. 12 (1) 46 - 60, March, 1984. https://doi.org/10.1214/aos/1176346391