summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorBenjamin Franzke <benjaminfranzke@googlemail.com>2012-04-10 22:44:12 +0200
committerBenjamin Franzke <benjaminfranzke@googlemail.com>2012-04-10 22:44:43 +0200
commit266d15cbc2dcc897702e77a8f5d6cbc30c15fad7 (patch)
tree9b2ab47e14e3971990490d2917da259a4ac84c4e
downloaddmc-266d15cbc2dcc897702e77a8f5d6cbc30c15fad7.tar.gz
dmc-266d15cbc2dcc897702e77a8f5d6cbc30c15fad7.tar.bz2
dmc-266d15cbc2dcc897702e77a8f5d6cbc30c15fad7.zip
Add octave script to print plots of products
-rw-r--r--display_data.m29
-rwxr-xr-xsplit.sh8
2 files changed, 37 insertions, 0 deletions
diff --git a/display_data.m b/display_data.m
new file mode 100644
index 0000000..4affa96
--- /dev/null
+++ b/display_data.m
@@ -0,0 +1,29 @@
+for i = [30 45]
+ p = dlmread(sprintf("products/p%03d", i), " ");
+
+ day = p(:,1);
+ quantity = p(:,4);
+
+ middle = repmat(mean(p(:,4)), 42, 1);
+ middle_20 = repmat(mean(p(1:20, 4)), 42, 1);
+
+ [m, n] = regress(day, quantity);
+ f = @(x) m .* x + n;
+
+ [m, n] = regress(day(1:20), quantity(1:20));
+ g = @(x) m .* x + n;
+
+ figure
+ plot(day, p(:,4), "-;quantity;",
+ day, middle, "-;Mittelwert;",
+ day, middle_20, "-;Mittelwert 1-20;",
+ day, f(day), "-;Regression 1-42;",
+ day, g(day), "-;Regression 1-20;")
+
+endfor
+
+function [m, n] = regress(x, y)
+ mn = [ x ones(length(x), 1) ] \ y;
+ m = mn(1);
+ n = mn(2);
+endfunction
diff --git a/split.sh b/split.sh
new file mode 100755
index 0000000..ee9458d
--- /dev/null
+++ b/split.sh
@@ -0,0 +1,8 @@
+#!/bin/bash
+
+sed 's/|/ /g' train.txt > train.csv
+
+for i in {1..570}
+do
+ grep "^[0-9]* $i " train.csv > products/p$(printf "%03d" $i)
+done