fit<- Tps( BD[,1:4], BD$lnya) # fit surface to data
# surface of variables 2 and 3
# holding 1 and 4 fixed at their median levels
out.p<-predict.surface(fit, xy=c(2,3))
plot.surface(out.p, type="C") # surface plot
Run the code above in your browser using DataLab