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

Question

I do apologise in advance but I'm not an R expert (at all). I used R rarely and the last time, I had to convert my data as tibble in order to plot them correctly, and now I did the same but I believe it would be better to convert them into POSIXct (but I don't know how to do it:confused: ) I have this 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))

Now... So far I managed to get them into the plot style I need, using the following code:

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()

Which gives this:

I know, it's absolutely horrible, but it was the only way I remember. What I do really need is to add a shaded rectangle in every plot (at the same date); something like this, for instance:

Could you kindly help me with this matter, please? I don't mind changing the data format to POSIXct (if required for ggplot), anything really, as long as it looks like this.

Thanks a lot. PS I don't need the data with the time (ie,: 01:01:01 can be removed if needed!)

Answer 1

Add the line geom_rect(aes(xmin=as.POSIXct("2006/01/01"), xmax=as.POSIXct("2006/06/01"), ymin=0, ymax=1), alpha=.01) to your ggplot code.