示例数据
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()
这里出了什么问题?
答案 0 :(得分:2)
问题在于,当你在审美中使用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()
答案 1 :(得分:0)
你的cumsum不是每年这里
使用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()