R语言绘制冲积图

R语言绘图

第一篇 R语言绘制会议图

前言

冲积图，也叫桑基图。它是一种特定类型的流程图，延伸的分支的宽度对应数据流量的大小，充分展现了数据的流动情况。

---

一、加载R包和数据集

本文中绘制冲击图需要用到的R包如下所示，请在使用前确保安装完成：

install.packages("ggalluvial")
install.packages("ggplot2")
library(ggalluvial)
library(ggplot2)

本文利用的数据集为Titanic,使用前先将其转化为数据框格式：

df<-as.data.frame(Titanic)
df

二、绘制冲击图

1.函数解释

• geom_alluvium()绘制冲积图
• geom_stratum()获取冲积地块地层数据集，为给定宽度的地层绘制矩形
• geom_text()向图形添加文本
• scale_x_discrete()给x轴添加标签

---

2.图形输出

ggplot(as.data.frame(Titanic),
       aes(y = Freq,
      axis1 = Class, axis2 = Sex,axis3=Age)) +
  geom_alluvium(aes(fill=Survived))+
  geom_stratum() +
  geom_text(stat = "stratum", aes(label = after_stat(stratum))) +
  scale_x_discrete(limits = c("Class", "Sex","Age"),expand = c(.3, .05))+
  xlab("Alluvia across strata")+
  theme_minimal() +
  ggtitle("Variable:",
          "Class,Sex,Age,Survived")

3.分组绘制

如果要根据Age变量分别绘制冲积图，则只要加上一行参数
facet_wrap(~ Age, scales = “free_y”)，对其他参数作相应修改

ggplot(as.data.frame(Titanic),
       aes(y = Freq,
      axis1 = Class, axis2 = Sex)) +
  geom_alluvium(aes(fill=Survived))+
  geom_stratum() +
  geom_text(stat = "stratum", aes(label = after_stat(stratum))) +
  scale_x_discrete(limits = c("Class", "Sex"),expand = c(.3, .05))+
  xlab("Alluvia across strata")+
  theme_minimal() +
  ggtitle("Variable:",
          "Class,Sex,Age,Survived")+
  facet_wrap(~ Age, scales = "free_y")

4.其它冲积图

• 五次流的时间序列凸点图

gg <- ggplot(alluvial::Refugees,
             aes(y = refugees, x = year, alluvium = country))
gg + geom_alluvium(aes(fill = country, colour = country),
                   width = 1/4, alpha = 2/3, decreasing = FALSE,
                   curve_type = "sigmoid")

• 基于国家分类的时间序列图

gg + geom_alluvium(aes(fill = country, colour = country),
                   decreasing = NA, width = 0, knot.pos = 0)

• 变量之间间距不规则

refugees_sub <- subset(alluvial::Refugees, year %in% c(2003, 2005,2007, 2010, 2013))
gg <- ggplot(data = refugees_sub,
             aes(x = year, y = refugees, alluvium = country)) +
  theme_bw() +
  scale_fill_brewer(type = "qual", palette = "Set3")
gg +
  geom_alluvium(aes(fill = country),
                alpha = .75, decreasing = FALSE, width = 1/2) +
  geom_stratum(aes(stratum = country), decreasing = FALSE, width = 1/2)

---