def get_canva_genome(self) :
# plot_canva
figure, ax = plt.subplots()
addValue = 0
xticks_positions = []
chronames = []
for chromosome, dataPosi in sorted(self.data.iteritems(), key = lambda x : (len(x[0]), x[0])) :
chronames.append(chromosome)
x, y = [], []
maxposi = 0
for posi, snp in sorted(dataPosi.iteritems(), key = lambda x : x[0]) :
genome_posi = posi + addValue
snp.set_genome_posi(genome_posi)
x.append(genome_posi)
y.append(snp.abhet)
maxposi = posi if posi > maxposi else posi
ax.plot(x, y, 'o', markersize=4)
xticks_positions.append((addValue * 2 + maxposi) / 2.0)
addValue += maxposi
ax.set_xlim([0, addValue])
ax.set_xticks(xticks_positions)
ax.set_xticklabels(chronames)
canva = MplCanva(self, None, figure, ax)
canva.mpl_connect('button_press_event', self.click_pressed_genome)
self.adjust_layout()
return canva
评论列表
文章目录