minor bugfix in cross-listing of results

bartongroup · Oct 5, 2016 · ddc04a1 · ddc04a1
1 parent 8d3e0b5
commit ddc04a1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/rats.R b/R/rats.R
@@ -259,7 +259,7 @@ call_DTU <- function(annot= NULL, TARGET_COL= "target_id", PARENT_COL= "parent_i
   with(resobj, {
     # Cross-display the DTU calls.
     if (test_transc)
-      Genes[, transc_DTU := Transcripts[, any(DTU), by = parent_id][, V1] ]
+      Genes[, transc_DTU := Transcripts[, any(DTU, na.rm=TRUE), by = parent_id][, V1] ]
     if (test_genes)
       Transcripts[, gene_DTU := merge(Genes[, .(parent_id, DTU)], Transcripts[, .(parent_id)])[, DTU] ]