Scatterplot со слишком большим количеством точек

Я пытаюсь построить две переменные, где N = 700K. Проблема в том, что слишком много перекрытий, так что сюжет становится в основном сплошным блоком черного. Есть ли способ иметь "облако" в оттенках серого, где темнота сюжета зависит от количества точек в регионе? Другими словами, вместо того, чтобы показывать отдельные точки, я хочу, чтобы сюжет был "облаком", и чем больше количество точек в регионе, тем темнее эта область.

Ответ 1

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

Это легко сделать в ggplot2:

df <- data.frame(x = rnorm(5000),y=rnorm(5000))
ggplot(df,aes(x=x,y=y)) + geom_point(alpha = 0.3)

Другим удобным способом справиться с этим (и, вероятно, более подходящим для количества точек) является гексагональное бининг:

ggplot(df,aes(x=x,y=y)) + stat_binhex()

И есть также обычное старое прямоугольное биннинг (изображение опущено), которое больше похоже на вашу традиционную тепловую карту:

ggplot(df,aes(x=x,y=y)) + geom_bin2d()

Ответ 2

Вы также можете посмотреть пакет ggsubplot. Этот пакет реализует функции, которые были представлены Хэдли Уикхэмом в 2011 году (http://blog.revolutionanalytics.com/2011/10/ggplot2-for-big-data.html).

(В дальнейшем я включаю слой "points" для иллюстрации).

library(ggplot2)
library(ggsubplot)

# Make up some data
set.seed(955)
dat <- data.frame(cond = rep(c("A", "B"), each=5000),
                  xvar = c(rep(1:20,250) + rnorm(5000,sd=5),rep(16:35,250) + rnorm(5000,sd=5)),
                  yvar = c(rep(1:20,250) + rnorm(5000,sd=5),rep(16:35,250) + rnorm(5000,sd=5)))


# Scatterplot with subplots (simple)
ggplot(dat, aes(x=xvar, y=yvar)) +
  geom_point(shape=1) +
  geom_subplot2d(aes(xvar, yvar,
                     subplot = geom_bar(aes(rep("dummy", length(xvar)), ..count..))), bins = c(15,15), ref = NULL, width = rel(0.8), ply.aes = FALSE)

Тем не менее, у этого объекта есть камни, если у вас есть третья переменная для управления.

# Scatterplot with subplots (including a third variable) 

ggplot(dat, aes(x=xvar, y=yvar)) +
  geom_point(shape=1, aes(color = factor(cond))) +
  geom_subplot2d(aes(xvar, yvar,
                     subplot = geom_bar(aes(cond, ..count.., fill = cond))),
                 bins = c(15,15), ref = NULL, width = rel(0.8), ply.aes = FALSE)

Или другой подход заключается в использовании smoothScatter():

smoothScatter(dat[2:3])

Ответ 3

Смешивание альфой легко связано с базовой графикой.

df <- data.frame(x = rnorm(5000),y=rnorm(5000))
with(df, plot(x, y, col="#00000033"))

Первые шесть чисел после # - это цвет в шестнадцатеричном формате RGB, а последние два - непрозрачность, снова в шестнадцатеричном виде, поэтому 33 ~ 3/16-й непрозрачный.

Ответ 4

Вы также можете использовать контуры плотности (ggplot2):

df <- data.frame(x = rnorm(15000),y=rnorm(15000))
ggplot(df,aes(x=x,y=y)) + geom_point() + geom_density2d()

Или объединить контуры плотности с альфа-смешиванием:

ggplot(df,aes(x=x,y=y)) + 
    geom_point(colour="blue", alpha=0.2) + 
    geom_density2d(colour="black")

Ответ 5

Обзор нескольких хороших вариантов в ggplot2:

library(ggplot2)
x <- rnorm(n = 10000)
y <- rnorm(n = 10000, sd=2) + x
df <- data.frame(x, y)

Вариант А: прозрачные точки

o1 <- ggplot(df, aes(x, y)) +
  geom_point(alpha = 0.05)

Вариант B: добавить контуры плотности

o2 <- ggplot(df, aes(x, y)) +
  geom_point(alpha = 0.05) +
  geom_density_2d()

Вариант C: добавить заполненные контуры плотности

o3 <- ggplot(df, aes(x, y)) +
  stat_density_2d(aes(fill = stat(level)), geom = 'polygon') +
  scale_fill_viridis_c(name = "density") +
  geom_point(shape = '.')

Вариант D: тепловая карта плотности

o4 <- ggplot(df, aes(x, y)) +
  stat_density_2d(aes(fill = stat(density)), geom = 'raster', contour = FALSE) +       
  scale_fill_viridis_c() +
  coord_cartesian(expand = FALSE) +
  geom_point(shape = '.', col = 'white')

Вариант E: гексбины

o5 <- ggplot(df, aes(x, y)) +
  geom_hex() +
  scale_fill_viridis_c() +
  geom_point(shape = '.', col = 'white')

Вариант F: коврики

o6 <- ggplot(df, aes(x, y)) +
  geom_point(alpha = 0.1) +
  geom_rug(alpha = 0.01)

Объединить в одну фигуру:

cowplot::plot_grid(
  o1, o2, o3, o4, o5, o6,
  ncol = 2, labels = 'AUTO', align = 'v', axis = 'lr'
)

Ответ 6

Вы можете найти полезный пакет hexbin. На странице справки hexbinplot:

library(hexbin)
mixdata <- data.frame(x = c(rnorm(5000),rnorm(5000,4,1.5)),
                      y = c(rnorm(5000),rnorm(5000,2,3)),
                      a = gl(2, 5000))
hexbinplot(y ~ x | a, mixdata)

Ответ 7

Мой любимый метод построения данных такого типа - это метод, описанный в этом вопросе, - график плотности рассеяния. Идея состоит в том, чтобы сделать точечную диаграмму, но раскрасить точки их плотностью (грубо говоря, количеством совпадений в этой области).

Это одновременно:

четко показывает расположение выбросов, и
выявляет любую структуру в плотной области участка.

Вот результат из верхнего ответа на связанный вопрос: