Почему неизвестен параметр bins для функции stat_de density2d? (ggmap)

Хорошо, добавлю это в качестве второго ответа, потому что я считаю, что описания и комментарии в первом ответе полезны, и мне не хочется их объединять. В принципе, я решил, что должен быть простой способ восстановить регрессированную функциональность. И через некоторое время, изучив некоторые основы ggplot2, я заставил это работать, переопределив некоторые ggplot2 функции:

library(ggmap)
library(ggplot2)

# -------------------------------
# start copy from stat-density-2d.R

stat_density_2d <- function(mapping = NULL, data = NULL, geom = "density_2d",
                            position = "identity", contour = TRUE,
                            n = 100, h = NULL, na.rm = FALSE,bins=0,
                            show.legend = NA, inherit.aes = TRUE, ...) {
  layer(
    data = data,
    mapping = mapping,
    stat = StatDensity2d,
    geom = geom,
    position = position,
    show.legend = show.legend,
    inherit.aes = inherit.aes,
    params = list(
      na.rm = na.rm,
      contour = contour,
      n = n,
      bins=bins,
      ...
    )
  )
}

stat_density2d <- stat_density_2d

StatDensity2d <- 
  ggproto("StatDensity2d", Stat,
          default_aes = aes(colour = "#3366FF", size = 0.5),

          required_aes = c("x", "y"),

          compute_group = function(data, scales, na.rm = FALSE, h = NULL,
                                   contour = TRUE, n = 100,bins=0) {
            if (is.null(h)) {
              h <- c(MASS::bandwidth.nrd(data$x), MASS::bandwidth.nrd(data$y))
            }

            dens <- MASS::kde2d(
              data$x, data$y, h = h, n = n,
              lims = c(scales$x$dimension(), scales$y$dimension())
            )
            df <- data.frame(expand.grid(x = dens$x, y = dens$y), z = as.vector(dens$z))
            df$group <- data$group[1]

            if (contour) {
              #  StatContour$compute_panel(df, scales,bins=bins,...) # bad dots...
              if (bins>0){
                StatContour$compute_panel(df, scales,bins)
              } else {
                StatContour$compute_panel(df, scales)
              }
            } else {
              names(df) <- c("x", "y", "density", "group")
              df$level <- 1
              df$piece <- 1
              df
            }
          }
  )

# end copy from stat-density-2d.R
# -------------------------------

set.seed(1)
n=100

df <- data.frame(x=rnorm(n, 0, 1), y=rnorm(n, 0, 1))

TestData <- ggplot (data = df) +
  stat_density2d(aes(x = x, y = y,fill = as.factor(..level..)),bins=5,geom = "polygon") +
  geom_point(aes(x = x, y = y)) +
  scale_fill_manual(values = c("yellow","red","green","royalblue", "black"))
print(TestData)

Что дает результат. Обратите внимание, что изменение параметра bins теперь дает желаемый эффект, который нельзя воспроизвести путем изменения параметра n.

Обновить:

После расширенного обсуждения с Роландом (см. Комментарии) он определил, что это, вероятно, ошибка регрессии, и отправил отчет об ошибке.

Поскольку вопрос: «Почему неизвестен параметр bins?», И я потратил немало времени на его изучение, я отвечу на него.

Ваш пример, очевидно, взят из этой ссылки от октября 2013 года, где используется параметр. Как правильно интерпретировать stat_de density2d ggplot

Однако этот параметр никогда не был задокументирован, и неясно, используется ли он там. Вероятно, это параметр, который передается другим библиотекам (например, MASS), которые используются stat_density2d.

Мы можем заставить код работать, избавившись от вызова scale_fill_manual и используя следующее:

library(ggmap)
set.seed(1)
n=100

df <- data.frame(x=rnorm(n, 0, 1), y=rnorm(n, 0, 1))

TestData <- ggplot (data = df) +
  stat_density2d(aes(x = x, y = y,fill = as.factor(..level..)), geom = "polygon",) +
  geom_point(aes(x = x, y = y)) 
#  scale_fill_manual(values = c("yellow","red","green","royalblue", "black"))
print(TestData)

что дает следующее:

Поскольку это выглядит несколько иначе, чем сюжет, опубликованный в исходной ссылке за октябрь 2013 года, я бы сказал, что с тех пор stat_density2d был сильно переписан, или, возможно, MASS:kde2d (или другая процедура MASS), и что параметр bin больше не принимается. Также непонятно, что этот параметр когда-либо делал (читайте ссылку).

Однако вы можете изменять параметры n и h - также недокументированные с точки зрения stat_density2d (насколько я могу судить).