In plot.hexbin, if the legend is not shown the previous plot is not cleared. A fix seems to be to add if (legend) before par(new = TRUE, mai=mai, fig = fig) about five lines from the bottom of the function. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlumley at u.washington.edu University of Washington, Seattle