def versus(data, x, y, xlabel="ratio", ylabels=("growth1", "growth2"), outfn="versus"): plt.figure(figsize=(12, 6)) re_data = reconfigure_data(data, x, y) # series1 = data[(data[x] > 1)] # series2 = data[(data[x] < 1)] # # growth1 = pd.concat([series1[y[0]], series2[y[1]]]) # growth2 = pd.concat([series1[y[1]], series2[y[0]]]) # ratio = pd.concat([series1[x], 1 / series2[x]]) # # indata = {} # indata[xlabel] = ratio # indata[ylabels[0]] = growth1 # indata[ylabels[1]] = growth2 # re_data = pd.DataFrame(indata) plt.subplot(131) sns.regplot(x="x", y="y1", data=re_data) shared.add_stats(re_data, "x", "y1") sns.despine() plt.subplot(132) sns.regplot(x="x", y="y2", data=re_data) shared.add_stats(re_data, "x", "y2") sns.despine() ax = plt.subplot(133) swarm(ax, re_data, ylabels[0], ylabels[1], "growth rate (\si{\per\hour})") plt.tight_layout() plt.savefig("ParA_inheritance/{0}.pdf".format(outfn))
def plot_representation(data, spnum): areas = pd.concat([data.child1_area, data.child2_area]) totals = pd.concat([data.child1_total, data.child2_total]) maxes = pd.concat([data.child1_max, data.child2_max]) d = pd.DataFrame({"area": areas, "total": totals, "max": maxes}) plt.subplot(spnum) sns.regplot(x="area", y="total", data=d) shared.add_stats(d, "area", "total", m=False) plt.xlabel("Cell Area (um$^2$)") plt.ylabel("Total ParA Intensity (AU)") sns.despine() plt.subplot(spnum + 1) sns.regplot(x="area", y="max", data=d) shared.add_stats(d, "area", "max", m=False) plt.xlabel("Cell Area (um$^2$)") plt.ylabel("Maximum ParA Intensity (AU)") sns.despine()
def plot_ratios(data, outfn): plt.figure() plt.subplot(221) plt.title("Total Intensity Ratio") sns.despine() sns.regplot(x="area_ratio", y="child_ratio", data=data) # plot regression line shared.add_stats(data, "area_ratio", "child_ratio") plt.xlabel("area\_ratio") plt.ylabel("child\_ratio") plt.legend(loc=2) # plt.plot(data.length_ratio, data.child_ratio, ls="none", marker=".") plt.subplot(222) plt.title("Max Intensity Ratio") sns.despine() sns.regplot(x="area_ratio", y="max_ratio", data=data) shared.add_stats(data, "area_ratio", "max_ratio") plt.xlabel("area\_ratio") plt.ylabel("max\_ratio") plt.legend(loc=4) plt.subplot(223) sns.despine() area_rel = data.child_ratio / data.area_ratio sns.distplot(area_rel) plt.axvline(1, color="k") plt.xlabel("child\_ratio / area\_ratio") plt.ylabel("Frequency") plt.subplot(224) sns.despine() length_rel = data.max_ratio / data.area_ratio sns.distplot(length_rel) plt.axvline(1, color="k") plt.xlabel("max\_ratio / area\_ratio") plt.ylabel("Frequency") plt.tight_layout() plt.savefig("ParA_inheritance/{0}.pdf".format(outfn))
def sanity_check(data): plt.figure() plt.subplot(221) plt.title("Total Intensity Ratio") sns.despine() sns.regplot(x="area_ratio", y="child_ratio", data=data) shared.add_stats(data, "area_ratio", "child_ratio") plt.xlabel("area\_ratio") plt.ylabel("child\_ratio") plt.legend(loc=2) plt.subplot(222) plt.title("Max Intensity Ratio") sns.despine() sns.regplot(x="area_ratio", y="max_ratio", data=data) shared.add_stats(data, "area_ratio", "max_ratio") plt.xlabel("area\_ratio") plt.ylabel("max\_ratio") plt.legend(loc=4) plt.subplot(223) sns.despine() area_rel = data.child_ratio / data.area_ratio sns.distplot(area_rel) plt.axvline(1, color="k") plt.xlabel("child\_ratio / area\_ratio") plt.ylabel("Frequency") plt.subplot(224) sns.despine() length_rel = data.max_ratio / data.area_ratio sns.distplot(length_rel) plt.axvline(1, color="k") plt.xlabel("max\_ratio / area\_ratio") plt.ylabel("Frequency") plt.tight_layout() plt.savefig("ParA_inheritance/sanity-check.pdf")
def diffusion(spot_data): mid = diffusion_data(spot_data, "x_mid") new = diffusion_data(spot_data, "x_new") old = diffusion_data(spot_data, "x_old") plt.figure() plt.subplot(131) plt.title("Mid") sns.regplot(x="delta_t", y="delta_x2", data=mid) shared.add_stats(mid, "delta_t", "delta_x2") plt.subplot(132) plt.title("New") sns.regplot(x="delta_t", y="delta_x2", data=new) shared.add_stats(new, "delta_t", "delta_x2") plt.subplot(133) plt.title("Old") sns.regplot(x="delta_t", y="delta_x2", data=old) shared.add_stats(old, "delta_t", "delta_x2") plt.tight_layout() plt.show()