相关文章推荐
乖乖的韭菜  ·  Electron Microscopy ...·  9 月前    · 
乖乖的韭菜  ·  nodejs gzip ...·  9 月前    · 
乖乖的韭菜  ·  HBuilderX.vue文件 以及 ...·  10 月前    · 
乖乖的韭菜  ·  CSS3变形 ...·  10 月前    · 
乖乖的韭菜  ·  Uploading to Steam ...·  10 月前    · 
神勇威武的长颈鹿  ·  python ...·  刚刚    · 
会开车的企鹅  ·  python ...·  刚刚    · 
挂过科的企鹅  ·  Java ...·  1小时前    · 

先生成随机的一组测试数据,包含两列,一列是代表group,一列是随机生成的y值。

y <- round(rnorm(1000), 1)

df <- data.frame(y = y,group = sample(c("G1", "G2", "G3", "G4", "G5"),size = 1000,replace = TRUE))

我们先来画一下普通散点图

通常,我们会首选绘制箱线图、提琴图或者柱形图等来展示整体概况,它们包含分位数或者均值等信息,以便初步观察和比较各药物治疗组的疗效的差异。但如果也想将每个个体的响应状态呈现出来,就需要在图中添加散点表示每个患者个体,参考以下示例。

ggplot(df, aes(x = group, y = y, color = group)) +

geom_point(size = 1.5, show.legend = FALSE) +  #普通散点图

scale_color_manual(values = c('#EA5C15', '#EEE938', '#E61F18', '#69B72A', '#002D8E')) +

theme(panel.grid = element_blank(), panel.background = element_blank(),

axis.line = element_line(color = 'black')) +  #去除默认的背景色、边框等

labs(x = 'Group', y = 'Response') +

stat_summary(fun.y = mean, fun.ymin = mean, fun.ymax = mean,

geom = 'crossbar', width = 0.3, size = 0.3, color = 'black') +  #计算各组均值并添加在图中

stat_summary(fun.data = function(x) median_hilow(x, 0.5),

geom = 'errorbar', width = 0.25, size = 0.2, color = 'black')  #计算各组标准差并添加在图中

下面我们测试一下抖动散点图。

抖动散点图(jitter plot)可以避免常规散点图中点过于重叠的情况,比如我们测试数据中一共1000个数据,分布比较集中,如果数据量再大,就会出现很多点重合的现象。这时候,我们就可以采用抖动散点图,它对散点添加随机的“抖动”效果,将散点适当地沿x轴刻度两侧拓展,在一定程度上表示出了局部分布。

ggplot(df, aes(x = group, y = y, color = group)) +

geom_jitter(size = 1.5, width = 0.2, show.legend = FALSE) +  #直接使用该函数绘制抖动点,width 可调整抖动幅度

#geom_point(position = 'jitter', size = 1.5, show.legend = FALSE) +  #或者在散点图中添加随机的抖动效果

scale_color_manual(values = c('#EA5C15', '#EEE938', '#E61F18', '#69B72A', '#002D8E')) +

theme(panel.grid = element_blank(), panel.background = element_blank(),

axis.line = element_line(color = 'black')) +

labs(x = 'Group', y = 'Response') +

stat_summary(fun.y = mean, fun.ymin = mean, fun.ymax = mean,

geom = 'crossbar', width = 0.3, size = 0.3, color = 'black') +  #计算各组均值并添加在图中

stat_summary(fun.data = function(x) median_hilow(x, 0.5),

geom = 'errorbar', width = 0.25, size = 0.2, color = 'black')  #计算各组标准差并添加在图中

而蜂群图((beeswarm plot)则和抖动散点图类似,也是对散点添加沿x轴刻度两侧的分散效果,但其特别之处在于,点的分散效果不是随机的。蜂群图采用了一种逻辑,以确保所绘制的点彼此靠近且不会重叠,并能更好地呈现出点分布的局部密度信息,直观而不失优雅。我们自己的感受是看上去更为紧凑。

ggplot(df, aes(x = group, y = y, color = group)) +

geom_beeswarm(cex = 1.5, show.legend = FALSE) +  #蜂群图,cex 与点的大小和分散程度有关

scale_color_manual(values = c('#EA5C15', '#EEE938', '#E61F18', '#69B72A', '#002D8E')) +  #颜色自定义赋值

theme(panel.grid = element_blank(), panel.background = element_blank(),

axis.line = element_line(color = 'black')) +  #去除默认的背景色、边框等

labs(x = 'Group', y = 'Response') +

stat_summary(fun.y = mean, fun.ymin = mean, fun.ymax = mean,

geom = 'crossbar', width = 0.3, size = 0.3, color = 'black') +  #计算各组均值并添加在图中

stat_summary(fun.data = function(x) median_hilow(x, 0.5),

geom = 'errorbar', width = 0.25, size = 0.2, color = 'black')  #计算各组标准差并添加在图中

其实,蜂群图的整理看起来就和箱型图或者小提琴图有点类似了,也就是你如果侧重整体规律就可以选择箱型图或者小提琴图,但是如果想看个体的规律,就可以选择蜂群图。另外目前箱型图或者小提琴图也是可以加入抖动效果的。

 
推荐文章