R + ggplot：有事件的时间序列

虽然我很喜欢@JD Long的答案，但我会在R/ggplot2中放一个.

该方法是创建第二个事件数据集，并使用该数据集确定位置.从@Angelo的经历开始:

library(ggplot2)
data(presidential)
data(economics)

取出事件(总统)数据，并对其进行转换.将baseline和offset计算为其将绘制的经济数据的分数.将底部(ymin)设置为基线.这就是棘手的部分.如果标签靠得太近，我们需要能够错开标签.因此，确定相邻标签之间的间距(假设事件已排序).如果小于某个数量(我 Select 了大约4年的数据)，那么请注意，该标签需要更高.但是它必须比后面的高，所以使用rle来获得TRUE的长度(也就是说，必须更高)，并使用它来计算偏移向量(TRUE的每个字符串必须从其长度倒计时到2，FALSE正好在偏移量1处).使用该值确定杆(ymax)的顶部.

events <- presidential[-(1:3),]
baseline = min(economics$unemploy)
delta = 0.05 * diff(range(economics$unemploy))
events$ymin = baseline
events$timelapse = c(diff(events$start),Inf)
events$bump = events$timelapse < 4*370 # ~4 years
offsets <- rle(events$bump)
events$offset <- unlist(mapply(function(l,v) {if(v){(l:1)+1}else{rep(1,l)}}, l=offsets$lengths, v=offsets$values, USE.NAMES=FALSE))
events$ymax <- events$ymin + events$offset * delta

把这些放在一起形成一个情节:

ggplot() +
    geom_line(mapping=aes(x=date, y=unemploy), data=economics , size=3, alpha=0.5) +
    geom_segment(data = events, mapping=aes(x=start, y=ymin, xend=start, yend=ymax)) +
    geom_point(data = events, mapping=aes(x=start,y=ymax), size=3) +
    geom_text(data = events, mapping=aes(x=start, y=ymax, label=name), hjust=-0.1, vjust=0.1, size=6) +
    scale_x_date("time") +  
    scale_y_continuous(name="unemployed \[1000's\]")

你可以 Select 分面，但不同的尺度很难做到.另一种方法是合成两个图.必须进行一些额外的调整，以确保绘图具有相同的x范围，使标签都适合下部绘图，并消除上部绘图中的x轴.

xrange = range(c(economics$date, events$start))

p1 <- ggplot(data=economics, mapping=aes(x=date, y=unemploy)) +
    geom_line(size=3, alpha=0.5) +
    scale_x_date("", limits=xrange) +  
    scale_y_continuous(name="unemployed [1000's]") +
    opts(axis.text.x = theme_blank(), axis.title.x = theme_blank())

ylims <- c(0, (max(events$offset)+1)*delta) + baseline
p2 <- ggplot(data = events, mapping=aes(x=start)) +
    geom_segment(mapping=aes(y=ymin, xend=start, yend=ymax)) +
    geom_point(mapping=aes(y=ymax), size=3) +
    geom_text(mapping=aes(y=ymax, label=name), hjust=-0.1, vjust=0.1, size=6) +
    scale_x_date("time", limits=xrange) +
    scale_y_continuous("", breaks=NA, limits=ylims)

#install.packages("ggExtra", repos="http://R-Forge.R-project.org")
library(ggExtra)

align.plots(p1, p2, heights=c(3,1))