使用公共阴影矩形 (geom_rect) 创建多面 plot

Question

我确实提前道歉，但我不是 R 专家（根本）。 我很少使用 R ，最后一次，我必须将我的数据转换为 tibble 以便正确地 plot 它们，现在我做了同样的事情，但我相信将它们转换为 POSIXct 会更好（但我不知道如何这样做：困惑：）我有这个 csv：

structure(list(Date = c("01/02/2003 01:01:01", "01/03/2003 01:01:01", 
"01/04/2003 01:01:01", "01/05/2003 01:01:01", "01/06/2003 01:01:01", 
"01/07/2003 01:01:01", "01/08/2003 01:01:01", "01/09/2003 01:01:01", 
"01/10/2003 01:01:01", "01/11/2003 01:01:01", "01/12/2003 01:01:01", 
"01/01/2004 01:01:01", "01/02/2004 01:01:01", "01/03/2004 01:01:01", 
"01/04/2004 01:01:01", "01/05/2004 01:01:01", "01/06/2004 01:01:01", 
"01/07/2004 01:01:01", "01/08/2004 01:01:01", "01/09/2004 01:01:01", 
"01/10/2004 01:01:01", "01/11/2004 01:01:01", "01/12/2004 01:01:01", 
"01/01/2005 01:01:01", "01/02/2005 01:01:01", "01/03/2005 01:01:01", 
"01/04/2005 01:01:01", "01/05/2005 01:01:01", "01/06/2005 01:01:01", 
"01/07/2005 01:01:01", "01/08/2005 01:01:01", "01/09/2005 01:01:01", 
"01/10/2005 01:01:01", "01/11/2005 01:01:01", "01/12/2005 01:01:01", 
"01/01/2006 01:01:01", "01/02/2006 01:01:01", "01/03/2006 01:01:01", 
"01/04/2006 01:01:01", "01/05/2006 01:01:01", "01/06/2006 01:01:01", 
"01/07/2006 01:01:01", "01/08/2006 01:01:01", "01/09/2006 01:01:01", 
"01/10/2006 01:01:01", "01/11/2006 01:01:01"), Pz1 = c(0.53, 
0.25, 0.3, 0.51, 0.23, 0.52, 0.48, 0.36, 0.56, 0.27, 0.44, 0.28, 
0.79, 0.15, 0.73, 0.44, 0.5, 0.26, 0.1, 0.26, 0.69, 0.38, 0.51, 
0.39, 0.42, 0.29, 0.68, 0.62, 0.5, 0.06, 0.29, 0.13, 0.6, 0.21, 
0.34, 0.17, 0.39, 0.21, 0.89, 0.19, 0.44, 0.53, 0.55, 0.89, 0.55, 
0.65), Pz2 = c(0.62, 0.99, 0.87, 0.77, 0.51, 0.66, 0.4, 0.68, 
0.87, 0.13, 0.81, 0.29, 0.11, 0.11, 0.23, 0.71, 0.85, 0.05, 0.78, 
0.32, 0.16, 0.54, 0.65, 0.09, 0.97, 0.81, 0.49, 0.36, 0.37, 0.78, 
0.04, 0.67, 0.91, 0.12, 0.34, 0.3, 0.71, 0.04, 0.73, 0.33, 0.59, 
0.23, 0.82, 0.04, 0.04, 0.82), Pz3 = c(0.2, 0.76, 0.04, 0.39, 
0.58, 0.49, 0.44, 0.12, 0.16, 0.12, 0.95, 0.95, 0.08, 0.68, 0.57, 
0.49, 0.58, 0.46, 0.39, 0.51, 0.69, 0.09, 0.68, 0.18, 0.3, 0.75, 
0.76, 0.85, 0.17, 0.6, 0.45, 0.26, 0.65, 0.07, 0.7, 0.71, 0.47, 
0.79, 0.58, 0.08, 0.37, 0.86, 0.23, 0.31, 0.06, 0.1), Pz4 = c(0.21, 
0.65, 0.67, 0.45, 0.32, 0.79, 0.94, 0.78, 0.73, 0.83, 0.79, 0.46, 
0.07, 0.84, 0.25, 0.27, 0.77, 0.37, 0.16, 0.67, 0.88, 0.67, 0.87, 
0.95, 0.63, 0.61, 0.21, 0.21, 0.4, 0.74, 0.62, 0.22, 0.08, 0.67, 
0.2, 0.18, 0.83, 0.3, 0.15, 0.7, 0.5, 0.43, 0.81, 0.17, 0.31, 
0.66), Pz5 = c(0.57, 0.78, 1, 0.87, 0.88, 0.5, 0.24, 0.71, 0.11, 
0.4, 0.08, 0.2, 0.67, 0.41, 0.28, 0.45, 0.6, 0.18, 0.27, 0.02, 
0.96, 0.48, 0.95, 0.01, 0.8, 0.07, 0.34, 0.09, 0.19, 0.59, 0.34, 
0.66, 0.48, 0.86, 0.97, 0.76, 0.93, 0.21, 0.5, 0.93, 0.41, 0.33, 
0.32, 0.12, 0.42, 0.94), Pz6 = c(0.42, 0.34, 0.34, 0.73, 0.7, 
0.67, 0.09, 0.45, 0.55, 0.88, 0.05, 0.15, 0.85, 0.02, 0.42, 0.14, 
0.68, 0.71, 0.57, 0.14, 0.85, 0.81, 0.2, 0.97, 0.42, 0.59, 0.23, 
0.39, 0.5, 0.87, 0.37, 0.63, 0.7, 0.3, 0.33, 0.29, 0.9, 0.75, 
0.38, 0.17, 0.87, 0.45, 0.79, 0.74, 0.21, 0.05)), class = "data.frame", row.names = c(NA, 
-46L))

现在......到目前为止，我设法使用以下代码将它们放入我需要的 plot 样式中：

a=read.csv("C:/Users/simon/Desktop/4.csv")
b=a
c=b %>% 
    mutate(Date = dmy_hms(Date)) %>% 
    arrange(Date) %>% 
    as_tbl_time(index = Date)
df=c
df_melt = melt(df, id= "Date")
d=ggplot(df_melt, aes(x = Date, y = value)) + 
    geom_line() + 
    facet_wrap(~ variable, scales = 'free_y', ncol = 2)
d+theme_bw()

这给出了这个：

我知道，这绝对是可怕的，但这是我记得的唯一方式。 我真正需要的是在每个 plot （在同一日期）中添加一个阴影矩形； 像这样的东西，例如：

你能帮我解决这个问题吗？ 我不介意将数据格式更改为 POSIXct（如果 ggplot 需要），只要看起来像这样。

非常感谢。 PS 我不需要带时间的数据（即：01:01:01 如果需要可以去掉！）

Answer 1

添加行geom_rect(aes(xmin=as.POSIXct("2006/01/01"), xmax=as.POSIXct("2006/06/01"), ymin=0, ymax=1), alpha=.01)到您的ggplot代码。