@follows( overlappedGeneCAPseqProfile, controlGeneCAPseqProfile, mkdir("plots") )
@merge("overlapped_genes/"+PARAMS["species"]+"_"+PARAMS["plots_fig4_tissue"]+"*.capseq_profile.counts.tsv.gz", "plots/"+PARAMS["species"]+"_"+PARAMS["plots_fig4_tissue"]+"_overlapped_genes_capseq_profile.pdf")
def plotFigure4a( infiles, outfile):
    overlapped,control = infiles
    species = overlapped[0:2]
    scriptsdir = PARAMS["scriptsdir"]
    R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals() )
    R('''overlappedGenesProfilePlot(overlapped="%(overlapped)s", control="%(control)s", outfile="%(outfile)s", ylabel="NMIs")''' % locals() )

############################################################
Start the service
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/plotFigure4a--crossover with--dataset/ETHPy150Open phensley/gstatsd/gstatsd/service.py/StatsDaemon.start