Use lattice version of hex bin - hexbinplot()
. With panel
you can add your line, and with style
you can choose different ways of visualizing hexagons. Check help for hexbinplot
for more.
library(hexbin)
library(lattice)
x <- rnorm(1e6)
y <- rnorm(1e6)
hexbinplot(x ~ y, aspect = 1, bins=50,
xlab = expression(alpha), ylab = expression(beta),
style = "nested.centroids",
panel = function(...) {
panel.hexbinplot(...)
panel.abline(h=0)
})