R语言绘制二维密度图

1047521767

R语言绘制二维密度图二维密度图显示了两个数值变量之间的关系，一个在x轴上表示，另一个在Y轴上表示，与散点图类似，然后计算二维空间中特定区域内的观测数，并用颜色梯度表示。二维密度图有几种类型，以下主要展示二维直方图和Hexbin图的绘图代码。
0 A' g9 }6 M8 k
R语言|绘制二维直方图
( @9 r8 J- C6 z9 k$ d对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。
. B: S$ u  N1 s
1.先绘制基础的散点图；
##二维直方图
# Y* C( H. a$ Y" ~' h) ]( Q# Library
; x2 l  g) x& U( H6 T- j/ ~: Plibrary(tidyverse)
( M& L4 x: [# u$ l" S% I" Plibrary(ggplot2)
2 c# I$ O: q* t  t8 f8 n5 ?
- f: s+ Q6 V* r  f: u. E, }7 `# 读取数据
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) )
6 n0 n$ Y9 [+ Sc <- data.frame( x=rnorm(20000, 9.5, 1.9), y=rnorm(20000, 15.5, 1.9) )
data <- rbind(a,b,c)

! R* x9 H" P: F# R#绘制基础的散点图
3 c6 S, a! V8 k5 S- S' ^* Kggplot(data, aes(x=x, y=y) ) +
# l* ~: I8 _9 i; h1 `  geom_point()
! D1 z# h0 I: o# s) d6 M+ P; m3 G! X


4 a1 N! p+ V) Y% q) ^- M; J" T8 @2.绘制二维直方图。
####二维直方图可使用geom_bin2d()函数绘制
% Z, P6 l  t+ r$ V+ b7 K#二维直方图默认选项
ggplot(data, aes(x=x, y=y) ) +
  geom_bin2d() +
  theme_bw()

* R7 u" f8 C  Y$ o; C5 ?' X$ E# Bin大小控制与调色板
, P& n- {3 v1 e4 A7 w! L! P& Rggplot(data, aes(x=x, y=y) ) +
  geom_bin2d(bins = 70) +
/ H2 Q3 R# y- s  scale_fill_continuous(type = "viridis") +
  theme_bw()

R语言|绘制Hexbin图

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

# Packages
% }- C+ ?1 ]' {0 ?8 [) v' alibrary(hexbin)
library(RColorBrewer)
+ h' `$ Z3 v9 S% o" t) s0 x3 }  e
# Create data
x <- rnorm(mean=1.5, 5000)
y <- rnorm(mean=1.6, 5000)

4 ?# |- u# ]  z9 B# Make the plot
bin<-hexbin(x, y, xbins=40)
my_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
plot(bin, main="" , colramp=my_colors , legend=F )



. f! c8 G6 W9 @( T9 K! N