生成堆叠累积平滑频率分布图

Question

我的数据中记录了两种类型的出现：type_a 和 type_b 以及它们的出现年份。

这是生成我的数据示例的一种方法：

set.seed(1)
years <- 1991:2010
type_a_years <- 20
type_b_years <- 10
type_a <- round(runif(type_a_years, 0, 5))
type_b <- c(rep(0, type_a_years-type_b_years),round(runif(type_b_years, 5, 7)))

df <- data.frame(year = unlist(sapply(1:length(years), function(x) c(rep(years[x], type_a[x]),rep(years[x], type_b[x])))),
                 type = unlist(sapply(1:length(years), function(x) c(rep("type_a", type_a[x]),rep("type_b", type_b[x])))))

head(df)
  year   type
1 1991 type_a
2 1992 type_a
3 1992 type_a
4 1993 type_a
5 1993 type_a
6 1993 type_a

我想生成按年份堆叠在 type_a 事件之上的 type_b 事件的累积频率分布图，并且我希望分布显示为曲线而不是条形。

我猜这应该是一些操纵：

library(ggplot2)
ggplot(df, aes(year)) + stat_ecdf()

我将得到两条曲线并按类型堆叠，其中每种类型下的区域将填充不同的颜色。即type_a曲线与x轴之间的区域为一种颜色，type_b曲线与type_a曲线之间的区域为另一种颜色。

Answer 1

您可以在dplyr 和tidyr 中进行一些聚合，然后使用geom_area。

library(tidyr)
library(dplyr)
df1 <- df %>% group_by(type, year) %>%
              summarise(total = n()) %>%
              mutate(total = cumsum(total)) %>%
              ungroup %>%
              complete(type, year, fill = list(total = 0))

现在剧情：

library(ggplot2)
ggplot(df1, aes(x = year, y = total, fill = type)) + geom_area()