Pseudo Bayesian Estimation of One-way ANOVA Model in Complex Surveys

Abstract

We devise survey-weighted pseudo posterior distribution estimators under 2-stage informative sampling of both primary clusters and secondary nested units for a one-way ANOVA population generating model as a simple canonical case where population model random effects are defined to be coincident with the primary clusters. We consider estimation on an observed informative sample under both an augmented pseudo likelihood that co-samples random effects, as well as an integrated likelihood that marginalizes out the random effects from the survey-weighted augmented pseudo likelihood. This paper includes a theoretical exposition that enumerates easily verified conditions for which estimation under the augmented pseudo posterior is guaranteed to be consistent at the true generating parameters. We reveal in simulation that both approaches produce asymptotically unbiased estimation of the generating hyperparameters for the random effects when a key condition on the sum of within cluster weighted residuals is met. We present a comparison with frequentist EM and a methods that requires pairwise sampling weights.Comment: 46 pages, 9 figure