@follows( mkdir("plots"), getReplicatedTranscriptTSSProfileCapseq )
@transform("tss-profile/*.replicated.gene.tss-profile.capseq.counts.tsv.gz", regex(r"tss-profile/(\S+).replicated.gene.tss-profile.capseq.counts.tsv.gz"), r"plots/\1.combined.gene.tss-profile.pdf")
def plotFigure2bGene( infile, outfile):
    capseq = infile
    scriptsdir = PARAMS["scriptsdir"]
    nocapseq = capseq.replace("capseq", "nocapseq")
    R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals() )
    R('''combinedTSSPlot(capseqfile="%(capseq)s", nocapseqfile="%(nocapseq)s", outfile="%(outfile)s", ylimit=c(0,10), scale=1)''' % locals() )

############################################################
Figure 2b: TSS profiles for CAPseq and non CAPseq genes
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/plotFigure2bGene