-
Notifications
You must be signed in to change notification settings - Fork 1
/
whitepaperplots.R
94 lines (64 loc) · 3.85 KB
/
whitepaperplots.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
#Graphs for the short-term white paper
library(tidyverse)
library(lubridate)
library(DroughtData)
load("data/NutrientsChlorophyll.RData")
yrs = read_csv("data/yearassignments.csv") %>%
rename(YearAdj = Year) %>%
select(YearAdj, Whitepaper)
#nitrate graph
Nitrate = left_join(Nitrate, yrs) %>%
mutate(Whitepaper = factor(Whitepaper, levels = c("Critical", "Dry", "Below Normal", "Above Normal", "Wet", "2020", "2021", "2022"),
ordered = T))
ggplot(Nitrate, aes(x = Whitepaper, y = DissNitrateNitrite, fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() + theme_bw() + ylab("Dissolved Nitrate+Nitrite (mg/L)")+xlab(NULL)
ggsave("plots/whitepaper/NitrateWhitepaper.tiff", device = "tiff", width = 8, height = 6)
#Phosphorus graph
Phos = left_join(Phos, yrs) %>%
mutate(Whitepaper = factor(Whitepaper, levels = c("Critical", "Dry", "Below Normal", "Above Normal", "Wet", "2020", "2021", "2022"), ordered = T))
ggplot(Phos, aes(x = Whitepaper, y = DissOrthophos, fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() + theme_bw() + ylab("Dissolved Orthophosphate (mg/L)")+xlab(NULL)
ggsave("plots/whitepaper/PHosWhitepaper.tiff", device = "tiff", width = 8, height = 6)
#Chlorophyll
ChlaA = left_join(ChlaA, yrs) %>%
mutate(Whitepaper = factor(Whitepaper, levels = c("Critical", "Dry", "Below Normal", "Above Normal", "Wet", "2020", "2021", "2022"), ordered = T))
ggplot(ChlaA, aes(x = Whitepaper, y = log(Chlorophyll), fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() + theme_bw() + ylab("log Chlorophyll-a (ug/L)")+xlab(NULL)
ggplot(ChlaA, aes(x = Whitepaper, y = Chlorophyll, fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() + theme_bw() + ylab("Chlorophyll-a (ug/L)")+xlab(NULL) +
facet_wrap(~Region, nrow = 4, strip.position = "right", scales = "free_y")+
scale_x_discrete(labels = c("Critical", "Dry", "Below\nNormal", "Above\nNormal", "Wet", "2020", "2021"))+
theme(legend.position ="none")
ggsave("plots/whitepaper/ChlWhitepaper.tiff", device = "tiff", width = 5, height = 6)
############################################
# Temperature
yrs = read_csv("data/yearassignments.csv") %>%
rename(YearAdj = Year) %>%
dplyr::select(YearAdj, Whitepaper)
wq = lt_avg_wq %>%
left_join(yrs)%>%
mutate(Whitepaper = factor(Whitepaper, levels = c("Critical", "Dry", "Below Normal", "Above Normal", "Wet", "2020", "2021", "2022"),
ordered = T),
Season = factor(Season, levels = c("Winter", "Spring", "Summer", "Fall")))
#secchi depth
ggplot(wq, aes(x = Whitepaper, y = Secchi, fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() + theme_bw() +
ylab("Secchi Depth (cm)")+xlab(NULL)+
scale_x_discrete(labels = c("Critical", "Dry", "Below\nNormal", "Above\nNormal", "Wet", "2020", "2021"))+
theme(legend.position ="none")
ggsave("plots/whitepaper/SecchiWhitepaper.tiff", device = "tiff", width = 5, height = 4)
#Temperature
ggplot(wq, aes(x = Whitepaper, y = Temperature, fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() +
facet_wrap(~Season, nrow = 4, scales = "free_y", strip.position = "right")+
theme_bw() +
scale_x_discrete(labels = c("Critical", "Dry", "Below\nNormal", "Above\nNormal", "Wet", "2020", "2021"))+
ylab("Water Temperature (C)")+xlab(NULL)+
theme(legend.position ="none")
ggsave("plots/whitepaper/TempWhitepaper.tiff", device = "tiff", width = 5, height = 6)
#Salinity
ggplot(wq, aes(x = Whitepaper, y = Salinity, fill = YearType)) + geom_boxplot()+
drt_color_pal_yrtype() + facet_wrap(~Region, nrow = 5, scales = "free_y", strip.position = "right")+
theme_bw() + scale_x_discrete(labels = c("Critical", "Dry", "Below\nNormal", "Above\nNormal", "Wet", "2020", "2021"))+
ylab("Salinity (PSU)")+xlab(NULL)+ theme(legend.position ="none")
ggsave("plots/whitepaper/SalWhitepaper.tiff", device = "tiff", width = 5, height = 7)