@@ -0,0 +1,1095 @@
+%22%22%22%0AInvestigate data in various ways%0A%22%22%22%0A%0Aimport sys%0A%0Aimport numpy as np%0Aimport matplotlib.pylab as plt%0A%0Afrom utils import extract_sig_entries%0Afrom plotter import plot_histogram%0A%0A%0Adef plot_correlation_hist(data):%0A    %22%22%22 Plot histogram of all correlations%0A    %22%22%22%0A    # gather data%0A    corrs = %5B%5D%0A    for raw_res, enh_list in data:%0A        _, raw_mat, _ = raw_res%0A%0A        if not raw_mat is None:%0A            raw_vec = extract_sig_entries(raw_mat)%0A            corrs.extend(raw_vec)%0A%0A        for enh_res in enh_list:%0A            _, enh_mat, _ = enh_res%0A%0A            if not enh_mat is None:%0A                enh_vec = extract_sig_entries(enh_mat)%0A                corrs.extend(enh_vec)%0A%0A    # plot result%0A    fig = plt.figure()%0A%0A    plot_histogram(corrs, plt.gca())%0A    plt.xlabel('simulated correlations')%0A%0A    fig.savefig('images/all_sim_corrs.pdf')%0A%0Adef main(data):%0A    %22%22%22 Analyse data%0A    %22%22%22%0A    plot_correlation_hist(data)%0A%0Aif __name__ == '__main__':%0A    if len(sys.argv) != 2:%0A        print('Usage: %25s %3Cdata file%3E' %25 sys.argv%5B0%5D)%0A        sys.exit(-1)%0A%0A    main(np.load(sys.argv%5B1%5D)%5B'data'%5D)%0A

Add further data investigations