显示ggplot2 hexbin图中每个bin的计数数

Question

我正在研究一个模拟失踪船只移动的项目.我使用以下代码制作了分布图:

library("ggplot2")

a <- rnorm(1000, 30.2, 2)
b <- rnorm(1000, 10, 5)
y <- (x + a + b) * 0.6

df <- data.frame(x,y)
p <- ggplot(df,aes(x=x,y=y)) +  
  ggtitle("A Priori Map") + xlab("Longtitude") + ylab("Latitude") +
  scale_fill_gradientn(colors = topo.colors(10))

p + stat_binhex(show.legend = T, bins = 20)

这会生成如下地图:

一个hexbin地图

但是,我想用一个点来显示实际的计数,而不是使用颜色显示计数的数量.因此,如果程序"落在"特定点3次,它将显示"3".

如何在R中完成？

Answer 1

以下是如何向现有图表添加计数:

library(ggplot2)
theme_set(theme_bw())

set.seed(2)
a <- rnorm(1000, 30.2, 2)
b <- rnorm(1000, 10, 5)
x = rnorm(1000)
y <- (x + a + b) * 0.6

df <- data.frame(x,y)

p <- ggplot(df,aes(x=x,y=y)) +  
  ggtitle("A Priori Map") + 
  xlab("Longtitude") + ylab("Latitude") +
  scale_fill_gradientn(colors = topo.colors(10)) +
  stat_binhex(show.legend = T, bins = 20)

p + geom_text(stat="binhex", bins=20, aes(label=..count..), show.legend=FALSE,
              colour=hcl(15,100,60), fontface="bold", size=3.5)

要删除填充颜色,您可以执行以下操作:

ggplot(df,aes(x=x,y=y)) +  
  ggtitle("A Priori Map") + 
  xlab("Longtitude") + ylab("Latitude") +
  stat_binhex(bins = 20, fill=NA, colour="black") + 
  geom_text(stat="binhex", bins=20, aes(label=..count..), colour="red")

您还可以使用文本大小突出显示最高密度的区域:

ggplot(df,aes(x=x,y=y)) +  
  ggtitle("A Priori Map") + 
  xlab("Longtitude") + ylab("Latitude") +
  stat_binhex(show.legend = T, bins = 20, fill=NA, colour="grey70") + 
  geom_text(stat="binhex", bins=20, aes(label=..count.., size=..count..), colour="red") +
  scale_size_continuous(range=c(3,6)) +
  guides(size=FALSE)

哪个也适用于没有十六进制网格:

ggplot(df,aes(x=x,y=y)) +  
  ggtitle("A Priori Map") + 
  xlab("Longtitude") + ylab("Latitude") +
  geom_text(stat="binhex", bins=20, aes(label=..count.., size=..count..), colour="red") +
  scale_size_continuous(range=c(3,6)) +
  guides(size=FALSE)