Попробуйте ggalt
пакет
library(ggalt)
#> Loading required package: ggplot2
m <- ggplot(faithful, aes(x = eruptions, y = waiting)) +
geom_point() +
xlim(0.5, 6) +
ylim(40, 110)
m + geom_bkde2d(bandwidth=c(0.5, 4))
![](https://i.imgur.com/uvqWMPO.png)
m + stat_bkde2d(bandwidth=c(0.5, 4), aes(fill = ..level..), geom = "polygon")
![](https://i.imgur.com/gqxjo0q.png)
Или используйте значение по умолчанию из ggplot2
пакета
m + geom_density_2d()
![](https://i.imgur.com/1GdyjjE.png)
m + stat_density_2d(aes(fill = calc(level)), geom = "polygon")
![](https://i.imgur.com/00sGEv0.png)
Создано в 2018-05-10 пакетом Представ. (v0.2.0).