def on_plot_results(self, event): self.frame.Show(True) result = event.result[0][0] filename = result.data.sas_data.filename current_time = datetime.datetime.now().strftime("%I:%M%p, %B %d, %Y") self.parent.SetTitle(self.window_name + " - " + filename + " - " + current_time) if hasattr(result, 'convergence'): best, pop = result.convergence[:, 0], result.convergence[:, 1:] self._get_view(ConvergenceView).update(best, pop) else: self._del_view(ConvergenceView) if hasattr(result, 'uncertainty_state'): stats = var_stats(result.uncertainty_state.draw()) msg = format_vars(stats) self._get_view(CorrelationView).update(result.uncertainty_state) self._get_view(UncertaintyView).update( (result.uncertainty_state, stats)) self._get_view(TraceView).update(result.uncertainty_state) # TODO: stats should be stored in result rather than computed in bumps UncertaintyView wx.PostEvent(self.frame.parent, StatusEvent(status=msg, info="info")) else: for view in (CorrelationView, UncertaintyView, TraceView): self._del_view(view)
def process_vars(T, draw): import matplotlib.pyplot as plt from bumps.dream import stats, views title = "T=%g" % T vstats = stats.var_stats(draw) print("=== %s ===" % title) print(stats.format_vars(vstats)) plt.figure() views.plot_vars(draw, vstats) plt.suptitle(title) plt.figure() views.plot_corrmatrix(draw) plt.suptitle(title)
def process_vars(title, draw, nwalkers, plot=True, file=None): import matplotlib.pyplot as plt vstats = stats.var_stats(draw) print("=== %s ==="%title, file=file) print(stats.format_vars(vstats), file=file) if plot: plt.figure() views.plot_vars(draw, vstats) plt.suptitle(title) plt.figure() views.plot_corrmatrix(draw) plt.suptitle(title) state = State(draw, nwalkers, title) plt.figure() views.plot_logp(state)
def on_plot_results(self, event): self.frame.Show() result = event.result[0][0] if hasattr(result, 'convergence'): from bumps.gui.convergence_view import ConvergenceView best, pop = result.convergence[:, 0], result.convergence[:, 1:] self.get_panel(ConvergenceView).update(best, pop) if hasattr(result, 'uncertainty_state'): from bumps.gui.uncertainty_view import UncertaintyView, CorrelationView, TraceView from bumps.dream.stats import var_stats, format_vars stats = var_stats(result.uncertainty_state.draw()) msg = format_vars(stats) self.get_panel(CorrelationView).update(result.uncertainty_state) self.get_panel(UncertaintyView).update((result.uncertainty_state, stats)) self.get_panel(TraceView).update(result.uncertainty_state) # TODO: stats should be stored in result rather than computed in bumps UncertaintyView wx.PostEvent(self.frame.parent, StatusEvent(status=msg, info="info")) print
def on_plot_results(self, event): self.frame.Show() result = event.result[0][0] if hasattr(result, 'convergence'): best, pop = result.convergence[:, 0], result.convergence[:, 1:] self._get_view(ConvergenceView).update(best, pop) else: self._del_view(ConvergenceView) if hasattr(result, 'uncertainty_state'): stats = var_stats(result.uncertainty_state.draw()) msg = format_vars(stats) self._get_view(CorrelationView).update(result.uncertainty_state) self._get_view(UncertaintyView).update((result.uncertainty_state, stats)) self._get_view(TraceView).update(result.uncertainty_state) # TODO: stats should be stored in result rather than computed in bumps UncertaintyView wx.PostEvent(self.frame.parent, StatusEvent(status=msg, info="info")) else: for view in (CorrelationView, UncertaintyView, TraceView): self._del_view(view)
def on_plot_results(self, event): self.frame.Show(True) result = event.result[0][0] filename = result.data.sas_data.filename current_time = datetime.datetime.now().strftime("%I:%M%p, %B %d, %Y") self.parent.SetTitle(self.window_name + " - " + filename + " - " + current_time) if hasattr(result, 'convergence'): best, pop = result.convergence[:, 0], result.convergence[:, 1:] self._get_view(ConvergenceView).update(best, pop) else: self._del_view(ConvergenceView) if hasattr(result, 'uncertainty_state'): stats = var_stats(result.uncertainty_state.draw()) msg = format_vars(stats) self._get_view(CorrelationView).update(result.uncertainty_state) self._get_view(UncertaintyView).update((result.uncertainty_state, stats)) self._get_view(TraceView).update(result.uncertainty_state) # TODO: stats should be stored in result rather than computed in bumps UncertaintyView wx.PostEvent(self.frame.parent, StatusEvent(status=msg, info="info")) else: for view in (CorrelationView, UncertaintyView, TraceView): self._del_view(view)