“…The use of specialized sample survey software necessary for analysis that properly incorporates the survey structure can often be delayed until initial models (which for linear models may then for simplicity of implementation include stepwise fitting of regressors) have been checked and the most important predictor variables and their interactions have been found, at least provisionally. The use of contextual variables as fixed effects that reduce variation in random effects may have technical ramifications (Haslett et al [17] ), but in the context of SAE modelling, where it is the predictions rather than the parameter estimates that are important, these issues are minor. Similarly, because they are usually swamped by aggregation or averaging, outliers in unit-level data (where they apply at the individual or household level) are much less of an issue than outliers in area-level models.…”