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)