R语言绘制二维密度图

1047521767

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

R语言|绘制二维直方图
对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。
, U! N1 J( H; @( G8 ~# |9 Y8 A
$ K- s) d- r: ^; U( \8 K1.先绘制基础的散点图；
7 {2 l7 d1 B, Q4 m) V% W##二维直方图
5 k$ G4 L1 n. ]* y# Library
library(tidyverse)
library(ggplot2)
; C$ v2 m( T; p1 Q9 Q$ T4 O
& U( U$ k$ h2 p/ E# 读取数据
4 C2 g6 q+ N5 L) i6 Ma <- data.frame( x=rnorm(20000, 10, 1.9), y=rnorm(20000, 10, 1.2) )
2 F) W0 Y/ _' U/ [1 d2 [, Fb <- 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) )
  ?* _! u( C. ?* a. d  Q+ J9 [5 xdata <- rbind(a,b,c)

#绘制基础的散点图
ggplot(data, aes(x=x, y=y) ) +
  geom_point()
  e9 k& E$ y! ?  e8 V+ {! }& @- v# E
; `* i+ D, }% h
9 u4 o  V" @; Q. ]4 A! M5 f
2.绘制二维直方图。
####二维直方图可使用geom_bin2d()函数绘制
5 h( e+ i! K. e! v6 w#二维直方图默认选项
! T6 f; Z& V$ O2 {! Y& nggplot(data, aes(x=x, y=y) ) +
4 Q; _0 F8 L. u& ^  geom_bin2d() +
  theme_bw()
, E! M5 D- e# V
# Bin大小控制与调色板
ggplot(data, aes(x=x, y=y) ) +
  geom_bin2d(bins = 70) +
  scale_fill_continuous(type = "viridis") +
  theme_bw()
4 Y! r/ o! m0 |: i

R语言|绘制Hexbin图

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

0 e% t, n0 w- ~$ q+ }# Packages
library(hexbin)
! R: M9 r6 U) \# glibrary(RColorBrewer)
; M3 i4 n. ]2 Q: K, ~% U& O
# Create data
# e. \* ~( Z  T, L7 I- Yx <- rnorm(mean=1.5, 5000)
' q0 \3 f" D0 ?6 ^4 t* @  vy <- rnorm(mean=1.6, 5000)
1 j* p# S+ k* w' T, a, Y- l
( Z% C+ \* R; x( A) t5 E1 q: ]$ d- F1 g  K# Make the plot
bin<-hexbin(x, y, xbins=40)
6 u- e& U* |1 y6 t- Cmy_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
) E  k( M, `- l! |( k+ y4 A# |/ Iplot(bin, main="" , colramp=my_colors , legend=F )

# z4 K+ m4 a$ p' d6 \% l* }
. K  R$ i- R4 C
) L, ]7 i! ?/ L; v