Наложение двух графиков ggplot2 stat_density2d с альфа-каналами

Я хочу наложить два графика ggplot2 с альфа-каналами таким образом, чтобы полученное изображение отображало оба набора данных. Это мои тестовые данные:

data = read.table(text="P1 -1 0 4\nP2 0 0 2\nP3 2 1 8\nP4 -2 -2 6\nP5 0.5 2 12")
data2 = read.table(text="Q1 1 1 3\nQ2 1 -1 2\nQ3 -1 1 8")
colnames(data) = c("name","x","y","score")
colnames(data2) = c("name","x","y","score")

И вот как я рисую эти данные:

ggplot(data, aes(x=x,y=y)) + 
  stat_density2d(data=data,geom="tile", aes(fill = ..density..,alpha=..density..), contour=FALSE) + 
  theme(legend.position="none") + scale_fill_gradient (low = "#FFFFFF", high = "#FF0000") + 
  xlim(-3,3) + ylim(-3,3) + 
  geom_point()

ggplot(data2, aes(x=x,y=y)) + 
  stat_density2d(data=data2,geom="tile", aes(fill = ..density..,alpha=..density..), contour=FALSE) + 
  theme(legend.position="none") + 
  scale_fill_gradient (low = "#FFFFFF", high = "#00FF00") + 
  xlim(-3,3) + ylim(-3,3) + 
  geom_point()

Первый график показывает данные, второй график данных2:

Теперь я хочу комбинацию обоих графиков. Следующее изображение - это то, что я хочу получить. Я создал его с моей программой редактирования изображений на рабочем столе, умножив оба изображения на слои.

Я попытался построить один набор данных поверх другого, но это не умножает оба уровня, а второй цвет перезаписывает первый.

ggplot(data, aes(x=x,y=y)) + 
  stat_density2d(data=data,geom="tile", aes(fill = ..density..,alpha=..density..), contour=FALSE) + 
  theme(legend.position="none") + scale_fill_gradient (low = "#FFFFFF", high = "#FF0000") + 
  xlim(-3,3) + ylim(-3,3) + 
  stat_density2d(data=data2,geom="tile", aes(fill = ..density..,alpha=..density..), contour=FALSE) + 
  scale_fill_gradient (low = "#FFFFFF", high = "#00FF00")

Кроме того, я получаю это предупреждение: Шкала для 'fill' уже присутствует. Добавление другой шкалы для "заполнения", которая заменит существующий масштаб.

Есть ли способ сделать это в R? Или есть другой способ (используя другие функции, например, smoothScatter), чтобы получить этот или аналогичный результат? Как своего рода обходной путь, я думаю, что я получу аналогичный результат с помощью ImageMagick на сервере, но я бы предпочел сделать все это в R.

Обновление 1

Умножение двух слоев выполняется в ImageMagick таким образом;

composite -compose multiply data-red.png data-green.png im-multiply.png

Это дает тот же результат, что и выше.

Обновление 2

@Роланд научила меня в своем ответе, как построить два набора данных в одном и том же сюжете. Хотя это здорово, остается одна проблема: изображение зависит от того, как вы подаете данные на график.

ggplot(rbind(data.frame(data, group="a"), data.frame(data2, group="b")), aes(x=x,y=y)) + 
  stat_density2d(geom="tile", aes(fill = group, alpha=..density..), contour=FALSE) + 
  scale_fill_manual(values=c("a"="#FF0000", "b"="#00FF00")) + 
  geom_point() + 
  theme_minimal() + 
  xlim(-3.3, 3.3) + ylim(-3.3, 3.3) +
  coord_cartesian(xlim = c(-3.2, 3.2), ylim = c(-3.2, 3.2))

дает следующий результат:

При замене порядка обоих наборов данных (теперь набор данных "b" aka data2 приходит первым, а затем данные данных aka "a" ), вы получаете аналогичный результат, но теперь доминирует красный цвет, потому что он получается позже и, следовательно, вид перезаписывает зеленые данные.

ggplot(rbind(data.frame(data2, group="a"), data.frame(data, group="b")), aes(x=x,y=y)) + 
  stat_density2d(geom="tile", aes(fill = group, alpha=..density..), contour=FALSE) + 
  scale_fill_manual(values=c("b"="#FF0000", "a"="#00FF00")) +
  geom_point() + theme_minimal() + 
  xlim(-3.3, 3.3) + ylim(-3.3, 3.3) + 
  coord_cartesian(xlim = c(-3.2, 3.2), ylim = c(-3.2, 3.2))

Мне нужны решения, которые не зависят от порядка наборов данных.

Ответ 1

Вот точно такое же решение, что и @Roland, за исключением того, что я предлагаю строку controur. Это позволит вам оценить перекрытие. Я не вижу, как geom_tile и ваша идея "умножения" могут помочь вам оценить это. Возможно, если вы используете синий и красный цвета для неперекрывающейся области и "взвешенный" фиолетовый цвет для перекрывающейся области. Но я думаю, вы должны были бы вычислить его на предыдущем шаге, прежде чем рисовать, я думаю.

ggplot(rbind(data.frame(data, group="a"), data.frame(data2, group="b")), 
       aes(x=x,y=y)) + 
  stat_density2d(geom="density2d", aes(color = group,alpha=..level..),
                 size=2,
                 contour=TRUE) + 
  #scale_color_manual(values=c("a"="#FF0000", "b"="#00FF00")) +
  geom_point() +
  theme_minimal() +
  xlim(-3.3, 3.3) + ylim(-3.3, 3.3) +
  coord_cartesian(xlim = c(-3.2, 3.2), ylim = c(-3.2, 3.2))

Ответ 2

Вы должны отображать обе плотности в одном масштабе:

ggplot(rbind(data.frame(data, group="a"), data.frame(data2, group="b")), 
       aes(x=x,y=y)) + 
  stat_density2d(geom="tile", aes(fill = group, alpha=..density..), 
                 contour=FALSE) + 
  scale_fill_manual(values=c("a"="#FF0000", "b"="#00FF00")) +
  geom_point() +
  theme_minimal() +
  xlim(-3.3, 3.3) + ylim(-3.3, 3.3) +
  coord_cartesian(xlim = c(-3.2, 3.2), ylim = c(-3.2, 3.2))

В противном случае вы увидите искаженное изображение своих данных.