inst/shiny/plotXY1.R

In bbreaker-usgs/MISTEuv: Web app for estimation of unit value streamflow data

output$plotXY1 <- renderPlot({
  
  estDat <- getEstDat()
  
  regDat <- getRegDat()
  
  datP <- allMISTEdat(estDat, regDat)
  
  if (input$Method == 1) {
    
    p <- ggplot(data = datP, aes(x = Flow.x, y = Flow.y)) +
      geom_point(size = 3) +
      scale_y_log10() +
      scale_x_log10() +
      annotation_logticks(sides = "trbl") +
      stat_smooth(method = "lm") + 
      labs(x = paste0("Discharge, in cubic feet per second \n for ", input$xID), 
           y = paste0("Discharge, in cubic feet per second \n for ", input$yID)) + 
      theme_bw() + theme(legend.title = element_blank())
    
    p
    
  }
  
  else if (input$Method == 2) {
    
    if (input$adjKnots == FALSE) {
      
      p <- ggplot(data = datP, aes(x = Flow.x, y = Flow.y)) +
        geom_point() +
        scale_y_log10() +
        scale_x_log10() +
        annotation_logticks(sides = "trbl") +
        stat_smooth(method = "gam", formula = y ~ s(x, bs = "cs")) +
        labs(x = paste0("Discharge, in cubic feet per second \n for ", input$xID), 
             y = paste0("Discharge, in cubic feet per second \n for ", input$yID)) + 
        theme_bw() + theme(legend.title = element_blank())
      
      p
      
    }
    
    else if (input$adjKnots == TRUE) {
      
      p <- ggplot(data = datP, aes(x = Flow.x, y = Flow.y)) +
        geom_point() +
        scale_y_log10() +
        scale_x_log10() +
        annotation_logticks(sides = "trbl") +
        stat_smooth(method = "gam", formula = y ~ s(x, bs = "cr", k = input$knots)) +
        labs(x = paste0("Discharge, in cubic feet per second \n for ", input$xID), 
             y = paste0("Discharge, in cubic feet per second \n for ", input$yID)) + 
        theme_bw() + theme(legend.title = element_blank())
      
      p
      
    }
    
  }
  
})