def plot_seismic_regression(station, phase): data = util.get_component(DATA[1:], regression.seismic_shard, [station, phase]) look_wrap = lookupWrapper.LookupWrapper(util.k_nearest_neighbors_lookup, regression.eval_seismic_metric, 50, util.nn_reduce) wrap = evalWrapper.EvalWrapper([station, phase], regression.residual_regression_setup, data, look_wrap.lookup) x, y = util.k_fold_cross_validation(10, wrap, data, regression.residual_err) filename = station + phase + 'Data.txt' util.record_eval(x, y, filename) plt.plot([i for i in range(len(x))], x, 'ro') plt.plot([i for i in range(len(y))], y, 'bo') plt.show()
def plot_plane_regression_validation(): data = [] for i in range(1000): x = 100*random.uniform(0, 1) fx = pow(x, 2) data = data + [(x, fx)] wrap = evalWrapper.EvalWrapper([], util.plane_regression_setup) x, y = util.k_fold_cross_validation(5, wrap, data, util.plane_err) plt.plot([i for i in range(len(x))], x, 'ro') plt.plot([i for i in range(len(y))], y, 'bo') plt.show()
def plot_plane_regression_validation():
data = []
for i in range(1000):
x = 100 * random.uniform(0, 1)
fx = pow(x, 2)
data = data + [(x, fx)]
wrap = evalWrapper.EvalWrapper([], util.plane_regression_setup)
x, y = util.k_fold_cross_validation(5, wrap, data, util.plane_err)
plt.plot([i for i in range(len(x))], x, 'ro')
plt.plot([i for i in range(len(y))], y, 'bo')
plt.show()
def plot_seismic_regression(station, phase):
data = util.get_component(DATA[1:], regression.seismic_shard,
[station, phase])
look_wrap = lookupWrapper.LookupWrapper(util.k_nearest_neighbors_lookup,
regression.eval_seismic_metric, 50,
util.nn_reduce)
wrap = evalWrapper.EvalWrapper([station, phase],
regression.residual_regression_setup, data,
look_wrap.lookup)
x, y = util.k_fold_cross_validation(10, wrap, data,
regression.residual_err)
filename = station + phase + 'Data.txt'
util.record_eval(x, y, filename)
plt.plot([i for i in range(len(x))], x, 'ro')
plt.plot([i for i in range(len(y))], y, 'bo')
plt.show()