If you are familiar with surveyreg, I think it is the right choice.  You have cluster variables, strata variables, a weight variable, and what looks like an interesting model.  The data were collected in a survey, and the inference space is the universe (population) where the survey was applied--not the typical mixed model inference space of all possible populations.  Consequently, the tests and confidence bounds produced by surveyreg are more likely to be correct for that particular population. Nickel says it runs without a memory problem, as well. Steve Denham ... View more