r ggplot 错误：美学长度匹配颜色和形状

Question

我正在尝试开发一个使用 ggplot 的函数，但我遇到了 aes() 函数和我将传递给 geom_point() 的参数的问题，特别是 color 和 shape。

这是一个非常简化的代码版本，它开始出现问题。了解代码试图实现的目标的一些背景可能会有所帮助。假设我们有带有变量 Dose 和 Response 的数据框 x，并且有一个名为 OutlierDetect(Dose,Response) 的函数返回 x 的索引，这似乎是某个拟合曲线的异常值。目标是将这些索引绘制成与其余数据不同的形状。

a <- ggplot(NULL, aes(x = Dose, y = Response)
shape.vec <- rep(19, nrow(x))
out.index <- OutlierDetect(Dose=x$Dose,Response=x$Response)
shape.vec[out.index] <- 17
a <- a + geom_point(data = x, size = 5, alpha = 0.8, aes(color = "group1"), shape = shape.vec)

我想避免在 x 中放置一个因子列这一看似显而易见的答案，因为如果可能的话，我不想修改 x。并且代码是这样写的，所以它足够灵活，可以添加一组新的数据 y，以便将它们绘制在同一个图上，并且 ggplot 会自动生成颜色图例。

我收到错误Error: Aesthetics must be either length 1 or the same as the data (1): shape, size, alpha。在我自己的研究中，制作图例时似乎是 ggplot 中的一个错误。出于某种原因，这两个参数需要具有相同的长度才能正常工作（参见下面的代码块）。

供参考，此代码块工作正常。

a <- ggplot(NULL, aes(x = Dose, y = Response)
shape.vec <- rep(19, nrow(x))
color.vec <- rep("blue", nrow(x)) #Need to manually specify new colors/color pallete which is more user input
out.index <- OutlierDetect(Dose=x$Dose,Response=x$Response)
shape.vec[out.index] <- 17
a <- a + geom_point(data = x, size = 5, alpha = 0.8, color = color.vec, shape = shape.vec)

这个块也是如此

a <- ggplot(NULL, aes(x = Dose, y = Response)
a <- a + geom_point(data = x, size = 5, alpha = 0.8, aes(color = "group1"), shape = 19) 
# can add more data frames y, z, ect, so long
# as the aes color parameter has a different name

总之，我想要的是这两个代码块的混合。有什么建议么？

可重现的示例代码

x <- data.frame(Dose = c(1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10),
Response = seq(20))
a <- ggplot(NULL, aes(x = Dose, y = Response))
shape.vec <- rep(19, nrow(x))
out.index <- nrow(x) # lets say the OutlierDetect function always says the last index is an outlier
shape.vec[out.index] <- 17
a <- a + geom_point(data = x, size = 5, alpha = 0.8, aes(color = "group1"),     shape = shape.vec)

Answer 1

您的代码有nrow(x$Dose)，当您尝试从一维向量调用行号时它不起作用。您可以将其更改为nrow(x) 或length(x$Dose)。

a <- ggplot(NULL, aes(x = Dose, y = Response))
shape.vec <- rep(19, nrow(x))
out.index <- nrow(x) # lets say the OutlierDetect function always says the last index is an outlier
shape.vec[out.index] <- 17
a <- a + geom_point(data = x, size = 5, alpha = 0.8, aes(color = "group1"),shape = shape.vec)