# NOT RUN {
##load data
data(mesa.model)
data(est.mesa.model)
##find regression parameters using GLS
x.reg <- predict(mesa.model, est.mesa.model, only.pars = TRUE)
str(x.reg$pars)
# }
# NOT RUN {
##compute predictions at all locations, including beta-fields
pred.mesa.model <- predict(mesa.model, est.mesa.model,
pred.var=TRUE)
# }
# NOT RUN {
##Let's load precomputed results instead.
data(pred.mesa.model)
##study results
print(pred.mesa.model)
# }
Run the code above in your browser using DataLab