From e5dada202c34521618bf82a086093c342841e5e8 Mon Sep 17 00:00:00 2001
From: Thibaut Horel <thibaut.horel@gmail.com>
Date: Sun, 7 Jun 2015 22:23:39 -0700
Subject: Project cleanup before handing it to Ben

---
 experiments/plot3d.py | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 experiments/plot3d.py

(limited to 'experiments/plot3d.py')

diff --git a/experiments/plot3d.py b/experiments/plot3d.py
new file mode 100644
index 0000000..64c144d
--- /dev/null
+++ b/experiments/plot3d.py
@@ -0,0 +1,23 @@
+from mpl_toolkits.mplot3d import Axes3D
+from matplotlib import cm
+import matplotlib.pyplot as plt
+import numpy as np
+
+with open("out.log") 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[4])
+    am[0] = 1./am[0]
+    print am
+    alpha, delta, beta, _ , l = zip(*values)
+    alpha = 1./np.array(alpha)
+
+    fig = plt.figure(figsize=(12, 8))
+    ax = fig.gca(projection='3d')
+
+    ax.plot_trisurf(alpha, delta, l, cmap=cm.jet, linewidth=0.001)
+    plt.xlabel("alpha")
+    plt.ylabel("delta")
+    ax.set_zlabel('Likelihood')
+    #plt.savefig("ll.pdf")
+    plt.show()
-- 
cgit v1.2.3-70-g09d2