data <- read.table("asyncTable.norm.txt", sep="\t", header=TRUE, row.names=1)

numPlots <- 30
ht <- ceiling(numPlots/5)*3;
pdf(file="compareG1eAsynReps.pdf",width=15,height=ht,pointsize=10)
par(mfrow=c(7,5), cex.lab=1.5, col="darkblue", cex.axis=1.2, mar=c(5.01, 4.5, 4.01, 1.5))
names=colnames(data)
x=.1
y=.1
library("LSD")

#loop through G1E
for (i in 1:6) {
   for (k in 1:6) {
   if (i != k) {
   heatscatter(data[,i], data[,k], xlab=names[i], ylab=names[k], xlim=c(0,x), ylim=c(0,y))
   abline(a=0, b=1, col="red")
   }
   }
}
dev.off()

pdf(file="compareEr4AsynReps.pdf",width=15,height=ht,pointsize=10)
par(mfrow=c(7,5), cex.lab=1.5, col="darkblue", cex.axis=1.2, mar=c(5.01, 4.5, 4.01, 1.5))
#loop through ER4
for (m in 7:12) {
   for (n in 7:12) {
   if (m != n) {
   heatscatter(data[,m], data[,n], xlab=names[m], ylab=names[n], xlim=c(0,x), ylim=c(0,y))
   abline(a=0, b=1, col="red")
   }
   }
}
dev.off()

pdf(file="compareG1eEr4.pdf", pointsize=10)
heatscatter(data[,6], data[,12], xlab=names[6], ylab=names[12], xlim=c(0,x), ylim=c(0,y))
abline(a=0, b=1, col="red")

dev.off()