0

我有两个 geom_line 时间序列,每个图上有一些 geom_text 和 geom_points。

两个地块之一还具有 geom_vline。我想知道是否可以合并两个,但我没有得到解决方案。

这是两个地块:

require(zoo)
require(ggplot2)

set.seed(10)

# plot 1:
tmp1 <- xts(cumsum(rnorm(5000,1,10)), Sys.Date()-5000:1)
data.tmp1 = data.frame(date=as.Date(index(tmp1)),
                      value=drop(coredata(tmp1)))
data.tmp1.year.end = data.frame(date=as.Date(index(tmp1[endpoints(tmp1, "years", 1)])),
                               value= drop(coredata(tmp1[endpoints(tmp1, "years", 1)])))

plot1 = 
  ggplot(data.tmp1, aes(x=date, y=value)) +
  geom_line(aes(y=value), size=1) +
  geom_point(data=data.tmp1.year.end, col="red") +
  geom_text(data=data.tmp1.year.end, label=data.tmp1.year.end$value, vjust=0, hjust=1) 

# plot 2:
tmp2 <- xts(cumsum(rnorm(5000,1,100)), Sys.Date()-5000:1)
data.tmp2 = data.frame(date=as.Date(index(tmp2)),
                       value=drop(coredata(tmp2)))
data.tmp2.year.end = data.frame(date=as.Date(index(tmp2[endpoints(tmp2, "years", 1)])),
                                value= drop(coredata(tmp2[endpoints(tmp2, "years", 1)])))
tmp2.date =as.Date(c("2008-01-01"))

plot2 = 
  ggplot(data.tmp2, aes(x=date, y=value)) +
  geom_line(aes(y=value), size=1) +
  geom_point(data=data.tmp2.year.end, col="red") +
  geom_vline(xintercept=as.numeric(tmp2.date), linetype="dotted") +
  geom_text(data=data.tmp2.year.end, label=data.tmp2.year.end$value, vjust=0, hjust=1)

现在的目标是 plot1 和 plot2 共享一个 x 轴,并且各个图的所有特征都保存在相应的图中。

结果应如下所示:

在此处输入图像描述

4

1 回答 1

1

您可以尝试将每日数据集和年终数据集组合成单个数据框,然后使用 ggplot 的分面显示在单个日期轴上。代码可能如下所示:

data.tmp1 <- cbind(data.tmp1, data_name="tmp1")
data.tmp1.year.end <- cbind(data.tmp1.year.end, data_name="tmp1")
data.tmp2 <- cbind(data.tmp2, data_name="tmp2")
data.tmp2.year.end <- cbind(data.tmp2.year.end, data_name="tmp2")

data.tmp <- rbind(data.tmp1,data.tmp2)
data.tmp.year.end <- rbind(data.tmp1.year.end, data.tmp2.year.end)

ggplot(data.tmp, aes(x=date, y=value)) +
  geom_line(aes(y=value), size=1) +
  geom_point(data=data.tmp.year.end, col="red") +
  geom_text(data=data.tmp.year.end, aes(label=data.tmp.year.end$value), vjust=0, hjust=1) +
  geom_vline(xintercept=as.numeric(tmp2.date), linetype="dotted") +
  facet_grid( data_name ~ . , scales="free_y")

这给出了图表

在此处输入图像描述

于 2015-03-19T14:02:54.313 回答