Non‐parametric small area estimation using penalized spline regression

Summary.  The paper proposes a small area estimation approach that combines small area random effects with a smooth, non‐parametrically specified trend. By using penalized splines as the representation for the non‐parametric trend, it is possible to express the non‐parametric small area estimation problem as a mixed effect model regression. The resulting model is readily fitted by using existing model fitting approaches such as restricted maximum likelihood. We present theoretical results on the prediction mean‐squared error of the estimator proposed and on likelihood ratio tests for random effects, and we propose a simple non‐parametric bootstrap approach for model inference and estimation of the small area prediction mean‐squared error. The applicability of the method is demonstrated on a survey of lakes in north‐eastern USA.

[1]  H. Chernoff On the Distribution of the Likelihood Ratio , 1954 .

[2]  H. D. Patterson,et al.  Recovery of inter-block information when block sizes are unequal , 1971 .

[3]  R. Fay,et al.  Estimates of Income for Small Places: An Application of James-Stein Procedures to Census Data , 1979 .

[4]  R. N. Kackar,et al.  Approximations for Standard Errors of Estimators of Fixed and Random Effects in Mixed Linear Models , 1984 .

[5]  K. Liang,et al.  Asymptotic Properties of Maximum Likelihood Estimators and Likelihood Ratio Tests under Nonstandard Conditions , 1987 .

[6]  D. Clayton,et al.  Empirical Bayes estimates of age-standardized relative risks for use in disease mapping. , 1987, Biometrics.

[7]  Rachel M. Harter,et al.  An Error-Components Model for Prediction of County Crop Areas Using Survey and Satellite Data , 1988 .

[8]  J. Rao,et al.  The estimation of the mean squared error of small-area estimators , 1990 .

[9]  W. Scott Overton,et al.  An EPA program for monitoring ecological status and trends , 1991, Environmental monitoring and assessment.

[10]  Malay Ghosh,et al.  Small Area Estimation: An Appraisal , 1994 .

[11]  D. Stram,et al.  Variance components testing in the longitudinal mixed effects model. , 1994, Biometrics.

[12]  T. Ferguson A Course in Large Sample Theory , 1996 .

[13]  Paul H. C. Eilers,et al.  Flexible smoothing with B-splines and penalties , 1996 .

[14]  H. T. V. Vu,et al.  Generalization of likelihood ratio tests under nonstandard conditions , 1997 .

[15]  Jiming Jiang,et al.  ASYMPTOTIC PROPERTIES OF THE EMPIRICAL BLUP AND BLUE IN MIXED LINEAR MODELS , 1998 .

[16]  Bradley P. Carlin,et al.  Generalized Linear Models for Small-Area Estimation , 1998 .

[17]  P. Lahiri,et al.  A UNIFIED MEASURE OF UNCERTAINTY OF ESTIMATED BEST LINEAR UNBIASED PREDICTORS IN SMALL AREA ESTIMATION PROBLEMS , 2000 .

[18]  M. Wand,et al.  Respiratory health and air pollution: additive mixed model analyses. , 2001, Biostatistics.

[19]  M. Wand,et al.  Incorporation of historical controls using semiparametric mixed models , 2001 .

[20]  N. S. Urquhart,et al.  Designs for Evaluating Local and Regional Scale Trends , 2001 .

[21]  M. Wand,et al.  Simple Incorporation of Interactions into Additive Models , 2001, Biometrics.

[22]  D. Ruppert Selecting the Number of Knots for Penalized Splines , 2002 .

[23]  Pranab Kumar Sen,et al.  An appraisal of some aspects of statistical inference under inequality constraints , 2002 .

[24]  P. Lahiri,et al.  On measures of uncertainty of empirical Bayes small-area estimators , 2003 .

[25]  D. Ruppert,et al.  Likelihood ratio tests in linear mixed models with one variance component , 2003 .

[26]  Partha Lahiri,et al.  On the Impact of Bootstrap in Survey Sampling and Small-Area Estimation , 2003 .

[27]  Matt P. Wand,et al.  Smoothing and mixed models , 2003, Comput. Stat..

[28]  R. Little,et al.  Penalized Spline Nonparametric Mixed Models for Inference About a Finite Population Mean from Two-Stage Samples , 2003 .

[29]  Mark Von Tress,et al.  Generalized, Linear, and Mixed Models , 2003, Technometrics.

[30]  B. Ripley,et al.  Semiparametric Regression: Preface , 2003 .

[31]  Jiming Jiang,et al.  Mean squared error of empirical predictor , 2004, math/0406455.

[32]  G. Claeskens Restricted likelihood ratio lack‐of‐fit tests using mixed spline models , 2004 .

[33]  M. Wand,et al.  Exact likelihood ratio tests for penalised splines , 2005 .

[34]  Tapabrata Maiti,et al.  On parametric bootstrap methods for small area prediction , 2006 .

[35]  Jiming Jiang,et al.  Mixed model prediction and small area estimation , 2006 .

[36]  Douglas W. Nychka,et al.  FUNFITS: data analysis and statistical tools for estimating functions , 2008 .

[37]  María Dolores Ugarte,et al.  Spline smoothing in small area trend estimation and forecasting , 2009, Comput. Stat. Data Anal..