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