diff options
-rw-r--r-- | calc_error.m | 13 | ||||
-rw-r--r-- | run_tests.m | 19 |
2 files changed, 19 insertions, 13 deletions
diff --git a/calc_error.m b/calc_error.m new file mode 100644 index 0000000..1f525ba --- /dev/null +++ b/calc_error.m @@ -0,0 +1,13 @@ +function calc_error(data, pred) + + quad = quad_error(data, pred); + total = total_error(data, pred); + + printf('quad mean : %9d\n', quad); + printf('total mean : %9d\n', total); + + plot((sum(abs(data - pred))), 'b'); + hold on; + plot(sort(sum(abs(data - pred))), 'r'); + hold off; +end diff --git a/run_tests.m b/run_tests.m index 4f56094..031dd82 100644 --- a/run_tests.m +++ b/run_tests.m @@ -1,16 +1,9 @@ [p, q] = get_products(); -quad = quad_error(q(29:42, :), mean_pred(q(1:28, :), 14)); -total = total_error(q(29:42, :), mean_pred(q(1:28, :), 14)); -printf('quad mean : %9d\n', quad); -printf('total mean : %9d\n', total); +% quantity sum over a day +plot(sum(q')); -quad = quad_error(p(29:42, :), regress_pred(p, q(1:28, :))); -total = total_error(p(29:42, :), regress_pred(p, q(1:28, :))); -printf('quad regress : %9d\n', quad); -printf('total regress : %9d\n', total); - -quad = quad_error(q(29:42, :), repmat(sevenday_pred(q(1:28, :), 4), 2, 1)); -total = total_error(q(29:42, :), repmat(sevenday_pred(q(1:28, :), 4), 2, 1)); -printf('quad sevenday : %9d\n', quad); -printf('total sevenday: %9d\n', total); +% plot prediction quallity +calc_error(q(29:42, :), mean_pred(q(1:28, :), 14)); +calc_error(q(29:42, :), regress_pred(p, q(1:28, :))); +calc_error(q(29:42, :), repmat(sevenday_pred(q(1:28, :), 4), 2, 1)); |