R语言 | 峰峦图 / 山脊图

目的:为展示不同数据分布的差异。

1. ggplot2 实现

复制代码
# 准备数据
dat=mtcars[, c("mpg", "cyl")]
colnames(dat)=c("value", "type")
head(dat)
#                  value type
#Mazda RX4         21.0   6
#Mazda RX4 Wag     21.0   6
#Datsun 710        22.8   4

cols=c("#F71480", "#76069A", "#FF8000")
#
p1=ggplot(dat, aes(x = value, fill = as.factor(type) ) ) +
  geom_density(alpha = 0.8) +
  scale_fill_manual(values = cols)+
  facet_wrap(~type, ncol=1) +  # 按气缸数分面
  labs(title = "Density of MPG by Cylinder Count-A",
       x = "Miles Per Gallon (MPG)",
       y = "Density",
       fill = "Cylinders") +
  theme_classic(base_size = 14)+
  theme(strip.background = element_blank(),  # 去掉小标题背景
        strip.placement = "outside");p1  # 小标题外部显示
#
p2=ggplot(dat, aes(x = value, fill = as.factor(type) ) ) +
  geom_density(alpha = 0.8) +
  scale_fill_manual(values = cols)+
  facet_wrap(~type, ncol=1, scales="free_y") +  # 按气缸数分面
  labs(title = "Density of MPG by Cylinder Count-B",
       x = "Miles Per Gallon (MPG)",
       y = "Density",
       fill = "Cylinders") +
  theme_classic(base_size = 14)+
  theme(strip.background = element_blank(),  # 去掉小标题背景
        strip.placement = "outside"); p2  # 小标题外部显示
#

2. 使用R包 ggridges

图放这里,方便和上图类似。

复制代码
library(ggridges)
pB=ggplot(dat, aes(x = value, y = type, fill = factor(type, levels = c("4", "6", "8")) )) + 
  ggridges::geom_density_ridges(alpha = 0.7, show.legend = T) +
  scale_fill_manual(values = cols)+
  #scale_y_continuous( expand = c(0,0) )+
  labs(title = "Density of MPG by Cylinder Count-C",
       x = "Miles Per Gallon (MPG)",
       y = "Density",
       fill = "Cylinders") +
  theme_classic(base_size = 14); pB
#
pB2=ggplot(dat, aes(x = value, y = type, fill = factor(type, levels = c("4", "6", "8")) )) + 
  ggridges::geom_density_ridges(alpha = 0.7, show.legend = T, 
                                stat="binline", bins=25) +
  scale_fill_manual(values = cols)+
  #scale_y_continuous( expand = c(0,0) )+
  labs(title = "Density of MPG by Cylinder Count-D",
       x = "Miles Per Gallon (MPG)",
       y = "Density",
       fill = "Cylinders") +
  theme_classic(base_size = 14); pB2
#

3. 去掉底部的空隙

复制代码
pB3=ggplot(dat, aes(x = value, y = type, fill = factor(type, levels = c("4", "6", "8")) )) + 
  ggridges::geom_density_ridges(alpha = 0.7, show.legend = T, 
                                scale = 2) +
  scale_fill_manual(values = cols)+
  #scale_y_continuous( expand = c(0,0) )+
  labs(title = "Density of MPG by Cylinder Count-E\nset scale=2",
       x = "Miles Per Gallon (MPG)",
       y = "Density",
       fill = "Cylinders") +
  # 去掉底部
  scale_y_discrete(expand = c(0, 0)) +     # will generally have to set the `expand` option
  scale_x_continuous(expand = c(0, 0)) +   # for both axes to remove unneeded padding
  coord_cartesian(clip = "on") + # to avoid clipping of the very top of the top ridgeline
  theme_classic(base_size = 14); pB3

Ref

相关推荐
iCxhust2 小时前
c# U盘映像生成工具
开发语言·单片机·c#
yangzhi_emo3 小时前
ES6笔记2
开发语言·前端·javascript
emplace_back4 小时前
C# 集合表达式和展开运算符 (..) 详解
开发语言·windows·c#
jz_ddk4 小时前
[学习] C语言数学库函数背后的故事:`double erf(double x)`
c语言·开发语言·学习
萧曵 丶4 小时前
Rust 所有权系统:深入浅出指南
开发语言·后端·rust
xiaolang_8616_wjl4 小时前
c++文字游戏_闯关打怪2.0(开源)
开发语言·c++·开源
收破烂的小熊猫~4 小时前
《Java修仙传:从凡胎到码帝》第四章:设计模式破万法
java·开发语言·设计模式
nananaij5 小时前
【Python进阶篇 面向对象程序设计(3) 继承】
开发语言·python·神经网络·pycharm
阿蒙Amon5 小时前
为什么 12 版仍封神?《C# 高级编程》:从.NET 5 到实战架构,进阶者绕不开的必修课
开发语言·c#