[p, q] = get_products();

% quantity sum over a day
%plot(sum(q'));

% plot prediction quallity
calc_error('mean', q(29:42, :), mean_pred(q(1:28, :), 14));
calc_error('regress', q(29:42, :), regress_pred(p, q(1:28, :)));
% quadratic just for reference, it sucks more than mean-predicition
calc_error('quad regress', q(29:42, :), quad_regress_pred(p, q(1:28, :)));
calc_error('log regress', q(29:42, :), log_regress_pred(p, q(1:28, :)));
calc_error('sevenday', q(29:42, :), repmat(sevenday_pred(q(1:28, :), 4), 2, 1));
calc_error('random', q(29:42, :), rand_pred(q(1:28, :), 14));

[q2, removed] = remove_sevenday_frequency(q(1:14, :));
tmp = regress_pred(p([1:14 29:42],:), q2);
tmp = fft(tmp);
%tmp = (tmp .* (removed == 0)) + removed;
tmp = tmp + removed;
tmp = real(ifft(tmp));
calc_error('regress2', q(29:42, :), tmp);