Browse code

Added check for presence of exon annotation

shians authored on 27/05/2024 07:45:34
Showing 1 changed files

... ...
@@ -21,7 +21,9 @@ plot_agg_genes <- function(
21 21
     span = 0.05,
22 22
     palette = ggplot2::scale_colour_brewer(palette = "Set1")
23 23
 ) {
24
-    gene_regions <- exons_to_genes(exons(x))
24
+    assertthat::assert_that(nrow(exons(x)) > 0, msg = "no exon annotations found in object")
25
+
26
+    gene_regions <- exons_to_genes(x)
25 27
     if (!is.null(genes)) {
26 28
         gene_regions <- gene_regions %>%
27 29
             filter(.data$symbol %in% genes)
... ...
@@ -49,7 +51,9 @@ plot_agg_tss <- function(
49 51
     span = 0.05,
50 52
     palette = ggplot2::scale_colour_brewer(palette = "Set1")
51 53
 ) {
52
-    gene_regions <- exons_to_genes(exons(x))
54
+    assertthat::assert_that(nrow(exons(x)) > 0, msg = "no exon annotations found in object")
55
+
56
+    gene_regions <- exons_to_genes(x)
53 57
     if (!is.null(genes)) {
54 58
         gene_regions <- gene_regions %>%
55 59
             filter(.data$symbol %in% genes)
... ...
@@ -103,7 +107,9 @@ plot_agg_tes <- function(
103 107
     span = 0.05,
104 108
     palette = ggplot2::scale_colour_brewer(palette = "Set1")
105 109
 ) {
106
-    gene_regions <- exons_to_genes(exons(x))
110
+    assertthat::assert_that(nrow(exons(x)) > 0, msg = "no exon annotations found in object")
111
+
112
+    gene_regions <- exons_to_genes(x)
107 113
     if (!is.null(genes)) {
108 114
         gene_regions <- gene_regions %>%
109 115
             filter(.data$symbol %in% genes)