R语言绘制二维密度图

1047521767

R语言绘制二维密度图二维密度图显示了两个数值变量之间的关系，一个在x轴上表示，另一个在Y轴上表示，与散点图类似，然后计算二维空间中特定区域内的观测数，并用颜色梯度表示。二维密度图有几种类型，以下主要展示二维直方图和Hexbin图的绘图代码。
7 x" o( `' w4 g. S0 W- n. V) F
R语言|绘制二维直方图
对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。
/ I* v; Z5 r! e" o2 p! P' r8 t
1.先绘制基础的散点图；
##二维直方图
1 f2 I" a" n/ M) n% ~3 U, p9 j7 Y# Library
4 m/ D* Z8 N' o6 Jlibrary(tidyverse)
library(ggplot2)

/ ^. K9 B( }: h8 J5 E# 读取数据
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)

#绘制基础的散点图
+ `# Y) b! N/ i8 ]4 K; Nggplot(data, aes(x=x, y=y) ) +
4 J  a. O. b  x7 s2 e% ]# b  geom_point()
5 N" `7 U; o! ^0 O$ d( A


2.绘制二维直方图。
: M" G' g" ]( I7 Z% o; B$ k, f+ Y####二维直方图可使用geom_bin2d()函数绘制
#二维直方图默认选项
, B' d0 V! {6 F, e5 E" j  g" kggplot(data, aes(x=x, y=y) ) +
8 \9 r; j' c1 A9 u  geom_bin2d() +
$ c  e: C$ G5 y1 Q, L" |  theme_bw()
! t- q  ~9 b: u* w6 M# M
. M% k1 R' E, T4 z6 s: @2 d# Bin大小控制与调色板
9 X+ ~+ B% V+ b* Yggplot(data, aes(x=x, y=y) ) +
  geom_bin2d(bins = 70) +
  scale_fill_continuous(type = "viridis") +
  theme_bw()
: c! Z9 N8 G! w  l

R语言|绘制Hexbin图

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

# Packages
library(hexbin)
3 [3 H- O( Y* D4 `9 D$ I* wlibrary(RColorBrewer)
9 M2 H& S$ R2 C4 |7 T, H
4 B$ g2 S* F9 Z  E, F( _' L, I# Create data
% \( ]( T+ z. P+ J9 B: dx <- rnorm(mean=1.5, 5000)
y <- rnorm(mean=1.6, 5000)

# Make the plot
2 f" v. y7 g7 a$ Qbin<-hexbin(x, y, xbins=40)
my_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
( j( k1 g1 K) S' B9 q; ~2 c1 _7 p1 Xplot(bin, main="" , colramp=my_colors , legend=F )

! i$ b3 ?# ~! E/ Q