@follows( mkdir("plots"),exportLiverTestesUniqueGC,exportLiverTestesSharedGC )
@files("liver_vs_testes/*.gc_content", "plots/"+PARAMS["species"]+".liver.testes.gc_content.pdf")
def plotFigure3GC( infiles, outfile):
    liver_shared, liver_uniq, testes_shared, testes_uniq = infiles
    scriptsdir = PARAMS["scriptsdir"]
    R('''source("%(scriptsdir)s/R/proj007/proj007.R") ''' % locals() )
    R('''sharesVsUniqueCpgPlot(liver_shared="%(liver_shared)s", liver_unique="%(liver_uniq)s", testes_shared="%(testes_shared)s", testes_unique="%(testes_uniq)s", outfile="%(outfile)s", xlabel="GC content", xlimit=c(0,1))''' % locals() )
    
############################################################
.. todo::
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_proj007.py/plotFigure3GC--crossover with--dataset/ETHPy150Open lisa-lab/pylearn2/pylearn2/models/dbm/ising.py/BoltzmannIsingVisible.set_biases