Browse code

Exported `plot_agg_genes()`

Shians authored on 10/10/2022 06:34:21
Showing 1 changed files

... ...
@@ -5,6 +5,8 @@
5 5
 #'
6 6
 #' @return a ggplot object containing the aggregate methylation trend of genes.
7 7
 #'
8
+#' @export
9
+#'
8 10
 #' @examples
9 11
 #' nmr <- load_example_nanomethresult()
10 12
 #' plot_agg_genes(nmr)
... ...
@@ -55,13 +57,13 @@ plot_agg_tss <- function(
55 57
     tss_regions <- gene_regions %>%
56 58
         mutate(
57 59
             start = case_when(
58
-                strand == "+" ~ start,
59
-                strand == "-" ~ end,
60
-                TRUE ~ start)
60
+                strand == "+" ~ .data$start,
61
+                strand == "-" ~ .data$end,
62
+                TRUE ~ .data$start)
61 63
         ) %>%
62 64
         mutate(
63
-            start = start - flank,
64
-            end = start + 2*flank
65
+            start = .data$start - flank,
66
+            end = .data$start + 2*flank
65 67
         )
66 68
     kb_marker <- round(flank / 1000, 1)
67 69
     labels <- c(
... ...
@@ -109,13 +111,13 @@ plot_agg_tes <- function(
109 111
     tes_regions <- gene_regions %>%
110 112
         mutate(
111 113
             start = case_when(
112
-                strand == "+" ~ end,
113
-                strand == "-" ~ start,
114
-                TRUE ~ end)
114
+                strand == "+" ~ .data$end,
115
+                strand == "-" ~ .data$start,
116
+                TRUE ~ .data$end)
115 117
         ) %>%
116 118
         mutate(
117
-            start = start - flank,
118
-            end = start + 2*flank
119
+            start = .data$start - flank,
120
+            end = .data$start + 2*flank
119 121
         )
120 122
     kb_marker <- round(flank / 1000, 1)
121 123
     labels <- c(