R语言绘制二维密度图

1047521767

R语言绘制二维密度图二维密度图显示了两个数值变量之间的关系，一个在x轴上表示，另一个在Y轴上表示，与散点图类似，然后计算二维空间中特定区域内的观测数，并用颜色梯度表示。二维密度图有几种类型，以下主要展示二维直方图和Hexbin图的绘图代码。
; t3 e" o5 w; f
' ^+ x! D3 P# S: }9 I* L. oR语言|绘制二维直方图
对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。

: L# A8 r6 n8 a  D# z3 [1.先绘制基础的散点图；
##二维直方图
# Library
library(tidyverse)
library(ggplot2)
) [/ [2 q% {5 T
# 读取数据
a <- data.frame( x=rnorm(20000, 10, 1.9), y=rnorm(20000, 10, 1.2) )
b <- data.frame( x=rnorm(20000, 14.5, 1.9), y=rnorm(20000, 14.5, 1.9) )
c <- data.frame( x=rnorm(20000, 9.5, 1.9), y=rnorm(20000, 15.5, 1.9) )
data <- rbind(a,b,c)

#绘制基础的散点图
7 c9 \' m: w2 z% Nggplot(data, aes(x=x, y=y) ) +
  geom_point()
: X" U# R; X4 p/ y8 [
/ @7 x( ~5 P* q% U, H3 ^

1 a9 H6 c  f/ ?0 c! _' q0 B) v2.绘制二维直方图。
####二维直方图可使用geom_bin2d()函数绘制
- Q! u4 n5 s  x. [#二维直方图默认选项
. ^: |7 J; V* j6 }* W6 |$ P- lggplot(data, aes(x=x, y=y) ) +
  geom_bin2d() +
  theme_bw()
  b  h# p6 @* w: A& ]1 X
# Bin大小控制与调色板
* o+ G7 T0 e: A, Lggplot(data, aes(x=x, y=y) ) +
: B4 b; @" T6 N% _$ w9 ^: I  geom_bin2d(bins = 70) +
  scale_fill_continuous(type = "viridis") +
  theme_bw()
, f, u" r/ a$ I" T* ], [) s
1 H' b' c1 D4 _8 ]

R语言|绘制Hexbin图

当显示大型数据集时，散点图很难解释，因为点不可避免地会被覆盖，并且可能不被单独识别出来。hexbin图中将地块划分为多个六边形，图中阴影部分代替了条形图的高度。图形是通过R的hexbin包绘制的。

# Packages
& [  b% H7 X( c6 h. g8 ?6 m7 vlibrary(hexbin)
) i1 x( ]: y2 ~, ~/ I$ Slibrary(RColorBrewer)
3 U* ~  O% {2 l' ~0 H7 u" b. P
# Create data
x <- rnorm(mean=1.5, 5000)
7 s' t: }( y. Zy <- rnorm(mean=1.6, 5000)

5 a8 h4 \* f+ P) V9 M. e# Make the plot
bin<-hexbin(x, y, xbins=40)
/ J0 q  x+ {- Mmy_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
plot(bin, main="" , colramp=my_colors , legend=F )
; [& d  q9 m# R" P  Z