Empirical best linear unbiased prediction (EBLUP) and
robust prediction of the area-level means under the basic unit-level
model. The model can be fitted by maximum likelihood or a (robust)
M-estimator. Mean square prediction error is computed by a parametric
bootstrap.