# HG changeset patch
# User davidvanzessen
# Date 1385125973 18000
# Node ID 87fb14480352dea5880976305a24af19a858efe5
# Parent  2b1d3d6843e86209b7f4f02fc6aad8840f05590a
Uploaded

diff -r 2b1d3d6843e8 -r 87fb14480352 RScript.r
--- a/RScript.r	Fri Nov 22 08:12:46 2013 -0500
+++ b/RScript.r	Fri Nov 22 08:12:53 2013 -0500
@@ -5,6 +5,7 @@
 inFile = args[1]
 outFile = args[2]
 outDir = args[3]
+clonalType = args[4]
 
 if (!("gridExtra" %in% rownames(installed.packages()))) {
 	install.packages("gridExtra", repos="http://cran.xl-mirror.nl/") 
@@ -33,7 +34,8 @@
 test$Top.D.Gene = gsub("[*]([0-9]+)", "", test$Top.D.Gene)
 test$Top.J.Gene = gsub("[*]([0-9]+)", "", test$Top.J.Gene)
 
-test$VDJCDR3 = do.call(paste, c(test[c("Top.V.Gene", "Top.D.Gene", "Top.J.Gene","CDR3.Seq.DNA")], sep = ":"))
+#test$VDJCDR3 = do.call(paste, c(test[c("Top.V.Gene", "Top.D.Gene", "Top.J.Gene","CDR3.Seq.DNA")], sep = ":"))
+test$VDJCDR3 = do.call(paste, c(test[unlist(strsplit(clonalType, ","))], sep = ":"))
 
 PROD = test[test$VDJ.Frame != "In-frame with stop codon" & test$VDJ.Frame != "Out-of-frame" & test$CDR3.Found.How != "NOT_FOUND" , ]