fit<-Krig(ozone$x,ozone$y, Covariance="Matern",theta=50, smoothness=1)
predict.se(fit) # std errors of predictions
#
# create a grid of points
xg<- make.surface.grid(
list(East.West=seq(-15,15,,20),North.South=seq(-20,20,,20) ) )
out<- predict.se(fit,xg)
image.plot( as.surface( xg, out))
Run the code above in your browser using DataLab