So here I am not able to replicate the exact plot as shown in the data but have been able to do what has been said in the text. Reason exact plot has not been replicated is that the data does not support it. Please find the code below
times = seq(as.Date("2017-01-01"), as.Date("2017-12-31"), "1 day") %>% format("%b-%d")
df1 = data.frame(date = times, id = rep(1,365),value = rnorm(n=365))
df2 = data.frame(date = times, id = rep(2,365),value = rnorm(n=365))
df3 = data.frame(date = times, id = rep(3,365),value = rnorm(n=365))
df4 = data.frame(date = times, id = rep(4,365),value = rnorm(n=365))
df5 = data.frame(date = times, id = rep(5,365),value = rnorm(n=365))
df6 = data.frame(date = times, id = rep(6,365),value = rnorm(n=365))
DF = do.call('rbind',list(df1,df2,df3,df4,df5,df6))
test <- DF %>% group_by(date) %>% summarize(Max= max(value),Min=min(value))
ggplot(test) + geom_ribbon(aes(x = 1:nrow(test),ymin = Min,ymax = Max),fill="darkorchid") +
scale_x_continuous(breaks = seq(1,nrow(test),5),labels = test$date[seq(1,nrow(test),5)],name = NULL) +
theme(axis.text.x = element_text(angle = 90)) +
labs(title = "Ribbon plot of desire")