blob: cbea48548d8dfcc425fabf55a9357155bd575bee (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
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()
|
