예제 #1
0
 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)
예제 #2
0
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)
예제 #3
0
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)
예제 #4
0
 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
예제 #5
0
 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)
예제 #6
0
 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)
예제 #7
0
 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)