R|散点图+边际图（柱形图，小提琴图）—科研工具箱

本文简单的介绍2种散点图添加边际图的方法。

一 载入数据，R包

使用经典数据集iris

library(ggplot2) #加载ggplot2包
library(ggExtra)
library(ggstatsplot)
data(iris)
head(iris)

二 ggplot2 + ggExtra绘制边际散点图

使用ggplot2绘制散点图，然后利用ggExtra包的函数添加边际柱形图

2.1 绘制基础散点图

p1 <- ggplot(iris, aes(Sepal.Length, Sepal.Width)) +
 geom_point(color = \"#00AFBB\")
p1

2.2 添加一点点细节

1）添加横轴，数轴线；

2）添加R2 和 P值

3）添加回归曲线

p2 <- ggplot(iris, aes(Sepal.Length, Sepal.Width)) +
 geom_point(color = \"#00AFBB\") +
 geom_smooth(method=\"lm\", se=T) +
 geom_hline(yintercept = 3, linetype = \"dashed\", color = \"blue\") +
 geom_vline(xintercept = 6, linetype = \"dashed\", color = \"red\") +
 annotate(\"text\", x=4.5, y=4.25, parse=TRUE,

          label=\"r^2 == 0.0138 * \' p-value = 0.1529\' \")

p2

既然是ggplot2绘制的，那更多细节还不是按照需求直接加就行嘛。

2.3 添加边际条形图

使用ggMarginal添加， Type 可选参数 histogram, density 和 boxplot.

ggMarginal(p2, type = \"histogram\", fill = \"#00AFBB\")

三 ggstatsplot绘制边际散点图

直接使用ggstatsplot包的ggscatterstats函数绘制

library(ggstatsplot)

ggscatterstats(

 data = iris,                                          
 x = Sepal.Length,                                                  
 y = Sepal.Width,
 xlab = \"Sepal Length\",
 ylab = \"Sepal Width\",
 marginal = TRUE,
 marginal.type = \"densigram\",
 margins = \"both\",
 xfill = \"blue\", # 分别设置颜色
 yfill = \"#009E73\",
 title = \"Relationship between Sepal Length and Sepal Width\",
 messages = FALSE
)

其中marginal.type可选 histograms，boxplots，density，violin，densigram (density + histogram)；可自行尝试效果。

OK，文献中常见的带边际图的散点图就绘制好了！更多参数设置详见参考资料。

参考资料：

https://www.r-graph-gallery.com/277-marginal-histogram-for-ggplot2

https://indrajeetpatil.github.io/ggstatsplot/