data(ozone2)
fit<- Tps( ozone2$lon.lat, ozone2$y[16,])
hold<- predict.derivative(fit)
set.panel(2,1)
quilt.plot( fit$x, hold[,1])
quilt.plot( fit$x, hold[,2])
# lazy plots -- better to create a grid of points and evaluate on grid.
set.panel()
Run the code above in your browser using DataLab