fix: redirects to single chr. plotting if only 2 chromosomes trans in…

…teractions are found
js2264 · Oct 3, 2023 · 1e6e8cb · 1e6e8cb
1 parent 85b00a4
commit 1e6e8cb
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/R/plotting.R b/R/plotting.R
@@ -330,7 +330,10 @@ setMethod("plotMatrix", "GInteractions", function(
         GenomicRanges::seqnames(InteractionSet::regions(gis))
     )))
 
-    if (nseqnames == 1) { ## SINGLE CHROMOSOME MAP
+    if (nseqnames == 1 | {nseqnames == 2 & all(
+        GenomicRanges::seqnames(InteractionSet::anchors(gis, "first")) != 
+        GenomicRanges::seqnames(InteractionSet::anchors(gis, "second"))
+    )}) { ## SINGLE CHROMOSOME MAP or 2 CHROMOSOMES TRANS INTERSECTION
 
         if (is.null(maxDistance)) { ##### REGULAR SQUARE MATRIX
             ## -- Convert gis to table and extract x/y