二维直方图 hexbin

install.packages("hexbin")
install.packages("RColorBrewer")
library(hexbin)
library(RColorBrewer)
x <- rnorm(mean=1.5, 5000)
y <- rnorm(mean=1.6, 5000)
bin<-hexbin(x, y, xbins=40)
my_colors=colorRampPalette(rev(brewer.pal(11,'Spectral')))
plot(bin, colramp=my_colors, legend=F)