@follows( mkdir("plots"),exportLiverTestesIntervalLengths )
@files(("liver_vs_testes/*nmi.liver.testes.shared.length","liver_vs_testes/*nmi.liver.testes.uniq.length"), "plots/"+PARAMS["species"]+".liver.testes.length.pdf")
def plotFigure3Length( infiles, outfile):
    liver_shared, testes_shared, liver_uniq, testes_uniq = infiles
    scriptsdir = PARAMS["scriptsdir"]
    R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals() )
    R('''sharesVsUniqueLengthPlot(liver_shared="%(liver_shared)s", liver_unique="%(liver_uniq)s", testes_shared="%(testes_shared)s", testes_unique="%(testes_uniq)s", outfile="%(outfile)s")''' % locals() )

############################################################
Figure 3 supplementary: length of liver and testes unique intervals compared to shared
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/plotFigure3Length