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
52
53
54
55
56
57
58
59
| rm(list = ls())
library(dplyr)
library(ggplot2)
library(ggsignif)
library(ggsci)
library(extrafont)
loadfonts()
df = data.frame(value = rep(c(20:25,25:27)),
group = rep(c('A','B'),c(6,3)))
mean = aggregate(df$value, by = list(df$group), FUN = mean)
colnames(mean) = c('group','mean')
df = merge(df,mean, by = 'group')
sd = aggregate(df$value, by = list(df$group), FUN = sd)
colnames(sd) = c('group','sd')
df = merge(sd,df, by = 'group')
df$se = ifelse(df$group == 'A',
df$sd/sqrt(6),
df$sd/sqrt(3))
complist = list(c('A','B'))
p = ggplot(df) +
geom_bar(aes(group, ifelse(group == 'A',mean/6, mean/3),fill = group),
stat = 'identity', width = 0.5) +
geom_point(aes(group,value)) +
geom_errorbar(aes(group,
ymin = mean - se,
ymax = mean + se),
width = 0.2) +
geom_signif(aes(group,value),
comparisons = complist,
test = 'wilcox.test',
step_increase = 0.1,
map_signif_level = TRUE,
vjust = 0.5,
y_position = 30,
tip_length = c(0.1,0.05),
linetype = 'dashed',
color = 'black') +
geom_hline(yintercept = 35, color = 'white') +
annotate('text',1,33, label = 'Wilcox.test',
color = 'black',size = 5, family = 'Arial', face = 'plain') +
scale_y_continuous(expand = c(0,0), breaks = seq(0,36,5)) +
scale_fill_aaas() +
labs(x = '', y = 'Mean value') +
theme_bw() +
theme(legend.position = 'none',
legend.title = element_blank(),
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_text(color = 'black',size = 10,
family = 'Arial', face = 'plain'))
p
ggsave(p, filename = '2.pdf',
width = 3, height = 3.5, device = cairo_pdf)
|
