R语言绘制二维密度图

1047521767

R语言绘制二维密度图二维密度图显示了两个数值变量之间的关系，一个在x轴上表示，另一个在Y轴上表示，与散点图类似，然后计算二维空间中特定区域内的观测数，并用颜色梯度表示。二维密度图有几种类型，以下主要展示二维直方图和Hexbin图的绘图代码。
7 c1 Z7 U# F$ @: g' \% m
R语言|绘制二维直方图
对于2d直方图，plot区域被划分为多个正方形，调用geombin2d（）函数，它提供了一个bins参数，用于控制要显示的bins数量。

% c: C# }5 d# w1.先绘制基础的散点图；
##二维直方图
# Library
library(tidyverse)
library(ggplot2)
, t# f, S+ ~6 Z9 F' F1 I& F* }
1 [, G0 ^5 L7 y4 Z# 读取数据
8 K5 C2 U1 ?# p8 U. i9 ]; {6 r0 [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) )
' f- i' E# @! s: Y; w- fc <- data.frame( x=rnorm(20000, 9.5, 1.9), y=rnorm(20000, 15.5, 1.9) )
data <- rbind(a,b,c)

9 P6 t  i. A# k* r: h#绘制基础的散点图
5 R  u; P: g, Y& T$ vggplot(data, aes(x=x, y=y) ) +
  geom_point()



/ w! e+ b& H, \' l. {0 m6 R2.绘制二维直方图。
' D) c2 z7 T5 k! J; y( a####二维直方图可使用geom_bin2d()函数绘制
#二维直方图默认选项
ggplot(data, aes(x=x, y=y) ) +
3 Z; k  K* L7 D  geom_bin2d() +
  theme_bw()

$ N0 D6 @/ [: y9 d3 G2 f% G- [# Bin大小控制与调色板
ggplot(data, aes(x=x, y=y) ) +
3 e' u5 t& M# Y  }  geom_bin2d(bins = 70) +
% H5 D+ F# l4 z/ L# L, q  scale_fill_continuous(type = "viridis") +
, s- G5 }) d# E0 a2 q  theme_bw()

2 t- _) V* Y4 n* D+ v

R语言|绘制Hexbin图

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

; c" w) \% o3 e/ o# Packages
library(hexbin)
library(RColorBrewer)

# Create data
x <- rnorm(mean=1.5, 5000)
4 A( U! i8 R* l8 \) g% by <- rnorm(mean=1.6, 5000)
! M$ _6 M5 r  `! o& E2 {' y' q1 q+ B, z
# Make the plot
bin<-hexbin(x, y, xbins=40)
my_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
; @. P2 V( n2 n4 I! [plot(bin, main="" , colramp=my_colors , legend=F )

6 d. B+ k4 s# ?2 H% r7 J