def main(data_dir, mtype_file):  
<EOL>stuff = load_neurite_features(data_dir)<EOL>sim_params = json.load(open(mtype_file))<EOL>_plots = []<EOL>for feat, d in stuff.items():<EOL><INDENT>for typ, data in d.items():<EOL><INDENT>dist = sim_params['<STR_LIT>'][typ].get(feat, None)<EOL>print('<STR_LIT>' % (typ, feat, dist))<EOL>if not data:<EOL><INDENT>print("<STR_LIT>" % (feat, typ))<EOL>continue<EOL><DEDENT>num_bins = <NUM_LIT:100><EOL>limits = calc_limits(data, dist)<EOL>bin_edges = np.linspace(limits[<NUM_LIT:0>], limits[<NUM_LIT:1>], num_bins + <NUM_LIT:1>)<EOL>histo = np.histogram(data, bin_edges, normed=True)<EOL>print('<STR_LIT>', limits)<EOL>plot = Plot(*view_utils.get_figure(new_fig=True, subplot=<NUM_LIT>))<EOL>plot.ax.set_xlim(*limits)<EOL>plot.ax.bar(histo[<NUM_LIT:1>][:-<NUM_LIT:1>], histo[<NUM_LIT:0>], width=bin_widths(histo[<NUM_LIT:1>]))<EOL>dp, bc = dist_points(histo[<NUM_LIT:1>], dist)<EOL>if dp is not None:<EOL><INDENT>plot.ax.plot(bc, dp, '<STR_LIT>')<EOL><DEDENT>plot.ax.set_title('<STR_LIT>' % (feat, typ))<EOL>_plots.append(plot)<EOL><DEDENT><DEDENT>return _plots<EOL>
Run the stuff