output$plotAllDat <- renderPlot({
estDat <- getEstDat()
regDat <- getRegDat()
datP <- allMISTEdat(regDat, estDat)
yvar <- paste0("y-", input$yID)
xvar <- paste0("x-", input$xID)
x2var <- paste0("x2-", input$xID2)
if(input$eventEst == FALSE) {
if (input$use2 == FALSE) {
cols <- c(yvar = "blue",
xvar = "green",
"Estimated" = "red")
if (input$log10 == FALSE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = -Inf, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x),], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper),], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower),], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
scale_y_continuous(labels = comma) +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual("", values = c("blue", "green", "red")) +
theme_bw() + theme(legend.title = element_blank())
}
else if (input$log10 == TRUE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = 0, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x),], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper),], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower),], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
scale_y_log10(labels = comma) +
annotation_logticks(sides = "rl") +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank(), panel.grid.minor = element_blank())
}
p
}
else if (input$use2 == TRUE) {
cols <- c(yvar = "blue",
xvar = "green",
x2var = "purple",
"Estimated" = "red")
if (input$log10 == FALSE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = -Inf, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x)], aes(x = dateTime, y = Flow.x, color = xvar), linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Flow.x2)], aes(x = dateTime, y = Flow.x2, color = x2var), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated)], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper)], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower)], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
scale_y_continuous(labels = comma) +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank())
}
else if (input$log10 == TRUE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = 0, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x)], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Flow.x2)], aes(x = dateTime, y = Flow.x2, color = x2var), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated)], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper)], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower)], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
scale_y_log10(labels = comma) +
annotation_logticks(sides = "rl") +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank(), panel.grid.minor = element_blank())
}
p
}
}
else if(input$eventEst == TRUE) {
if (input$giveObsQ == TRUE) {
addQs <- as.numeric(unlist(strsplit(input$ObsQ, ", ")))
addDTs <- unlist(strsplit(input$ObsQDate, ", "))
ObsDf <<- data.frame(dateTime = as.POSIXct(addDTs, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
Flow.obs = addQs)
cols <- c(yvar = "blue",
xvar = "green",
"Smoothed" = "orange",
"Estimated" = "red")
if (input$log10 == FALSE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = -Inf, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x),], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Smoothed, color = "Smoothed"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper),], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower),], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
geom_point(data = ObsDf, aes(x = dateTime, y = Flow.obs), alpha = 1, fill = "red", pch = 21, size = 3) +
scale_y_continuous(labels = comma) +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank())
}
else if (input$log10 == TRUE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = 0, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x),], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Smoothed, color = "Smoothed"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper),], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower),], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
geom_point(data = ObsDf, aes(x = dateTime, y = Flow.obs), alpha = 1, fill = "red", pch = 21, size = 3) +
scale_y_log10(labels = comma) +
annotation_logticks(sides = "rl") +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank(), panel.grid.minor = element_blank())
}
} else {
cols <- c(yvar = "blue",
xvar = "green",
"Estimated" = "red")
if (input$log10 == FALSE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = -Inf, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x),], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper),], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower),], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
scale_y_continuous(labels = comma) +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank())
}
else if (input$log10 == TRUE) {
p <- ggplot() +
geom_line(data = datP, aes(x = dateTime, y = Flow.y, color = yvar), size = 1) +
geom_rect(aes(xmin = as.POSIXct(input$estDateSt, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
xmax = as.POSIXct(input$estDateEn, format = "%Y-%m-%d %H:%M:%S", tz = "GMT"),
ymin = 0, ymax = Inf), alpha = 0.4, color = "grey80") +
geom_line(data = datP[!is.na(datP$Flow.x),], aes(x = dateTime, y = Flow.x, color = xvar), size = 1, linetype = "dotdash") +
geom_line(data = datP[!is.na(datP$Estimated),], aes(x = dateTime, y = Estimated, color = "Estimated"), size = 1, linetype = "longdash") +
geom_line(data = datP[!is.na(datP$fitUpper),], aes(x = dateTime, y = fitUpper), color = "black", size = 1, linetype = "dashed") +
geom_line(data = datP[!is.na(datP$fitLower),], aes(x = dateTime, y = fitLower), color = "black", size = 1, linetype = "dashed") +
scale_y_log10(labels = comma) +
annotation_logticks(sides = "rl") +
labs(x = "Date (UTC)", y = "Discharge, in cubic feet per second") +
#scale_color_manual(values = cols) +
theme_bw() + theme(legend.title = element_blank(), panel.grid.minor = element_blank())
}
}
p
}
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.