def casewise_hist():
"""DEPRECATED"""
if debug:
e = test_events()
else:
e = pip2015events()
n_cases = e.events.paper.count()
fd, axarrd = subplots(n_cases)
fm, axarrm = subplots(n_cases)
fn, axarrn = subplots(n_cases)
for i, c in enumerate(e.events.paper.values):
print(c)
data = hist_data(c)
plots(data, axarrd[i], axarrm[i], axarrn[i], title=c.dtstr(), y=0.85,
fontdict={'verticalalignment': 'top', 'fontsize': 10})
del(e)
gc.collect()
for f in (fd, fm, fn):
sf.remove_subplot_gaps(f, axis='col')
fd.savefig(path.join(savedir, 'd0_cases' + tld), **savekws)
fm.savefig(path.join(savedir, 'mu_cases' + tld), **savekws)
fn.savefig(path.join(savedir, 'nw_cases' + tld), **savekws)
stats.to_csv(path.join(paths['tables'], 'hist_stats.csv'), sep='\t')
for ax in (axarrdd[-1], axarrmd[-1], axarrnd[-1]):
ax.set_title('$\\rho > %s$' % (rhomin*read.RHO_SCALE), **titlekws)
#ax.set_title('$\\rho > %s$' % rhomin, **titlekws)
ax_pdf_labeled = [axarrdd[1], axarrmd[1]]
ax_pdf_labeled.extend(list(axarra[1,:-1]))
for ax in ax_pdf_labeled:
ax.set_ylabel('PDF')
for ax in (axarrnd[1], axarra[1,2]):
ax.set_ylabel('Frequency')
for f in (fdd, fmd, fnd):
sf.remove_subplot_gaps(f, axis='col')
sf.remove_subplot_gaps(fa, axis='col', axarr=axarra)
fnd.subplots_adjust(left=0.15) # prevent ylabel cutoff
tld = '.png'
savekws = {'dpi': 400}
fdd.savefig(path.join(savedir, 'd0_rho' + tld), **savekws)
fmd.savefig(path.join(savedir, 'mu_rho' + tld), **savekws)
fnd.savefig(path.join(savedir, 'nw_rho' + tld), **savekws)
fdd.savefig(path.join(paths['paper'], 'hist_d0' + tld), **savekws)
fmd.savefig(path.join(paths['paper'], 'hist_mu' + tld), **savekws)
fnd.savefig(path.join(paths['paper'], 'hist_nw' + tld), **savekws)
# Wat?
for axarr in (axarrdd, axarrmd, axarrnd):