# NOT RUN {
# Load data set
data(paddy)
# Fit nonstationary Fay-Herriot model using sample and non-sample part of paddy data
result <- eblupNSFH2(y ~ x1+x2, var, latitude, longitude, indicator , "REML", 100, 1e-04,paddy)
result
# }
Run the code above in your browser using DataLab