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()