Mercurial > repos > davidvanzessen > plotting_merged
changeset 42:af2834de221c draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Fri, 22 Nov 2013 07:51:17 -0500 |
| parents | 46c5d9d396a5 |
| children | 1b3ef7b2f4a0 |
| files | RScript.r |
| diffstat | 1 files changed, 3 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/RScript.r Fri Nov 22 07:51:11 2013 -0500 +++ b/RScript.r Fri Nov 22 07:51:17 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" , ]