@follows( mkdir("plots") )
@transform(exportPairsScoreMatrix2, regex(r"(\S+).matrix2"), r"plots/\1.score2.pdf")
def plotPairwiseScoreHeatmap2( infile, outfile ):
    R('''library(RColorBrewer)''')
    R('''library(gplots)''')
    R('''x=read.table("%(infile)s", header=TRUE)''' % locals() )
    R('''mat <- as.matrix(x[,2:dim(x)[2]])''')
    R('''rownames(mat) <- x[,1]''')
    R('''pdf("%(outfile)s", height=10, width=10)''' % locals() )
    R('''heatmap.2(mat, Rowv=TRUE, Colv=TRUE, dendrogram=c("none"), distfun=dist, hclustfun=hclust, xlab="", ylab="", key=TRUE, keysize=2, trace="none", density.info=c("none"), margins=c(10, 8), breaks=seq(0.64,1,0.02), col=rev(heat.colors(18)),colsep=seq(1,7,1),rowsep=seq(1,7,1),sepcolor="black", sepwidth=c(0.01,0.01) )''')
    R('''dev.off()''')
        
###################################################################
###################################################################
###################################################################
## Human, mouse zebrafish threeway orthologs
Join a MUC room by sending presence to it.
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_species_conservation.py/plotPairwiseScoreHeatmap2--crossover with--dataset/ETHPy150Open ralphm/wokkel/wokkel/muc.py/MUCClient.join