def plot_ld(gn, title): m = allel.rogers_huff_r(gn)**2 ax = allel.plot_pairwise_ld(m) ax.set_title(title)
def plot_ld(gn, title, filename): m = al.rogers_huff_r(gn) ** 2 ax = al.plot_pairwise_ld(m) ax.set_title(title) ax.figure.savefig(os.path.join(pcafP, filename), bbox_inches='tight')