1

我已经坚持了几个小时。当我运行这个：

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..)),bins=4, geom = "polygon",) +
  geom_point(aes(x = x, y = y)) +
  scale_fill_manual(values = c("yellow","red","green","royalblue", "black"))

我收到此错误消息：

Error: Unknown parameters: bins

有谁知道为什么？

4

2 回答 2

3

好的，添加这个作为第二个答案，因为我认为第一个答案中的描述和评论很有用，我不想合并它们。基本上我认为必须有一种简单的方法来恢复回归的功能。过了一会儿，学习了一些关于 的基础知识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参数来复制。

在此处输入图像描述

于 2015-12-29T14:56:52.430 回答
0

更新：

在与 Roland 进行了长时间的讨论（见评论）后，他确定这可能是一个回归错误并提交了一份错误报告。

由于问题是“为什么bins参数未知？”，我花了相当多的时间研究它，我会回答它。

您的示例显然来自 2013 年 10 月的此链接，其中使用了参数。如何正确解释 ggplot 的 stat_density2d

然而，它从来都不是一个记录在案的参数，也不清楚它是否在那里被使用。它可能是传递给其他库（如 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 年 10 月原始链接中发布的情节有很大不同，我想说stat_density2d从那时起已经被广泛重写，或者可能MASS:kde2d（或另一个 MASS 例程），并且该bin参数不再被接受。也不清楚该参数是否做过任何事情（阅读链接）。

但是，您可以更改参数n并且h- 从stat_density2d角度来看也没有记录（据我所知）。

于 2015-12-22T13:07:23.697 回答