如何在 R 中使用 Plot 实现 plot 多个时间序列（带垂直面的公共 x 轴）？

Question

I am trying to plot a Vertical stack plot for 4 time series in R. The following codes can show the figure, but the 'lty' and 'col' parameters are not working correctly. Meanwhile, although "main=NULL" is used, the title still shows. Is there any way to fix this?

In addition, is there any way to set the "xlim" for each small plot? Thank you.

data(UKgas)
ts1 = UKgas
ts2 = UKgas+100
ts3 = UKgas+500
ts4 = UKgas+1000

plot(ts.union(TS1 = ts1, TS2 = ts2,TS3 = ts3, TS4 = ts4), lty = 1:4, col = 1:4, main = NULL)

Answer 1

Plotting a ts.union object calls the method stats:::plot.ts , and from looking at the source code, this will only take a single color I'm afraid. You can get the result you want using a little data wrangling and ggplot2 , with the added benefit that the plot is fully customizable.

library(ggplot2)
library(tidyr)

ts_s  <- ts.union(TS1 = ts1, TS2 = ts2,TS3 = ts3, TS4 = ts4)
times <- attr(ts_s, "tsp")
df    <- cbind(as.data.frame(ts_s), year = seq(times[1], times[2], 1/times[3]))

ggplot(pivot_longer(df, 1:4), aes(year, value, colour = name)) +
  geom_line(aes(linetype = name)) +
  scale_x_continuous(breaks = seq(1960, 1985, 5)) +
  facet_grid(name~., switch = "y", scales = "free_y") +
  labs(title = NULL) +
  theme_classic() +
  theme(panel.background = element_rect(color = "black"),
        strip.placement = "outside",
        strip.background = element_blank(),
        panel.spacing.y = unit(0, "npc"),
        strip.text = element_text(size = 16, face = 2),
        legend.position = "none")

^{Created on 2022-03-13 by the reprex package (v2.0.1)}