1Load[name = ccd, file = 'ccd_lhs_4d_42p.spd', n_predictors = 4, n_responses = 1] 2CreateAxes[name = ax4d, bounds = '0.5 1.5 | 1.5 4.5 | 0.7 1.0 | 0.7'] 3CreateSample[name = test_data, axes = ax4d, grid_points = (11,11,11,1)] 4 5# First-order polynomial 6CreateSurface[name = poly1, data = ccd, type = polynomial, order =1] 7Save[surface = poly1, file = 'poly1.sps'] 8Evaluate[surface = poly1, data = test_data, label = poly1] 9 10# MARS 11CreateSurface[name = mars, data = ccd, type = mars] 12Save[surface = mars, file = 'mars.bsps'] 13Evaluate[surface = mars, data = test_data, label = mars] 14 15# RBF 16CreateSurface[name = rbf, data = ccd, type = rbf] 17Save[surface = rbf, file = 'rbf.sps'] 18Evaluate[surface = rbf, data = test_data, label = rbf] 19 20# ANN 21CreateSurface[name = ann, data = ccd, type = ann] 22Save[surface = ann, file = 'ann.sps'] 23Evaluate[surface = ann, data = test_data, label = ann] 24 25# MLS 26CreateSurface[name = mls, data = ccd, type = mls] 27Save[surface = mls, file = 'mls.bsps'] 28Evaluate[surface = mls, data = test_data, label = mls] 29 30# Kriging 31CreateSurface[name = kriging, data = ccd, type = kriging] 32Save[surface = kriging, file = 'kriging.sps'] 33Evaluate[surface = kriging, data = test_data, label = kriging] 34 35Save[data = test_data, file = 'test_data_save.spd'] 36 37