def drawNumsamVsClonesDist2(sample2dist, outfile):
dpi = 300
outformat = 'pdf'
fig, pdf = iseqlib.initImage2(10.0, 10.0, outformat, outfile, dpi)
axes = iseqlib.setAxes(fig)
drawDistData(axes, sample2dist)
iseqlib.writeImage2(fig, pdf, outformat, outfile, dpi)
def drawOverlapReads(exps, options):
sample = options.sample
options.out = os.path.join(options.outdir, "overlapPlot-%s" % sample)
fig, pdf = iseqlib.initImage( 10.0, 12.0, options )
axes = iseqlib.setAxes(fig)
drawOverlapReadsData(axes, exps, sample, options.sampleOrder)
iseqlib.writeImage(fig, pdf, options)
def drawPca(rownames, transformedM, outfile, options):
#Draw V usage distribution of each sample onto the same pdf
options.out = outfile
fig, pdf = immunoseqLib.initImage( 12.5, 8.0, options )
axes = immunoseqLib.setAxes(fig)
drawPcaData( axes, rownames, transformedM, options )
immunoseqLib.writeImage( fig, pdf, options )
def drawDist(samples, outfile, uniq):
dpi = 300
outformat = 'pdf'
fig, pdf = iseqlib.initImage2(10.0, 10.0, outformat, outfile, dpi)
axes = iseqlib.setAxes(fig)
sample2dist = getSharedSeqDist(samples, uniq)
drawDistData(axes, sample2dist)
iseqlib.writeImage2(fig, pdf, outformat, outfile, dpi)
def drawVJ( sample, vgenes, jgenes, options, minvj, maxvj ):
options.out = os.path.join( options.outdir, "%s-vjUsage" %sample.name )
fig, pdf = immunoseqLib.initImage( 8.0, 10.0, options )
axes = immunoseqLib.setAxes(fig)
#axesList = setCompareAxes( fig )
drawVJdata( fig, axes, sample, vgenes, jgenes, options, minvj, maxvj )
immunoseqLib.writeImage( fig, pdf, options )
def drawDist(sam2nt2aa, options):
options.out = os.path.join(options.outdir, 'nt2aa')
fig, pdf = iseqlib.initImage(10.0, 10.0, options)
axes = iseqlib.setAxes(fig)
drawDistData(axes, sam2nt2aa)
iseqlib.writeImage(fig, pdf, options)
def drawDist(samples, options):
options.out = os.path.join(options.outdir, "sharedSeqsDist")
fig, pdf = iseqlib.initImage(10.0, 10.0, options)
axes = iseqlib.setAxes(fig)
drawDistData(axes, samples)
iseqlib.writeImage(fig, pdf, options)