fit<-Tps(ozone$x,ozone$y)
predict.se(fit) # std errors of predictions
xg<- expand.grid(seq(-15,15,,40),seq(-20,20,,40))
out<- predict.se(fit,xg)
image.plot( matrix( out, 40))
Run the code above in your browser using DataLab