pyapemost.set_function(myloglike, myprior) if len(sys.argv) < 2: #print "SYNAPSIS: %s [calibrate|run|analyse]" % sys.argv[0] cmd = ["calibrate", "run", "analyse"] else: cmd = sys.argv[1:] if "calibrate" in cmd: pyapemost.calibrate() # run APEMoST if "run" in cmd: w = pyapemost.watch.ProgressWatcher() w.start() pyapemost.run(max_iterations = 10000, append = False) w.stop() # lets analyse the results if "analyse" in cmd: plotter = pyapemost.analyse.VisitedAllPlotter() plotter.plot() show("chain0.pdf") histograms = pyapemost.create_histograms() i = 1 plt.clf() plt.figure(figsize=(7, 4 * len(histograms))) for k,(v,stats) in histograms.iteritems(): plt.subplot(len(histograms), 1, i) plt.plot(v[:,0], v[:,2], ls='steps--', label=k) plt.legend()
pyapemost.set_function(myloglike, myprior) if len(sys.argv) < 2: #print "SYNAPSIS: %s [calibrate|run|analyse]" % sys.argv[0] cmd = ["calibrate", "run", "analyse"] else: cmd = sys.argv[1:] if "calibrate" in cmd: pyapemost.calibrate() # run APEMoST if "run" in cmd: w = pyapemost.watch.ProgressWatcher() w.start() pyapemost.run(max_iterations=10000, append=False) w.stop() # lets analyse the results if "analyse" in cmd: plotter = pyapemost.analyse.VisitedAllPlotter() plotter.plot() show("chain0.pdf") histograms = pyapemost.create_histograms() i = 1 plt.clf() plt.figure(figsize=(7, 4 * len(histograms))) for k, (v, stats) in histograms.items(): plt.subplot(len(histograms), 1, i) plt.plot(v[:, 0], v[:, 2], ls='steps--', label=k) plt.legend()
return prob pyapemost.set_function(myloglike, myprior) if len(sys.argv) < 2: # print "SYNAPSIS: %s [calibrate|run|analyse]" % sys.argv[0] cmd = ["calibrate", "run", "analyse"] else: cmd = sys.argv[1:] if "calibrate" in cmd: pyapemost.calibrate() # run APEMoST if "run" in cmd: pyapemost.run(max_iterations=100000) # lets analyse the results if "analyse" in cmd: import matplotlib.pyplot as plt histograms = pyapemost.create_histograms() i = 1 for k, (v, stats) in histograms.iteritems(): plt.subplot(len(histograms), 1, i) plt.plot(v[:, 0], v[:, 2], ls="steps--", label=k) plt.legend() print k, stats i = i + 1 plt.show()