## load observations:
library(sp)
library(rgdal)
loadMeuse()
## fit a model:
omm <- fit.gstatModel(meuse, om~dist, family=gaussian(link="log"), meuse.grid)
show(omm@regModel)
## produce SpatialPredictions:
om.rk <- predict(omm, predictionLocations = meuse.grid)
x = summary(om.rk)
str(x)Run the code above in your browser using DataLab