ggplot2灰度计划：澄清变点图的建议

Question

我可以使用ggplot2绘制多个同时进行变更点和制度的时间序列，并且我可以使用颜色来使制度清晰（使用geom_rect绘制不同颜色的不同部分）。我需要制作一个情节，清楚地表明政权在没有使用颜色的地方。通过三种制度，可以区分使用白色，黑色和灰色进行着色的制度，但如果存在三种以上的制度，则很难区分它们。

我举了一个我可以使用颜色制作的情节的例子，如果有人可以推荐一个不使用颜色传达相同信息的情节，我将非常感激。

library(ggplot2)
library(scales)
# generate 3 time series and store them in a data frame
generate_cp_ts <- function(tau, params) {
    ts(c(arima.sim(model = list(ar = 0.2), n = tau[1], rand.gen = function(n) params[1] * rnorm(n)), arima.sim(model = list(ar = 0.2), n = tau[2] - tau[1], rand.gen = function(n) params[2] * rnorm(n)), arima.sim(model = list(ar = 0.2), n = tau[3] - tau[2], rand.gen = function(n) params[3] * rnorm(n)), arima.sim(model = list(ar = 0.2), n = tau[4] - tau[3], rand.gen = function(n) params[4] * rnorm(n))))
}
tau <- 100 * (1:4)
ts1 <- generate_cp_ts(tau, c(1.7, 0.3, 1.7, 1.7))
ts2 <- generate_cp_ts(tau, c(0.3, 2, 0.3, 0.9))
ts3 <- generate_cp_ts(tau, c(2, 2, 0.1, 0.7))
tsframe <- data.frame(ts = c(ts1, ts2, ts3), ts_level = factor(paste("Time Series", rep(1:3, each = 400))), time = rep(1:400, 3))
# Work out which colors are needed to color the plot and store in a data frame
CPs <- c(0, tau)
colour.frame <- data.frame(regime.from = rep(CPs[-length(CPs)], each = 3), regime.to = rep(CPs[-1], each = 3), ts_level = factor(paste("Time Series", rep(c(1:3), length(CPs) - 1))), regime = factor(c(0,0,0, 1,1,0, 0,0,1, 0,2,2) + 1))
# Plotting
qplot(x = time, y = ts, data = tsframe, facets = ts_level ~ ., alpha = I(1), geom = "line", ylab = "Time Series", xlab = "Time") +
    geom_rect(aes(NULL, NULL, xmin = regime.from, xmax = regime.to, fill = regime), ymin = -6, ymax = 6, data = colour.frame) +
    scale_fill_manual(values = alpha(c("blue", "red", "green"), 0.2))

Plot generated by the above code

Answer 1

创建colour.frame后，您可以插入以下代码：

tdf <- colour.frame

tdf$xval <- (tdf$regime.from + tdf$regime.to)/2
tdf$yval <- max(tsframe$ts) * 0.8 # if 0.8 is higher (0.9) then the text is set higher 

ggplot(tsframe, aes(x = time, y = ts)) + 
  geom_line() + 
  facet_grid(ts_level~.) + 
  geom_vline(xintercept = CPs) + # maybe play around with linetype
  geom_text(aes(x = xval, y = yval, label = regime), data = tdf)

给出了这个情节：