我有一堆数据 - 三个时间序列(模型组均值),按组着色,标准偏差由geom_ribbon表示。默认情况下,他们有一个很好的传奇。我还有一个时间序列的观察,我想叠加在情节上(没有geom_ribbon),像这样:
df <- data.frame(year=1991:2010, group=c(rep('group1',20), rep('group2',20), rep('group3',20)), mean=c(cumsum(abs(rnorm(20))),cumsum(abs(rnorm(20))),cumsum(abs(rnorm(20)))),sd=3+rnorm(60))
obs_df <- data.frame(year=1991:2010, value=cumsum(abs(rnorm(20))))
ggplot(df, aes(x=year, y=mean)) + geom_line(aes(colour=group)) + geom_ribbon(aes(ymax=mean+sd, ymin=mean-sd, fill=group), alpha = 0.2) +geom_line(data=obs_df, aes(x=year, y=value))
但观察确实出现在传说中,因为它没有着色(我希望它是黑色的)。如何将遮挡添加到图例?
答案 0 :(得分:6)
首先,创建df
和obs_df
的组合数据框:
dat <- rbind(df, data.frame(year = obs_df$year,
group = "obs", mean = obs_df$value, sd = 0))
简介:
ggplot(dat, aes(x=year, y=mean)) +
geom_line(aes(colour=group)) +
geom_ribbon(aes(ymax=mean+sd, ymin=mean-sd, fill=group), alpha = 0.2) +
scale_colour_manual(values = c("red", "green", "blue", "black")) +
scale_fill_manual(values = c("red", "green", "blue", NA))
答案 1 :(得分:0)
我猜你在构造'obs_df'时犯了错误。如果使用year = 1991:2010
创建它,则在其余数据的上下文中更有意义,它会为您提供您希望ggplot
调用保持不变的情节。