Merge pull request #41 from imminfo/develop

add fix for parse.mixcr for the new MiXCR format (all.x.alignment -> …
imminfo · Nov 3, 2015 · c6a0e80 · c6a0e80
2 parents b568328 + e3c66d5
commit c6a0e80
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/R/input.R b/R/input.R
@@ -641,9 +641,9 @@ parse.mixcr <- function (.filename) {
   .reads <- 'clone.count'
   .barcodes <- 'clone.count'
   .sep = '\t'
-  .vend <- "all.v.alignment"
-  .jstart <- "all.j.alignment"
-  .dalignments <- "all.d.alignment"
+  .vend <- "all.v.alignments"
+  .jstart <- "all.j.alignments"
+  .dalignments <- "all.d.alignments"
   .vd.insertions <- "VD.insertions"
   .dj.insertions <- "DJ.insertions"
   .total.insertions <- "Total.insertions"
@@ -741,7 +741,7 @@ parse.mixcr <- function (.filename) {
   .total.insertions <- "Total.insertions"
   if (recomb_type == "VJ") {
     df$Total.insertions <- 
-      sapply(strsplit(df[[.jstart]], "|", T, F, T), "[[", 4) - sapply(strsplit(df[[.vend]], "|", T, F, T), "[[", 5) - 1
+      as.numeric(sapply(strsplit(df[[.jstart]], "|", T, F, T), "[[", 4)) - as.numeric(sapply(strsplit(df[[.vend]], "|", T, F, T), "[[", 5)) - 1
   } else if (recomb_type == "VDJ") {
     df$Total.insertions <- df[[.vd.insertions]] + df[[.dj.insertions]]
   } else {