changeset 3:d27d745d0556 draft

Uploaded
author davidvanzessen
date Mon, 14 Oct 2013 07:20:38 -0400
parents aac19ac7cb21
children 10cfa5e9186e
files RScript.r
diffstat 1 files changed, 6 insertions(+), 6 deletions(-) [+]
line wrap: on
line diff
--- a/RScript.r	Fri Oct 11 09:11:03 2013 -0400
+++ b/RScript.r	Mon Oct 14 07:20:38 2013 -0400
@@ -111,7 +111,7 @@
 	ggplot() + 
 	geom_tile(data=dat, aes(x=factor(Top.D.Gene), y=factor(Top.V.Gene), fill=log)) + 
 	theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
-	scale_fill_gradient(low="gold", high="blue") + 
+	scale_fill_gradient(low="gold", high="blue", na.value="white") + 
 	ggtitle(unique(dat$Sample)) + 
 	xlab("D genes") + 
 	ylab("V Genes")
@@ -124,7 +124,7 @@
 completeVD = merge(VandDCount, cartegianProductVD, all.y=TRUE)
 completeVD$Length = as.numeric(completeVD$Length)
 completeVD$log = log(completeVD$Length)
-completeVD$log[is.na(completeVD$log)] = 0
+#completeVD$log[is.na(completeVD$log)] = 0
 l = split(completeVD, f=completeVD[,"Sample"])
 png("HeatmapVD%d.png", width=150+(15*length(Dchain$v.name)), height=100+(15*length(Vchain$v.name)))
 lapply(l, FUN=plotVD)
@@ -135,7 +135,7 @@
 	ggplot() + 
 	geom_tile(data=dat, aes(x=factor(Top.J.Gene), y=factor(Top.V.Gene), fill=log)) + 
 	theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
-	scale_fill_gradient(low="gold", high="blue") + 
+	scale_fill_gradient(low="gold", high="blue", na.value="white") + 
 	ggtitle(unique(dat$Sample)) + 
 	xlab("J genes") + 
 	ylab("V Genes")
@@ -147,7 +147,7 @@
 completeVJ = merge(VandJCount, cartegianProductVJ, all.y=TRUE)
 completeVJ$Length = as.numeric(completeVJ$Length)
 completeVJ$log = log(completeVJ$Length)
-completeVJ$log[is.na(completeVJ$log)] = 0
+#completeVJ$log[is.na(completeVJ$log)] = 0
 l = split(completeVJ, f=completeVJ[,"Sample"])
 png("HeatmapVJ%d.png", width=150+(15*length(Jchain$v.name)), height=100+(15*length(Vchain$v.name)))
 lapply(l, FUN=plotVJ)
@@ -157,7 +157,7 @@
 	ggplot() + 
 	geom_tile(data=dat, aes(x=factor(Top.J.Gene), y=factor(Top.D.Gene), fill=log)) + 
 	theme(axis.text.x = element_text(angle = 90, hjust = 1)) + 
-	scale_fill_gradient(low="gold", high="blue") + 
+	scale_fill_gradient(low="gold", high="blue", na.value="white") + 
 	ggtitle(unique(dat$Sample)) + 
 	xlab("J genes") + 
 	ylab("D Genes")
@@ -169,7 +169,7 @@
 completeDJ = merge(DandJCount, cartegianProductDJ, all.y=TRUE)
 completeDJ$Length = as.numeric(completeDJ$Length)
 completeDJ$log = log(completeDJ$Length)
-completeDJ$log[is.na(completeDJ$log)] = 0
+#completeDJ$log[is.na(completeDJ$log)] = 0
 l = split(completeDJ, f=completeDJ[,"Sample"])
 png("HeatmapDJ%d.png", width=150+(15*length(Jchain$v.name)), height=100+(15*length(Dchain$v.name)))
 lapply(l, FUN=plotDJ)