Расчет тепловой карты плотности

Question

Я пытаюсь построить позиционную тепловую карту, используя ggplot2, с учетом дня недели. Мой набор данных довольно большой, а вычисления очень медленные. Я начал с попытки

ggplot(aes(x=x, y=y), data=df) + facet_wrap(~day) + 
stat_density_2d(geom="raster", aes(fill=..density.., alpha=..density..))

Однако я понял, что это займет много времени, поэтому вместо этого попытался

ggplot(aes(x=x, y=y), data=df) + facet_wrap(~day) + 
stat_density_2d(geom="polygon", n=10, aes(fill=..level.., alpha=..level))

Но это все еще занимает слишком много времени для производства.

Каков наилучший подход для этого с большими наборами данных?