-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path.Rhistory
51 lines (51 loc) · 1.26 KB
/
.Rhistory
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
devtools::document()
devtools::document()
devtools::document()
devtools::check()
devtools::check()
devtools::check()
usethis::use_package("stringr")
devtools::check()
devtools::check()
devtools::check()
devtools::document()
devtools::check()
devtools::document()
devtools::check()
?draw_volcano
library(tidyverse)
library(janitor)
library(magrittr)
library(edgeR)
library(ComplexHeatmap)
library(sva)
library(patchwork)
library(tidyverse)
draw_boxplot <- function(df, gene, groups) {
gene <- df %>%
tibble::rownames_to_column("gene_symbol") %>%
filter(gene_symbol == gene) %>%
pull(gene)
samples <- colnames(df)
df[gene, ] %>%
pivot_longer(cols = samples, names_to = "sample", values_to = "expression") %>%
mutate(groups = group_by) %>%
ggplot(aes(groups, expression)) +
geom_boxplot() +
geom_jitter(shape = 21, cex = 2.5, color = 'black', width = 0.2, height = 0, aes(fill = groups)) +
theme(legend.position = "none") +
scale_y_continuous(limits = c(0, NA)) +
labs(y = "Norm intensity") +
ggtitle(gene)
}
devtools::check()
devtools::check()
usethis::use_package("magrittr")
usethis::use_package("tidyverse")
usethis::use_package("dplyr")
usethis::use_package("tidyr")
usethis::use_package("tibble")
devtools::document()
rm(list = c("draw_boxplot"))
devtools::check()
devtools::check()