p1 <- GHI3 %>% arrange(max) %>% mutate(name = factor(name, levels = name)) %>% 
    ggplot(aes(x = name)) + geom_segment(aes(x = name, xend = name, y = min, 
    yend = max), color = "#c5c5c5", size = 2) + geom_point_interactive(aes(y = max), 
    color = "#1b7fbd", size = 2) + geom_point_interactive(aes(y = min), color = "#ee4540", 
    size = 2) + coord_flip() + bbc_style()


finalise_plot(p1, source_name = "Source: GaStat", save_filepath = "solar1.png", 
    logo_image_path = "FA_LOGO.png")

GHI_long <- GHI3 %>% gather(year, GHI, -name) %>% filter(!year %in% c("max", 
    "min"))
GHI_long$year <- as.numeric(GHI_long$year)

p2 <- ggplot(GHI_long, aes(x = year)) + stat_smooth(aes(y = GHI), se = FALSE, 
    color = "#c5c5c5") + geom_point_interactive(aes(y = GHI, color = as.factor(year), 
    tooltip = name), position = position_jitter(width = 0.3, height = 0), size = 2) + 
    # geom_jitter(aes(y=GHI, color=as.factor(year)), width = 0.2)+
# geom_boxplot(aes(y = GHI, color=as.factor(year)), alpha = 0.3)
geom_flat_violin(aes(y = GHI, fill = as.factor(year)), alpha = 0.3) + scale_fill_manual(values = c("#88A9D7", 
    "#a2a8d3", "#38598b", "#113f67", "#263859")) + scale_color_manual(values = c("#88A9D7", 
    "#a2a8d3", "#38598b", "#113f67", "#263859")) + 
stat_summary(aes(y = GHI), fun.data = mean_sdl, mult = 1, geom = "pointrange", 
    color = "#3B60B0", size = 1) + bbc_style() + theme(legend.position = "none") + 
    labs(subtitle = "GHI(Wh/m2/day) change over years")



finalise_plot(p2, source_name = "Source: GaStat", save_filepath = "solar2.png", 
    logo_image_path = "FA_LOGO.png")

station2[11, 2] <- "K.A.CARE location in Riyadh"
station2[41, 2] <- "Prince Sattam Bin Abdulaziz University (Kharj)"
station2[9, 2] <- "Hafar Al-Batin Technical Institute"
station2[17, 2] <- "King Faisal University (Ahsa)"
station2[15, 2] <- "King Abdullah University"
station2[39, 2] <- "Taibah University in Madenah"



GHI_final <- full_join(station2, GHI3, "name")

p3 <- ggplot(GHI_final, aes(x = elev, y = median)) + geom_point(color = "#c5c5c5") + 
    geom_point(data = GHI_final[c(31, 26, 36, 40), ], color = "blue") + stat_smooth(method = "lm") + 
    geom_text(x = 100, y = 6500, label = paste("R=", round(cor(GHI_final$elev, 
        GHI_final$median), 2))) + geom_label_repel(data = GHI_final[c(31, 26, 
    36, 40), ], color = "#263859", aes(elev, median, label = name)) + bbc_style()

finalise_plot(p3, source_name = "Source: GaStat", save_filepath = "solar3.png", 
    logo_image_path = "FA_LOGO.png")