library(fields) data( ozone2) # data 16 x<- ozone2$lon.lat y<- ozone2$y[16,] # look at data quilt.plot( x,y) US( add=TRUE) # fit Tps spline tps.obj<- Tps(x,y) surface( tps.obj) US( add=TRUE) # summary of fit summary( tps.obj) # diagnostic plots set.panel( 2,2) plot( tps.obj)