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 6Load[name = poly1, file = 'poly1.sps'] 7Evaluate[surface = poly1, data = test_data, label = poly1] 8 9# MARS 10Load[name = mars, file = 'mars.bsps'] 11Evaluate[surface = mars, data = test_data, label = mars] 12 13# RBF 14Load[name = rbf, file = 'rbf.sps'] 15Evaluate[surface = rbf, data = test_data, label = rbf] 16 17# ANN 18Load[name = ann, file = 'ann.sps'] 19Evaluate[surface = ann, data = test_data, label = ann] 20 21# MLS 22Load[name = mls, file = 'mls.bsps'] 23Evaluate[surface = mls, data = test_data, label = mls] 24 25# Kriging 26Load[name = kriging, file = 'kriging.sps'] 27Evaluate[surface = kriging, data = test_data, label = kriging] 28 29Save[data = test_data, file = 'test_data_load.spd'] 30