ggplot2 中的累积图

Cumulative plot in ggplot2

示例数据

dat <- data.frame(year = as.factor(rep(c(2012:2015),each = 6)),id.2wk = rep(c(18,19,20,21,22,23),times = 4), 
                    value = c(1.8,15.6,32.9,27.5,19.6,2.6,1,8,42,35,11,3,2,7,12,47,26,7,2,13,24,46,12,4))

with(dat, plot(id.2wk[year == 2012], cumsum(value[year == 2012]), type = "b"))
with(dat, points(id.2wk[year == 2013], cumsum(value[year == 2013]), type = "b"))
with(dat, points(id.2wk[year == 2014], cumsum(value[year == 2014]), type = "b"))
with(dat, points(id.2wk[year == 2015], cumsum(value[year == 2015]), type = "b"))

我想使用 ggplot2 创建相同的绘图。我这样做了:

  ggplot(dat, aes(x = id.2wk, y = cumsum(value), colour = factor(year))) + 
  geom_line(size = 1)+
  geom_point() 

这里出了什么问题?

这里你的月经不是每年的

有了 data.table 你可以做到

library(data.table)
dat <- setDT(dat)
dat[,cumulsum :=cumsum(value), by = year]
  ggplot(, aes(x = id.2wk, y = cumulsum, colour = factor(year))) + 
  geom_line(size = 1)+
  geom_point() 

问题是,当您在美学中使用 cumsum() 时,它适用于所有值,而不仅仅是特定年份内的值。

与其用ggplot做变换,不如先用dplyr做变换,然后绘制结果会更安全。例如

ggplot(dat %>% group_by(year) %>% mutate(cv=cumsum(value)), 
      aes(x = id.2wk, y = cv, colour = factor(year))) + 
  geom_line(size = 1)+
  geom_point()