From ab0b1f3cefedb35327a19ec1b6afd560bfdf802d Mon Sep 17 00:00:00 2001
From: Thibaut Horel <thibaut.horel@gmail.com>
Date: Mon, 14 Sep 2015 23:08:02 -0400
Subject: Import supplements and repo reorganization

---
 ic_experiments/old_stuff/plot.py | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)
 create mode 100644 ic_experiments/old_stuff/plot.py

(limited to 'ic_experiments/old_stuff/plot.py')

diff --git a/ic_experiments/old_stuff/plot.py b/ic_experiments/old_stuff/plot.py
new file mode 100644
index 0000000..cbea485
--- /dev/null
+++ b/ic_experiments/old_stuff/plot.py
@@ -0,0 +1,21 @@
+from mpl_toolkits.mplot3d import Axes3D
+from matplotlib import cm
+import matplotlib.pyplot as plt
+import numpy as np
+
+with open("results3.txt") as fh:
+    values = [map(float, line.strip().split()) for line in fh]
+    #values = [(b, a, l) for (b, a, l) in values if b >= 0.04]
+    am = max(values, key=lambda x: x[2])
+    print am
+    beta, alpha, l = zip(*values)
+
+    fig = plt.figure(figsize=(12, 8))
+    ax = fig.gca(projection='3d')
+
+    ax.plot_trisurf(beta, alpha, l, cmap=cm.jet, linewidth=0.001)
+    plt.xlabel("Beta")
+    plt.ylabel("Alpha")
+    ax.set_zlabel('Likelihood')
+    #plt.savefig("ll.pdf")
+    plt.show()
-- 
cgit v1.2.3-70-g09d2