R语言 分类线图上的自定义ggplot2阴影错误区域

3ks5zfa0  于 2023-11-14  发布在  其他
关注(0)|答案(1)|浏览(87)

我试图绘制一条线,用黄土平滑,但我试图弄清楚如何包括由现有变量定义的阴影误差区域,但也平滑。
此代码创建示例数据:

set.seed(12345)
data <- cbind(rep("A", 100), rnorm(100, 0, 1))
data <- rbind(data, cbind(rep("B", 100), rnorm(100, 5, 1)))
data <- rbind(data, cbind(rep("C", 100), rnorm(100, 10, 1)))
data <- rbind(data, cbind(rep("D", 100), rnorm(100, 15, 1)))
data <- cbind(rep(1:100, 4), data)
data <- data.frame(data)
names(data) <- c("num", "category", "value")
data$num <- as.numeric(data$num)
data$value <- as.numeric(data$value)
data$upper <- data$value+0.20
data$lower <- data$value-0.30

字符串
将下面的数据绘制出来,这就是我得到的:

ggplot(data, aes(x=num, y=value, colour=category)) +
  stat_smooth(method="loess", se=F)


的数据
我想要的是一个类似于下面的图,除了阴影区域的上界和下界由生成的数据中的“上”和“下”变量的平滑线限定。

任何帮助将不胜感激。

f3temu5u

f3temu5u1#

这里有一种添加upperlower的平滑版本的方法。我们将upperlower的LOESS预测添加到 Dataframe 中,然后使用geom_ribbon绘制这些预测。如果这一切都可以在对ggplot的调用中完成,那将更加优雅。这可能通过向stat_summary提供专用函数来实现,希望其他人也会用这种方法来回答。

# Expand the scale of the upper and lower values so that the difference
# is visible in the plot
data$upper = data$value + 10
data$lower = data$value - 10

# Order data by category and num
data = data[order(data$category, data$num),]

# Create LOESS predictions for the values of upper and lower 
# and add them to the data frame. I'm sure there's a better way to do this,
# but my attempts with dplyr and tapply both failed, so I've resorted to the clunky 
# method below.
data$upperLoess = unlist(lapply(LETTERS[1:4], 
                  function(x) predict(loess(data$upper[data$category==x] ~ 
                                                  data$num[data$category==x]))))
data$lowerLoess = unlist(lapply(LETTERS[1:4], 
                  function(x) predict(loess(data$lower[data$category==x] ~ 
                                                  data$num[data$category==x]))))

# Use geom_ribbon to add a prediction band bounded by the LOESS predictions for 
# upper and lower
ggplot(data, aes(num, value, colour=category, fill=category)) +
  geom_smooth(method="loess", se=FALSE) +
  geom_ribbon(aes(x=num, y=value, ymax=upperLoess, ymin=lowerLoess), 
              alpha=0.2)

字符串
结果是这样的:


的数据

相关问题