Mercurial > repos > davidvanzessen > combined_immune_repertoire_pipeline

comparison RScript.r @ 8:00d432c66fb8 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 25 Nov 2013 09:21:55 -0500
parents 87fb14480352
children 8d83319a0f3d
comparison
equal deleted inserted replaced
7:e6f28d0a0f14 8:00d432c66fb8
27 27
28 28
29 test = read.table(inFile, sep="\t", header=TRUE, fill=T) 29 test = read.table(inFile, sep="\t", header=TRUE, fill=T)
30 30
31 test = test[test$Sample != "",] 31 test = test[test$Sample != "",]
32
33 if("Replicate" %in% colnames(test))
34 {
35 test$SRID = do.call(paste, c(test[c("Sample", "Replicate")], sep = "-"))
36 }
32 37
33 test$Top.V.Gene = gsub("[*]([0-9]+)", "", test$Top.V.Gene) 38 test$Top.V.Gene = gsub("[*]([0-9]+)", "", test$Top.V.Gene)
34 test$Top.D.Gene = gsub("[*]([0-9]+)", "", test$Top.D.Gene) 39 test$Top.D.Gene = gsub("[*]([0-9]+)", "", test$Top.D.Gene)
35 test$Top.J.Gene = gsub("[*]([0-9]+)", "", test$Top.J.Gene) 40 test$Top.J.Gene = gsub("[*]([0-9]+)", "", test$Top.J.Gene)
36 41