image.png
未写完......
def show_map(self):
control,treat=self.samples
fig=plt.figure(figsize=(8,4))
x_start=0.1;awidth=0.7
fig.suptitle("{}:{:,}-{:,}".format(self.chr,self.start,self.end),fontsize=14)
#draw Gene annotation.....
ax=fig.add_axes([x_start,0.1,awidth,0.2],axisbg='w')
self.adjust_spines(ax,['top','left','right','bottom'],'Genes')
Regions=["{}:{}-{}".format(self.chr,self.start,self.end),]
gene_regions=Gene_annotation(self.gtfFile,Regions)
gene_regions.draw_all_region_genes(ax)
#draw RNA....
ax=fig.add_axes([x_start,0.3,awidth,0.2])
self.adjust_spines(ax,['top','left','right','bottom'],'RNA-seq')
ax.plot([self.start,self.end],[0,0],linestyle='-', color='k', linewidth=0.5,label=self.chr,gid=self.chr)
rna_max=self.cut_rna.describe().ix['max',0]
rna_min=self.cut_rna.describe().ix['min',0]
for i,r in self.cut_rna.iterrows():
left=r['start'];height=r['Log2FoldChange'];bottom=0;
color='#00A145' if r['updown']=='up' else 'blue'
width=r['end']-r['start']
ax.bar(left=left, height=height, width=width, bottom=bottom,color=color,align="edge",edgecolor=color)
#draw treat loop
ax=fig.add_axes([x_start,0.5,awidth,0.2])
self.adjust_spines(ax,['top','left','right','bottom'],"{} loop".format(treat))
print(self.cut_loop[treat])
for i,r in self.cut_loop[treat].iterrows():
x=(r['x1']+r['y2'])/2;y=0
width=r['width'];height=r['height']
ellipse = mpatches.Ellipse((x,y), width,height,fc='#F3DFEA',ec='#A03B95',linewidth=2,alpha=0.3)
ax.add_patch(ellipse)
ax.set_xlim([self.start,self.end])
ax.set_ylim([0,1])
#draw control loop
ax=fig.add_axes([x_start,0.7,awidth,0.2])
self.adjust_spines(ax,['top','left','right','bottom'],"{} loop".format(control))
for i,r in self.cut_loop[control].iterrows():
x=(r['x1']+r['y2'])/2;y=0
width=r['width'];height=r['height']
ellipse = mpatches.Ellipse((x,y), width,height,fc='#D8EADE',ec='#00A144',linewidth=2,alpha=0.3)
ax.add_patch(ellipse)
ax.set_ylim([0,1])