GNU R: график плотности числа точек

Иногда бывает нужно построить то, что называется диаграммой рассеяния (scatter-plot). Когда число точек на графике превышает некоторый предел, график становится почти бесполезен потому-что на нем появляется гигантская клякса (это точки (или другие символы, которые изображаются) многократно налезают друг на друга). В таких случаях предлагается строить график плотности числа точек. Сначала нужно вычислить её (хоть каким-нибудь образом), в GNU R для этого предусмотрена функция kde2d из пакета MASS. Потом полученный объект можно рисовать через contour.

Пример который рисует отличную разноцветную диаграмму:

density <- kde2d(x1,x2,lims=c(lower1,upper1,lower2,upper2))
mycol <- function(n) gray((n+1):0/(n+1))
filled.contour(density,color.palette=mycol)

filled.contour понимает title и все остальное.