R语言绘制二维密度图

1047521767

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

$ _( s; _- B4 N: l$ P+ [R语言|绘制二维直方图
对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。
+ I( W4 n7 {" ~8 N6 [8 x
7 C5 l* ]" E# G$ O" r. B' _1.先绘制基础的散点图；
) r, q3 Z+ y/ z6 S##二维直方图
# Library
library(tidyverse)
2 P! G- T$ w2 K% I* Y5 u1 Elibrary(ggplot2)
) l5 c  ]4 k! ?" x
2 M, Q( U# ?* ^  ?7 z- Y5 B# 读取数据
  h# z& `6 C" ~) E" aa <- 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) )
3 K" {& ~* j- [$ [c <- data.frame( x=rnorm(20000, 9.5, 1.9), y=rnorm(20000, 15.5, 1.9) )
data <- rbind(a,b,c)

2 [% D2 B! x0 ]; \, _1 M#绘制基础的散点图
- B3 b+ u  C- Y) I( B  Y  e5 uggplot(data, aes(x=x, y=y) ) +
  geom_point()
! q9 K6 V5 z8 M% s


2.绘制二维直方图。
####二维直方图可使用geom_bin2d()函数绘制
#二维直方图默认选项
) _1 z) t* v# Z: _ggplot(data, aes(x=x, y=y) ) +
9 N  L; F1 F8 ?  geom_bin2d() +
  theme_bw()

# Bin大小控制与调色板
9 T) C( \- H% ~+ y! ~ggplot(data, aes(x=x, y=y) ) +
) [) w: P+ e* u- L0 k, ?' q& v( P  geom_bin2d(bins = 70) +
  scale_fill_continuous(type = "viridis") +
. e9 L' o( K, w  theme_bw()

! O9 J* B' p' h" ?* u
( Y& V( ^5 m" V9 p7 Q: J# P

R语言|绘制Hexbin图

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

# Packages
6 C# u* j6 X# Xlibrary(hexbin)
library(RColorBrewer)

7 X! z  R4 X5 @! \  J: c5 g/ b# Create data
x <- rnorm(mean=1.5, 5000)
2 z; i* S  Z3 C- m3 n9 O4 @y <- rnorm(mean=1.6, 5000)
6 W0 p! ?3 y6 _/ h0 |
% G* u& {7 H5 x. _/ v- R5 C# Make the plot
bin<-hexbin(x, y, xbins=40)
! q2 Q0 F. }0 T8 lmy_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
plot(bin, main="" , colramp=my_colors , legend=F )
  i7 Q+ I- e6 F8 s4 k: j; s
' R0 i4 C! a2 S! h

8 v4 ?* k+ l: `" \