R语言绘制二维密度图

1047521767

R语言绘制二维密度图二维密度图显示了两个数值变量之间的关系，一个在x轴上表示，另一个在Y轴上表示，与散点图类似，然后计算二维空间中特定区域内的观测数，并用颜色梯度表示。二维密度图有几种类型，以下主要展示二维直方图和Hexbin图的绘图代码。

R语言|绘制二维直方图
7 e$ X3 l; ~4 t9 y$ w. m- c* U对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。
8 L" @+ ~: |; ~& ^9 H
1.先绘制基础的散点图；
##二维直方图
/ p3 I# M" m& u5 A: c# Library
library(tidyverse)
library(ggplot2)
- f6 m* \1 E& d  q
$ p" q7 W# D! v: R0 c# 读取数据
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) )
* a4 k+ Q! g4 n7 @2 `data <- rbind(a,b,c)
$ [2 P+ \6 M! g% G; S; A
. p" r- D/ L, T" t* Q#绘制基础的散点图
ggplot(data, aes(x=x, y=y) ) +
  geom_point()
" G% p# V" y* V, K& `5 Q& T/ f
7 C% r- x& @$ X% ?+ C+ M

2.绘制二维直方图。
+ E. F* z. f1 I' [####二维直方图可使用geom_bin2d()函数绘制
#二维直方图默认选项
ggplot(data, aes(x=x, y=y) ) +
2 ?2 A9 @7 Q* D  geom_bin2d() +
; g+ T0 ^3 n/ q% S! ^  theme_bw()
- x' z& k5 ^6 d* r
# Bin大小控制与调色板
" a- D4 t+ z. N! G1 k2 f) ?ggplot(data, aes(x=x, y=y) ) +
3 j4 d! D: t( h$ K9 q! R  geom_bin2d(bins = 70) +
  scale_fill_continuous(type = "viridis") +
# {  P7 S" @0 l0 n# V* L& G  theme_bw()

' {# A: r( i! F0 F
! ~) ^! n  Z* m$ \  d/ m& b

R语言|绘制Hexbin图

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

# Packages
library(hexbin)
library(RColorBrewer)

3 N9 }- r; p- P6 a. _; W# Create data
% s: X- W) K8 Ox <- rnorm(mean=1.5, 5000)
$ {1 X% L5 o' ~1 F! by <- rnorm(mean=1.6, 5000)

# Make the plot
bin<-hexbin(x, y, xbins=40)
4 J- U, D7 a9 }+ r- q% f3 qmy_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
plot(bin, main="" , colramp=my_colors , legend=F )
" {6 l0 R8 q/ Z% c) B2 M+ r

5 K0 _" ^, o5 N% h7 g' j- N