Model-based small area estimation under informative sampling 5. Concluding remarks

In this paper, we studied model-based small area estimation for different levels of design informativeness under a nested error linear regression model for the population units. Estimators considered were the EBLUP, the pseudo-EBLUP (You and Rao 2002) and an estimator given by Pfeffermann and Sverchkov (2007). The EBLUP and the pseudo-EBLUP were computed under two scenarios: (i) Ignore informative sampling and assume that the population model holds for the sample; (ii) Take account of informative sampling by using a suitable function of the unit selection probability p j | i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGWbWaaS baaSqaamaaeiaabaGaamOAaaGaayjcSdGaamyAaaqabaaaaa@3CF9@ as an additional auxiliary variable in the sample model.

Results from a simulation study showed that design informativeness can have a big impact on the bias and MSE of the EBLUP that ignores informative sampling (scenario (i)). Results under scenario (ii) showed that the EBLUP, based on the augmented model, performs extremely well in terms of bias and MSE, provided that the augmenting variable is chosen properly. The bias-adjusted estimator of Pfeffermann and Sverchkov (2007) also performed well under informative sampling in terms of bias but its MSE is significantly larger than the corresponding MSE of the EBLUP and the pseudo-EBLUP based on the augmented model. Pseudo-EBLUP under scenario (i) performed significantly better than the corresponding EBLUP. It can be significantly improved by using the augmented model, similar to the case of EBLUP.

An advantage of the augmented model approach is that no new theory is required for estimation and MSE estimation. However, the area mean G ¯ i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaaceWGhbGbae badaWgaaWcbaGaamyAaaqabaaaaa@3A63@ of the augmenting variable g i j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadMgacaWGQbaabeaaaaa@3B5A@ is required, unlike in the approach of Pfeffermann and Sverchkov (2007). For some choices of g i j , G ¯ i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadMgacaWGQbaabeaakiaacYcaceWGhbGbaebadaWgaaWc baGaamyAaaqabaaaaa@3E12@ is readily known; for example g i j = p j | i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadMgacaWGQbaabeaakiabg2da9iaadchadaWgaaWcbaWa aqGaaeaacaWGQbaacaGLiWoacaWGPbaabeaaaaa@40FE@ gives G ¯ i = 1 / N i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaaceWGhbGbae badaWgaaWcbaGaamyAaaqabaGccqGH9aqpdaWcgaqaaiaaigdaaeaa caWGobWaaSbaaSqaaiaadMgaaeqaaaaaaaa@3E31@ and g i j = n i w j | i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadMgacaWGQbaabeaakiabg2da9iaad6gadaWgaaWcbaGa amyAaaqabaGccaWG3bWaaSbaaSqaamaaeiaabaGaamOAaaGaayjcSd GaamyAaaqabaaaaa@431C@ gives G ¯ i = n i W ¯ i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaaceWGhbGbae badaWgaaWcbaGaamyAaaqabaGccqGH9aqpcaWGUbWaaSbaaSqaaiaa dMgaaeqaaOGabm4vayaaraWaaSbaaSqaaiaadMgaaeqaaaaa@3F98@ and W ¯ i MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaaceWGxbGbae badaWgaaWcbaGaamyAaaqabaaaaa@3A73@ is often known for some surveys. We have also given a method of choosing the augmenting variable g i j . MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadMgacaWGQbaabeaakiaac6caaaa@3C16@

In this paper, we focused on the special case where all the areas are sampled. Extension of the augmented model approach to handle non-sampled areas requires the knowledge of the area means G ¯ i , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaaceWGhbGbae badaWgaaWcbaGaamyAaaqabaGccaGGSaaaaa@3B1D@ as well as the area selection probabilities, p i , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrpipeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXddrpe0=1qpeea0=yrVue9 Fve9Fve8meaabaqaciaacaGaaeqabaWaaeaaeaaakeaacaWGWbWaaS baaSqaaiaadMgaaeqaaOGaaiilaaaa@3B2E@ for the non-sampled areas. This extension is currently under study.

Acknowledgements

We are thankful to the Associate Editor, the referees and M. Sverchkov for many constructive comments and suggestions.

References

Asparouhov, T. (2006). General multi-level modelling with sampling weights. Communication in Statistics, Theory and Methods, 439-460.

Battese, G.E., Harter, R.M. and Fuller, W.A. (1988). An error component model for prediction of county crop areas using survey and satellite data. Journal of the American Statistical Association, 83, 28-36.

Henderson, C.R. (1975). Best linear unbiased estimation and prediction under a selection model. Biometrics, 31, 423-447.

Pfeffermann, D., and Sverchkov, M. (2003). Small area estimation under informative sampling. Proceedings of the Survey Research Methods Section, American Statistical Association, 3284-3295.

Pfeffermann, D., and Sverchkov, M. (2007). Small-area estimation under informative probability sampling of areas and within the selected areas. Journal of the American Statistical Association, 102, 480, 1427-1439.

Prasad, N.G.N., and Rao, J.N.K. (1999). On robust small area estimation using a simple random effects model. Survey Methodology, 25, 1, 67-72.

Rao, J.N.K. (1965). On two simple schemes of unequal probability sampling without replacement. Journal of the Indian Statistical Association, 3, 173-180.

Rao, J.N.K. (2003). Small Area Estimation, New York: John Wiley & Sons, Inc.

Sampford, M.R. (1967). On sampling without replacement with unequal probabilities of selection. Biometrika, 54, 499-513.

Skinner, C.J. (1994). Sampling models and weights. Proceedings of the Survey Research Methods Section, American Statistical Association, 133-142.

Stefan, M. (2005). Contributions à l’estimation pour petits domaines. Ph.D. thesis, Université Libre de Bruxelles.

Torabi, M., and Rao, J.N.K. (2010). Mean squared error estimators of small area means using survey weights. The Canadian Journal of Statistics, 38, 4, 595-608.

You, Y., and Rao, J.N.K. (2002). A pseudo-empirical best linear unbiased prediction approach to small area estimation using survey weights. The Canadian Journal of Statistics, 30, 3, 431-439.

Date modified: