8. Conclusions

Isabel Molina, J.N.K. Rao and Gauri Sankar Datta

Previous

The following major conclusions may be drawn from the results of our simulation study on the estimation of small area means, based on the Fay-Herriot area-level model when the number of areas is modest in size ( say  m = 15 ) : MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaamaabmqaba Gaae4CaiaabggacaqG5bGaaeiiaiaad2gacqGH9aqpcaaIXaGaaGyn aaGaayjkaiaawMcaaiaacQdaaaa@4230@  1) Under the Fay-Herriot model with a value of random effects variance, A , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaadgeaca GGSaaaaa@3A73@  clearly away from zero, the PTE does not seem to noticeably improve efficiency relative to the usual EBLUP unless the significance level is taken small (α0.1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaacIcacq aHXoqycqGHKjYOcaaIWaGaaGOlaiaaigdaaaa@3F2A@  in our simulation study). 2) Our simulation results indicate that using the PT procedure with a moderate α , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiabeg7aHj aacYcaaaa@3B4C@  in particular α = 0.2 , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiabeg7aHj abg2da9iaaicdacaGGUaGaaGOmaiaacYcaaaa@3E7A@  to estimate the MSE of the usual EBLUP leads to a reduction in bias as compared with the usual MSE estimator. Hence, we recommend the use of mse PT ( θ ^ RE , i ) , MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaab2gaca qGZbGaaeyzamaaBaaaleaacaqGqbGaaeivaaqabaGcdaqadeqaaiqb eI7aXzaajaWaaSbaaSqaaiaabkfacaqGfbGaaGilaiaadMgaaeqaaa GccaGLOaGaayzkaaGaaiilaaaa@4522@  given by (4.1), to estimate the MSE of the EBLUP. 3) Among the estimators that attach a strictly positive weight to the direct estimator for all areas, we recommend the combined estimator REML-AML given by (6.2), because it achieves slightly higher efficiency than the EBLUP based on AML and the PT-AML given by (6.1). 4) For estimating the MSE of the recommended REML-AML estimator, the estimator mse PT ( θ ^ REAML , i ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiFu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9LqFf0x e9q8qqvqFr0dXdbrVc=b0P0xb9peuD0xXdbvk9qq=xd9qqaq=Jf9sr 0=vr0=vrWZqaaeaabiGaaiaacaqabeaadaqaaqaaaOqaaiaab2gaca qGZbGaaeyzamaaBaaaleaacaqGqbGaaeivaaqabaGcdaqadeqaaiqb eI7aXzaajaWaaSbaaSqaaiaabkfacaqGfbGaaeyqaiaab2eacaqGmb GaaGilaiaadMgaaeqaaaGccaGLOaGaayzkaaaaaa@46D5@  given by (6.5) performs better than the alternative ones. 5) Our results on prediction intervals, based on normal theory, indicate that the good performance of the proposed MSE estimators may not translate to coverage properties of these intervals. Construction of prediction intervals that lead to accurate coverages, using the proposed MSE estimates, appears to be a difficult task.

Smooth alternatives to the preliminary test estimates in the case of location parameters have been proposed in the literature using weighted means of the estimates obtained under the null and alternative hypotheses, with weights depending on the test statistic, see e.g., Saleh (2006). Mean squared error estimates of this kind have not been studied and we leave this subject for further research.

Acknowledgements

We would like to thank the editor for very constructive suggestions. Gauri S. Datta’s research was partially supported through the grant H98230-11-1-0208 from the National Security Agency, Isabel Molina’s research by grants ref. MTM2009-09473, MTM2012-37077-C02-01 and SEJ2007-64500 from the Spanish Ministerio de Educación y Ciencia and J.N.K. Rao’s research by the Natural Sciences and Engineering Research Council of Canada.

References

Bancroft, T.A. (1944). On biases in estimation due to the use of preliminary tests of significance. The Annals of Mathematical Statistics, 15, 190-204.

Chatterjee, S., Lahiri, P. and Li, H. (2008). Parametric bootstrap approximation to the distribution of EBLUP and related prediction intervals in linear mixed models. The Annals of Statistics, 36, 1221-1245.

Das, K., Jiang, J. and Rao, J.N.K. (2004). Mean squared error of empirical predictor. The Annals of Statistics, 32, 818-840.

Datta, G.S., and Lahiri, P. (2000). A unified measure of uncertainty of estimated best linear unbiased predictors in small area estimation problems. Statistica Sinica, 10, 613-627.

Datta, S., Hall, P. and Mandal, A. (2011). Model selection by testing for the presence of small-area effects, and application to area-level data. Journal of the American Statistical Association, 106, 362-374.

Diao, L., Smith, D.D., Datta, G.S., Maiti, T. and Opsomer, J.D. (2014). Accurate confidence interval estimation of small area parameters under the Fay-Herriot model. Scandinavian Journal of Statistics, to appear.

Fay, R.E., and Herriot, R.A. (1979). Estimation of income from small places: An application of James-Stein procedures to census data. Journal of the American Statistical Association, 74, 269-277.

Han, C.-P., and Bancroft, T.A. (1968). On pooling means when variance is unknown. Journal of the American Statistical Association, 63, 1333-1342.

Li, H., and Lahiri, P. (2010). An adjusted maximum likelihood method for solving small area estimation problems. Journal of Multivariate Analysis, 101, 882-892.

Prasad, N.G.N., and Rao, J.N.K. (1990). The estimation of the mean squared error of small-area estimators. Journal of the American Statistical Association, 85, 163-171.

Rao, J.N.K. (2003). Small Area Estimation. Hoboken, NJ: Wiley.

Rubin-Bleuer, S., and Yu, Y. (2013). A positive variance estimator for the Fay-Herriot small area model. SRID2-12-001E, Statistics Canada.

Saleh, A.K. Md. E. (2006). Theory of Preliminary Test and Stein-type Estimation with Applications. New York: John Wiley & Sons, Inc.

Previous

Date modified: