Fix error which occurs if no FIMO matches are detected

neurogenomics · Jun 5, 2024 · 5ff0b8f · 5ff0b8f
1 parent b538c15
commit 5ff0b8f
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/R/summit_to_motif.R b/R/summit_to_motif.R
@@ -64,6 +64,13 @@ summit_to_motif <- function(peak_input,
                             bfile = bfile,
                             thresh = fimo_threshold,
                             ...)
+
+  # Return NULL lists if FIMO matches are detected
+  if (is.null(fimo_df)) return(
+    list(peak_set = NULL,
+         distance_to_summit = NULL)
+  )
+
   index_to_repeat <- base::match(as.vector(GenomicRanges::seqnames(fimo_df)),
                                  names(peaks))
   expanded_peaks <- peaks[index_to_repeat]