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