I have a large dataset that is weighted due to:
- probability proportional to size sampling, &
- disproportionate stratification.
I would like to compare the means for two different domains within the same sample, using the weighted means and variance estimates. I believe I could proceed using the Delta approximation of variance of a function of random variables: $$ {\rm V}(\mu_1 - \mu_2) = {\rm V}(\mu_1) + {\rm V}(\mu_2) - 2{\rm Cov}(\mu_1, \mu_2) $$ but I cannot find any precedent, nor a different solution, to this problem in the literature. Nor could I find any other questions on this forum that address this specific problem.