Statistical inference for linear models has classically focused on either estimation or hypothesis testing of linear combinations of fixed effects or of variance components for random effects.A third form of inference.prediction of linear combinations of fixed and random effects.has important advantages over conventional estimators in many applications.None of these approaches will result in accurate inference if the data contain strong, unaccounted for local gradients, such as spatial trends in field-plot data.Nearest neighbor methods to adjust for such trends have been widely discussed in recent literature.So far, however, these methods have been developed exclusively for classical estimation and hypothesis testing.In this article a method of obtaining nearest neighbor adjusted (NNA) predictors, along the lines of .best linear unbiased prediction,. or BLUP, is developed.A simulation study comparing .NNABLUP. to conventional NNA methods and to non-NNA alternatives suggests considerable potential for improved efficiency.