def plot_all(state, portion=1.0, figfile=None): # Print/save uncertainty report before loading pylab or creating plots draw = state.draw(portion=portion) all_vstats = var_stats(draw) print(format_vars(all_vstats)) print("\nStatistics and plots based on {nsamp:d} samples " "({psamp:.1%} of total samples drawn)".format( \ nsamp=len(draw.points), psamp=portion)) if figfile is not None: save_vars(all_vstats, figfile + "-err.json") from pylab import figure, savefig, suptitle, rcParams figext = '.' + rcParams.get('savefig.format', 'png') # histograms figure(figsize=varplot.var_plot_size(len(all_vstats))) varplot.plot_vars(draw, all_vstats) if state.title: suptitle(state.title, x=0, y=1, va='top', ha='left') if figfile is not None: savefig(figfile + "-vars" + figext) # parameter traces figure() plot_traces(state, portion=portion) suptitle("Parameter history" + (" for " + state.title if state.title else "")) if figfile is not None: savefig(figfile + "-trace" + figext) # convergence plot figure() plot_logp(state, portion=portion) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile + "-logp" + figext) # correlation plot if draw.num_vars <= 25: figure() plot_corrmatrix(draw) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile + "-corr" + figext) # parallel coordinates plot if draw.num_vars > 1: from . import parcoord figure() parcoord.plot(draw, control_var=0) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile + "-parcor" + figext)
def plot_all(state, portion=1.0, figfile=None): # Print/save uncertainty report before loading pylab or creating plots draw = state.draw(portion=portion) all_vstats = var_stats(draw) print(format_vars(all_vstats)) if figfile is not None: save_vars(all_vstats, figfile+"-err.json") from pylab import figure, savefig, suptitle, rcParams figext = '.'+rcParams.get('savefig.format', 'png') # histograms figure(figsize=varplot.var_plot_size(len(all_vstats))) varplot.plot_vars(draw, all_vstats) if state.title: suptitle(state.title, x=0, y=1, va='top', ha='left') if figfile is not None: savefig(figfile+"-vars"+figext) # parameter traces figure() plot_traces(state, portion=portion) suptitle("Parameter history" + (" for " + state.title if state.title else "")) if figfile is not None: savefig(figfile+"-trace"+figext) # R stat plot #figure() #plot_R(state, portion=portion) #if state.title: # suptitle(state.title) #if figfile is not None: # savefig(figfile+"-R"+format) # convergence plot figure() plot_logp(state, portion=portion) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile+"-logp"+figext) # correlation plot if draw.num_vars <= 25: figure() plot_corrmatrix(draw) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile+"-corr"+figext)
def plot_all(state, portion=1.0, figfile=None): from pylab import figure, savefig, suptitle, rcParams figext = '.' + rcParams.get('savefig.format', 'png') draw = state.draw(portion=portion) all_vstats = var_stats(draw) figure() plot_vars(draw, all_vstats) if state.title: suptitle(state.title) print(format_vars(all_vstats)) if figfile is not None: savefig(figfile + "-vars" + figext) if figfile is not None: save_vars(all_vstats, figfile + "-err.json") figure() plot_traces(state, portion=portion) suptitle("Parameter history" + (" for " + state.title if state.title else "")) if figfile is not None: savefig(figfile + "-trace" + figext) # Suppress R stat for now #figure() #plot_R(state, portion=portion) #if state.title: # suptitle(state.title) #if figfile is not None: # savefig(figfile+"-R"+format) figure() plot_logp(state, portion=portion) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile + "-logp" + figext) if draw.num_vars <= 25: figure() plot_corrmatrix(draw) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile + "-corr" + figext)
def plot_all(state, portion=1.0, figfile=None): from pylab import figure, savefig, suptitle, rcParams figext = '.'+rcParams.get('savefig.format', 'png') draw = state.draw(portion=portion) all_vstats = var_stats(draw) figure() plot_vars(draw, all_vstats) if state.title: suptitle(state.title) print(format_vars(all_vstats)) if figfile is not None: savefig(figfile+"-vars"+figext) figure() plot_trace(state, portion=portion) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile+"-trace"+figext) # Suppress R stat for now #figure() #plot_R(state, portion=portion) #if state.title: # suptitle(state.title) #if figfile is not None: # savefig(figfile+"-R"+format) figure() plot_logp(state, portion=portion) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile+"-logp"+figext) if draw.num_vars <= 25: figure() plot_corrmatrix(draw) if state.title: suptitle(state.title) if figfile is not None: savefig(figfile+"-corr"+figext)