@follows( liverTestesUniqueChromatinProfile, mkdir("plots") )
@merge("liver_vs_testes/"+PARAMS["species"]+"_liver*H3K4Me3-1*profile.area.tsv.gz", r"plots/"+PARAMS["species"]+r"_liver_unique_intervals_H3K4Me3_profile.pdf")
def plotFigure3cH3K4Me3Liver( infiles, outfile):
    if len(infiles) == 2:
        inlist = "','".join(infiles)
        inlist = "'"+inlist+"'"
        scriptsdir = PARAMS["scriptsdir"]
        R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals() )
        R('''infiles <- c(%(inlist)s) ''' % locals() )
        R('''liverTestesChromatinPlot(infiles=infiles, outfile="%(outfile)s")''' % locals() )

############################################################
Figure 3c: Liver and testes H3K4Me3 reads over liver and testes unique intervals
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/plotFigure3cH3K4Me3Liver