Rely on hexbin() instead of plot()<br><br>hb &lt;- hexbin(your_independent_var, your_dependent_var, xbins = 40)<br>plot(hb)<br><br>Duccio<br><br>