示例#1
0
def plot_ld(gn, title):
    m = allel.rogers_huff_r(gn)**2
    ax = allel.plot_pairwise_ld(m)
    ax.set_title(title)
示例#2
0
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')