Skip to content
Snippets Groups Projects
Commit 4c5425fc authored by Ben Anderson's avatar Ben Anderson
Browse files

amende dplot labels for clarity

parent db4f5bba
No related branches found
No related tags found
No related merge requests found
...@@ -761,7 +761,7 @@ wcDT <- baselineWE_DR_DT[, .(nObs = .N, ...@@ -761,7 +761,7 @@ wcDT <- baselineWE_DR_DT[, .(nObs = .N,
] ]
makeComparisonPlot <- function(dt){ makeComparisonPlot <- function(dt){
myPlot <- ggplot(dt, aes(x = obsDate)) + myPlot <- ggplot(dt, aes(x = obsDate)) +
geom_line(aes(y = meanBaseline.madj, colour = "Baseline")) + geom_line(aes(y = meanBaseline.madj, colour = "Baseline (no WE)")) +
geom_line(aes(y = meanBaseline.madj.we, colour = "WE model")) + geom_line(aes(y = meanBaseline.madj.we, colour = "WE model")) +
geom_line(aes(y = meanBaseline.madj.we.dr, colour = "WE + Drought model")) + geom_line(aes(y = meanBaseline.madj.we.dr, colour = "WE + Drought model")) +
theme(legend.title = element_blank()) + theme(legend.title = element_blank()) +
...@@ -950,7 +950,7 @@ ggsave(paste0("plots_v2/sumDailyOverTimeByModel.pdf"), dpi = 400) ...@@ -950,7 +950,7 @@ ggsave(paste0("plots_v2/sumDailyOverTimeByModel.pdf"), dpi = 400)
Let's compare the baseline we with baseline we + dr results for each version of the model Let's compare the baseline we with baseline we + dr results for each version of the model
```{r selectData} ```{r selectData fpr plotTotalDT}
plotTotalDT <- baselineWE_DR_DT[, .( plotTotalDT <- baselineWE_DR_DT[, .(
meanBaselineWeSum = mean(sumDaily.baseline.madj.we, na.rm = TRUE), meanBaselineWeSum = mean(sumDaily.baseline.madj.we, na.rm = TRUE),
sdBaselineWeSum = sd(sumDaily.baseline.madj.we, na.rm = TRUE), sdBaselineWeSum = sd(sumDaily.baseline.madj.we, na.rm = TRUE),
...@@ -958,7 +958,7 @@ plotTotalDT <- baselineWE_DR_DT[, .( ...@@ -958,7 +958,7 @@ plotTotalDT <- baselineWE_DR_DT[, .(
sdBaselineDrSum = sd(sumDaily.baseline.madj.we.dr, na.rm = TRUE), sdBaselineDrSum = sd(sumDaily.baseline.madj.we.dr, na.rm = TRUE),
nObs = .N nObs = .N
) )
, by = .(obsDate, model, metered) , by = .(obsDate, model, metered, droughtStage = Colne)
] ]
# calculate CIs if we want them # calculate CIs if we want them
...@@ -966,6 +966,7 @@ plotTotalDT <- plotTotalDT[, yminWE := meanBaselineWeSum - qnorm(0.975)*(sdBasel ...@@ -966,6 +966,7 @@ plotTotalDT <- plotTotalDT[, yminWE := meanBaselineWeSum - qnorm(0.975)*(sdBasel
plotTotalDT <- plotTotalDT[, ymaxWE := meanBaselineWeSum + qnorm(0.975)*(sdBaselineWeSum/sqrt(nObs))] plotTotalDT <- plotTotalDT[, ymaxWE := meanBaselineWeSum + qnorm(0.975)*(sdBaselineWeSum/sqrt(nObs))]
plotTotalDT <- plotTotalDT[, yminDr := meanBaselineDrSum - qnorm(0.975)*(sdBaselineDrSum/sqrt(nObs))] plotTotalDT <- plotTotalDT[, yminDr := meanBaselineDrSum - qnorm(0.975)*(sdBaselineDrSum/sqrt(nObs))]
plotTotalDT <- plotTotalDT[, ymaxDr := meanBaselineDrSum + qnorm(0.975)*(sdBaselineDrSum/sqrt(nObs))] plotTotalDT <- plotTotalDT[, ymaxDr := meanBaselineDrSum + qnorm(0.975)*(sdBaselineDrSum/sqrt(nObs))]
``` ```
...@@ -1054,6 +1055,31 @@ myDrPlot + ...@@ -1054,6 +1055,31 @@ myDrPlot +
guides(linetype=guide_legend(ncol=1)) + guides(linetype=guide_legend(ncol=1)) +
ggsave(paste0("plots_v2/compareBaselineWEDroughtWE_v2_0.pdf"), height = 4, dpi = 480) ggsave(paste0("plots_v2/compareBaselineWEDroughtWE_v2_0.pdf"), height = 4, dpi = 480)
dt <- plotTotalDT[model == "v2_0"]
dt[, diff := meanBaselineDrSum - meanBaselineWeSum]
dt[, pc_diff := 100*(meanBaselineDrSum - meanBaselineWeSum)/meanBaselineWeSum]
t <- dt[, .(maxLitres_per_day_we = max(meanBaselineWeSum),
minLitres_per_day_we = min(meanBaselineWeSum)), keyby = .(droughtStage, metered)]
kable(t, caption = "Summary of Model v2 by drought stage", digits = 3)
min <- min(dt$pc_diff)
head(dt[pc_diff < -2.9, .(obsDate, droughtStage, metered, pc = round(pc_diff,2))])
twe <- dt[, .(maxL = max(meanBaselineWeSum),
minL = min(meanBaselineWeSum)), keyby = .(metered)]
twe[, Model := "water efficiency"]
twe[, pc_reduction := 100*(1 - (minL/maxL))]
tdr <- dt[, .(maxL = max(meanBaselineDrSum),
minL = min(meanBaselineDrSum)), keyby = .(metered)]
tdr[, pc_reduction := 100*(1 - (minL/maxL))]
tdr[, Model := "water efficiency + drought"]
t <- rbind(twe, tdr)
kable(t[, .(Model, metered, maxL, minL, pc_reduction)], digits = 2)
``` ```
...@@ -1086,6 +1112,26 @@ myDrPlot + ...@@ -1086,6 +1112,26 @@ myDrPlot +
guides(linetype=guide_legend(ncol=1)) + guides(linetype=guide_legend(ncol=1)) +
ggsave(paste0("plots_v2.1/compareBaselineWEDroughtWE_v2_1.pdf"), height = 4, dpi = 480) ggsave(paste0("plots_v2.1/compareBaselineWEDroughtWE_v2_1.pdf"), height = 4, dpi = 480)
dt <- plotTotalDT[model == "v2_1"]
dt[, diff := meanBaselineDrSum - meanBaselineWeSum]
dt[, pc_diff := 100*(meanBaselineDrSum - meanBaselineWeSum)/meanBaselineWeSum]
min <- min(dt$pc_diff)
head(dt[pc_diff < -3, .(obsDate, droughtStage, metered, pc = round(pc_diff,2))])
twe <- dt[, .(maxL = max(meanBaselineWeSum),
minL = min(meanBaselineWeSum)), keyby = .(metered)]
twe[, Model := "water efficiency"]
twe[, pc_reduction := 100*(1 - (minL/maxL))]
tdr <- dt[, .(maxL = max(meanBaselineDrSum),
minL = min(meanBaselineDrSum)), keyby = .(metered)]
tdr[, pc_reduction := 100*(1 - (minL/maxL))]
tdr[, Model := "water efficiency + drought"]
t <- rbind(twe, tdr)
kable(t[, .(Model, metered, maxL, minL, pc_reduction)], digits = 2)
``` ```
...@@ -1228,7 +1274,7 @@ ggsave(paste0("plots_v1/Fig4_5_DroughtPlotKey.pdf"), dpi = 400) ...@@ -1228,7 +1274,7 @@ ggsave(paste0("plots_v1/Fig4_5_DroughtPlotKey.pdf"), dpi = 400)
Compare 2_0 & 2_1 Compare 2_0 & 2_1
```{r} ```{r compare 2.0 and 2.1}
m2_0 <- plotTotalDT[model == "v2_0", .(obsDate, metered, m2_0_meanBaselineWeSum = meanBaselineWeSum, m2_0_meanBaselineDrSum = meanBaselineDrSum)] m2_0 <- plotTotalDT[model == "v2_0", .(obsDate, metered, m2_0_meanBaselineWeSum = meanBaselineWeSum, m2_0_meanBaselineDrSum = meanBaselineDrSum)]
m2_1 <- plotTotalDT[model == "v2_1", .(obsDate, metered, m2_1_meanBaselineWeSum = meanBaselineWeSum, m2_1_meanBaselineDrSum = meanBaselineDrSum)] m2_1 <- plotTotalDT[model == "v2_1", .(obsDate, metered, m2_1_meanBaselineWeSum = meanBaselineWeSum, m2_1_meanBaselineDrSum = meanBaselineDrSum)]
...@@ -1240,8 +1286,8 @@ models <- m2_0[m2_1] ...@@ -1240,8 +1286,8 @@ models <- m2_0[m2_1]
models[, weDiff := (m2_1_meanBaselineWeSum - m2_0_meanBaselineWeSum)/m2_0_meanBaselineWeSum] models[, weDiff := (m2_1_meanBaselineWeSum - m2_0_meanBaselineWeSum)/m2_0_meanBaselineWeSum]
models[, drDiff := (m2_1_meanBaselineDrSum - m2_0_meanBaselineDrSum)/m2_0_meanBaselineDrSum] models[, drDiff := (m2_1_meanBaselineDrSum - m2_0_meanBaselineDrSum)/m2_0_meanBaselineDrSum]
p <- ggplot2::ggplot(models, aes(x = obsDate, y = 100 * drDiff)) + p <- ggplot2::ggplot(models, aes(x = obsDate, colour = metered, y = 100 * drDiff)) +
geom_point() + geom_line() +
labs(y = "% difference", labs(y = "% difference",
x = "Date") x = "Date")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment