fix error when intron not on same chr as any genes in the txdb (issue #…

…74)
gagneurlab · May 2, 2024 · 3db74d9 · 3db74d9
1 parent e775832
commit 3db74d9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/resultAnnotations.R b/R/resultAnnotations.R
@@ -815,7 +815,7 @@ checkIntergenic <- function(junctions_gr, i, refseq.genes){
     # check if distance to nearest is > 1000 -> intergenic
     # otherwise up/downstream
     dist = min(distance(test_junction, refseq.genes), na.rm = TRUE)
-    if(dist > 0){
+    if(is.finite(dist) && dist > 0){
         # find nearest and compare starts
         if(start(refseq.genes[nearest(junctions_gr[i], 
                                         refseq.genes)]) > start){