library(toxEval)
library(shiny)
library(shinydashboard)
library(magrittr)
library(leaflet)

cleaned_ep <- clean_endPoint_info(end_point_info) %>%
  dplyr::rename(endPoint = assay_component_endpoint_name)

trimmed_ep <- cleaned_ep
choicesPerGroup <- apply(cleaned_ep, 2, function(x) length(unique(x[!is.na(x)])))
choicesPerGroup <- choicesPerGroup[which(as.numeric(choicesPerGroup) > 6)]

groupChoices <- paste0(names(choicesPerGroup)," (",choicesPerGroup,")")

named_choices <- names(choicesPerGroup)
names(named_choices) <- groupChoices

ep <- data.frame(cleaned_ep[,c("endPoint", "intended_target_family")])
ep <- ep[!is.na(ep[,2]),]
ep <- ep[ep[,2] != "NA",]

orderBy <- ep[,2]
orderNames <- names(table(orderBy))
nEndPoints <- as.integer(table(orderBy))

df <- data.frame(orderNames,nEndPoints,stringsAsFactors = FALSE) %>%
  dplyr::arrange(desc(nEndPoints))

dropDownHeader <- c(paste0(df$orderNames," (",df$nEndPoints,")"))

selChoices <- df$orderNames

flagsALL <- c("Borderline active",
              "Only highest conc above baseline, active" ,      
              "Only one conc above baseline, active",
              "Noisy data",                                 
              "Hit-call potentially confounded by overfitting",
              "Gain AC50 < lowest conc & loss AC50 < mean conc",
              "Biochemical assay with < 50% efficacy",
              "Less than 50% efficacy",
              "AC50 less than lowest concentration tested",
              "Cell viability assay fit with gnls winning model")

shortFlags <- c("Borderline",
                "OnlyHighest",
                "OneAbove",
                "Noisy",
                "HitCall",
                "GainAC50",
                "Biochemical",
                "LessThan50",
                "ACCLessThan",
                "GNLSmodel")

assay_names <- c("Apredica" = "APR",
                 "Attagene" = "ATG",
                 "BioSeek" = "BSK",
                 "NovaScreen" = "NVS",
                 "Odyssey Thera" = "OT",
                 "Toxicity Testing" = "TOX21",
                 "CEETOX" = "CEETOX",
                 "CLD" = "CLD",
                 "TANGUAY" = "TANGUAY",
                 "NHEERL_PADILLA" = "NHEERL_PADILLA",
                 "NCCT" = "NCCT",
                 "ACEA Biosciences" = "ACEA",
                 "NHEERL_HUNTER" = "NHEERL_HUNTER", 
                 "NHEERL_NIS" = "NHEERL_NIS" , 
                 "NHEERL_MED" = "NHEERL_MED", 
                 "UPITT" = "UPITT")

init_text <- "######################################
# Setup:
library(toxEval)
#NOTE: Add path to path_to_file!!!
path_to_file <- 'Choose File'"

names(shortFlags) <- flagsALL

sidebar <- dashboardSidebar(
  sidebarMenu(
   fileInput("data", "Load Data File(Excel Worksheet)",multiple = FALSE),
   actionButton("exampleData", label = "Load Example Data"),
   radioButtons("radioMaxGroup", label = NULL,
                choices = list("Group" = 1, "Chemical" = 2, "Class" = 3), 
                selected = 3),
   radioButtons("meanEAR", choices = list("Mean EAR"=TRUE, "Max EAR" = FALSE),
                inline = TRUE, label = NULL, selected = FALSE),
   conditionalPanel(
     condition = "output.isTox == false",
     radioButtons("sumEAR", choices = list("Sum"=TRUE, "No Sum" = FALSE),
                  inline = TRUE, label = NULL, selected = TRUE)
   ),
   downloadButton('downloadBenchmarks', 'Download Benchmarks', style='margin-left:13px; color: #444'),
   menuItem("Assay", icon = icon("th"), tabName = "assay",
        checkboxGroupInput("assay", "Assays:",
                           assay_names,
          selected= c("ACEA", "APR", "ATG", 
                      "NVS", "OT",            
                      "TOX21", "CEETOX", "CLD", "TANGUAY", "NHEERL_PADILLA", "NCCT",          
                      "NHEERL_HUNTER", "NHEERL_NIS", "NHEERL_MED", "UPITT")),
          actionButton("pickAssay", label="Switch Assays")),
   menuItem("Annotation", icon = icon("th"), tabName = "annotation",
            selectInput("groupCol", label = "Annotation (# Groups)", 
                        choices = named_choices,
                        selected = "intended_target_family", 
                        multiple = FALSE),
            actionButton("changeAnn", label="Switch Annotation")
            ),
   menuItem("Group", icon = icon("th"), tabName = "groupMenu",
            checkboxGroupInput("group", "Groups (# End Points)",
                               setNames(df$orderNames,dropDownHeader),
                               selected=df$orderNames[c(-3)]),
            actionButton("allGroup", label="Switch Groups")),
   menuItem("Sites", icon = icon("th"), tabName = "siteMenu",
            selectInput("sites", label = "Site", 
                        choices = c("All"),
                        selected = "All", multiple = FALSE)
   ),
   menuItem("Flags", icon = icon("th"), tabName = "flagMenu",
            checkboxGroupInput("flags", "Include Flags",choices = shortFlags, selected = c("Borderline",
                                                                                           "OnlyHighest",
                                                                                           "GainAC50",
                                                                                           "Biochemical", 
                                                                                           "ACCLessThan")),
            actionButton("pickFlags", label="Switch flags")),
   menuItem("Hit Threshold",icon = icon("th"), tabName = "hitThresTab",
            numericInput("hitThres",label = "Hit Threshold",value = 0.1,min = 0.0000001),
            actionButton("changeHit", label="Change Hit Threshold")
   ),
   menuItem("Source code", icon = icon("fa-file-code-o"), 
            href = "https://code.usgs.gov/water/toxEval/tree/main/inst/shiny")
  )
)

body <- dashboardBody(
  includeCSS("www/custom.css"),
  h4(textOutput("title_text", inline = TRUE)),
  h3(textOutput("siteText")),
  tabBox(width = 12, id="mainOut",
    tabPanel(title = tagList("Map", shiny::icon("map-marker")),
             value="map",
             leafletOutput("mymap",height = "500px"),
             htmlOutput("mapFooter"),
             h4("R Code:"),
             shinyAce::aceEditor(outputId = "mapCode_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Box Plots", shiny::icon("chart-bar")),
             value="summary",
             fluidRow(
               column(3, checkboxInput("plot_thres_box", "Consider Threshold?", FALSE)),
               column(3, checkboxInput("plot_ND", "Plot non-detects's?", TRUE))
             ),
             uiOutput("graphGroup.ui", width = "100%"),
             fluidRow(
               column(3, downloadButton('downloadBoxPlot', 'Download PNG')),
               column(3, downloadButton('downloadBoxPlot_csv', 'Download CSV'))
             ),
             h4("R Code:"),
             shinyAce::aceEditor(outputId = "boxCode_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Bar Charts", shiny::icon("chart-bar")),
             value="summaryBar",
             plotOutput("stackBarGroup", width = "100%", height = "750px"),
             fluidRow(
               column(3, downloadButton('downloadStackPlot', 'Download PNG')),
               column(3, downloadButton('downloadStackPlot_csv', 'Download CSV'))
             ),
             h4("R Code:"),
             shinyAce::aceEditor(outputId = "barCode_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Max EAR and Frequency", shiny::icon("bars")),
             value="maxEAR",
            h3(textOutput("meanText")),
            h3(textOutput("freqText")),
            DT::dataTableOutput('tableSumm'),
            downloadButton('downloadTable', 'Download CSV'),
            h4("R Code:"),
            shinyAce::aceEditor(outputId = "tableSumm_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Hit Counts", shiny::icon("bars")),
             value="maxHits",
             h3(textOutput("nGroup")),
            DT::dataTableOutput('tableGroupSumm'),
            downloadButton('downloadGroupTable', 'Download CSV'),
            h4("R Code:"),
            shinyAce::aceEditor(outputId = "tableGroup_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Site Hits", shiny::icon("barcode")),
            value="siteHits",
            h3(textOutput("siteHitText")),
            div(DT::dataTableOutput("hitsTable"), style="font-size:90%"),
            downloadButton('downloadSiteHitTable', 'Download CSV'),
            h4("R Code:"),
            shinyAce::aceEditor(outputId = "siteHit_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Endpoints Hits", shiny::icon("barcode")),
             value="endHits",
             h3(textOutput("epHitTitle")),
             div(DT::dataTableOutput("hitsTableEPs"), style="font-size:90%"),
             downloadButton('downloadHitTable', 'Download CSV'),
             h4("R Code:"),
             shinyAce::aceEditor(outputId = "hitsTable_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Endpoint", shiny::icon("chart-bar")),
             value="endpoint",
             fluidRow(
               column(3,selectInput("epGroup", label = "Choose Chemical",
                                    choices = "All",
                                    multiple = FALSE,
                                    selected = "All")),
               column(2, checkboxInput("plot_thres_ep", "Consider Threshold?", FALSE)),
               column(2,numericInput("topNum",label = "Number to Show",value = 10, max = 50, min = 2))
             ),
            uiOutput("endpointGraph.ui"),
            fluidRow(
             column(3, downloadButton('downloadEndpoint', 'Download PNG')),
             column(3, downloadButton('downloadEndpoint_csv', 'Download CSV'))
            ),
            h4("R Code:"),
            shinyAce::aceEditor(outputId = "epGraph_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tabPanel(title = tagList("Heat Map", shiny::icon("chart-bar")),
                   value="heat",
             checkboxInput("plot_ND_heat", "Plot non-detects's?", TRUE),
             uiOutput("graphHeat.ui"),
             fluidRow(
               column(3, downloadButton('downloadHeatPlot', 'Download PNG')),
               column(3, downloadButton('downloadHeatPlot_csv', 'Download CSV'))
             ),
             h4("R Code:"),
             shinyAce::aceEditor(outputId = "heat_out", value = init_text, mode = "r", theme = "chrome", readOnly = TRUE)
    ),
    tags$head(tags$link(rel="shortcut icon", href="favicon.ico"))
  ),

  fluidRow(
    column(1, HTML('<svg xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink" version="1.1" width="300" height="75"><path id="USGS" d="m234.95 15.44v85.037c0 17.938-10.132 36.871-40.691 36.871-27.569 0-40.859-14.281-40.859-36.871v-85.04h25.08v83.377c0 14.783 6.311 20.593 15.447 20.593 10.959 0 15.943-7.307 15.943-20.593v-83.377h25.08m40.79 121.91c-31.058 0-36.871-18.27-35.542-39.03h25.078c0 11.462 0.5 21.092 14.282 21.092 8.472 0 12.62-5.482 12.62-13.618 0-21.592-50.486-22.922-50.486-58.631 0-18.769 8.968-33.715 39.525-33.715 24.42 0 36.543 10.963 34.883 36.043h-24.419c0-8.974-1.492-18.106-11.627-18.106-8.136 0-12.953 4.486-12.953 12.787 0 22.757 50.493 20.763 50.493 58.465 0 31.06-22.75 34.72-41.85 34.72m168.6 0c-31.06 0-36.871-18.27-35.539-39.03h25.075c0 11.462 0.502 21.092 14.285 21.092 8.475 0 12.625-5.482 12.625-13.618 0-21.592-50.494-22.922-50.494-58.631 0-18.769 8.969-33.715 39.531-33.715 24.412 0 36.536 10.963 34.875 36.043h-24.412c0-8.974-1.494-18.106-11.625-18.106-8.144 0-12.955 4.486-12.955 12.787 0 22.757 50.486 20.763 50.486 58.465 0 31.06-22.75 34.72-41.85 34.72m-79.89-46.684h14.76v26.461l-1.229 0.454c-3.816 1.332-8.301 2.327-12.453 2.327-14.287 0-17.943-6.645-17.943-44.177 0-23.256 0-44.348 15.615-44.348 12.146 0 14.711 8.198 14.933 18.107h24.981c0.198-23.271-14.789-36.043-38.42-36.043-41.021 0-42.52 30.724-42.52 60.954 0 45.507 4.938 63.167 47.12 63.167 9.784 0 25.36-2.211 32.554-4.18 0.436-0.115 1.212-0.596 1.212-1.216v-59.598h-38.612v18.09" style="fill:rgb(40%,40%,40%); fill-opacity: 0.3" transform="scale(0.5)"/>
  <path id="waves" d="m48.736 55.595l0.419 0.403c11.752 9.844 24.431 8.886 34.092 2.464 6.088-4.049 33.633-22.367 49.202-32.718v-10.344h-116.03v27.309c7.071-1.224 18.47-0.022 32.316 12.886m43.651 45.425l-13.705-13.142c-1.926-1.753-3.571-3.04-3.927-3.313-11.204-7.867-21.646-5.476-26.149-3.802-1.362 0.544-2.665 1.287-3.586 1.869l-28.602 19.13v34.666h116.03v-24.95c-2.55 1.62-18.27 10.12-40.063-10.46m-44.677-42.322c-0.619-0.578-1.304-1.194-1.915-1.698-13.702-10.6-26.646-5.409-29.376-4.116v11.931l6.714-4.523s10.346-7.674 26.446 0.195l-1.869-1.789m16.028 15.409c-0.603-0.534-1.214-1.083-1.823-1.664-12.157-10.285-23.908-7.67-28.781-5.864-1.382 0.554-2.7 1.303-3.629 1.887l-13.086 8.754v12.288l21.888-14.748s10.228-7.589 26.166 0.054l-0.735-0.707m68.722 12.865c-4.563 3.078-9.203 6.203-11.048 7.441-4.128 2.765-13.678 9.614-29.577 2.015l1.869 1.797c0.699 0.63 1.554 1.362 2.481 2.077 11.418 8.53 23.62 7.303 32.769 1.243 1.267-0.838 2.424-1.609 3.507-2.334v-12.234m0-24.61c-10.02 6.738-23.546 15.833-26.085 17.536-4.127 2.765-13.82 9.708-29.379 2.273l1.804 1.729c0.205 0.19 0.409 0.375 0.612 0.571l-0.01 0.01 0.01-0.01c12.079 10.22 25.379 8.657 34.501 2.563 5.146-3.436 12.461-8.38 18.548-12.507l-0.01-12.165m0-24.481c-14.452 9.682-38.162 25.568-41.031 27.493-4.162 2.789-13.974 9.836-29.335 2.5l1.864 1.796c1.111 1.004 2.605 2.259 4.192 3.295 10.632 6.792 21.759 5.591 30.817-0.455 6.512-4.351 22.528-14.998 33.493-22.285v-12.344" style="fill:rgb(40%,40%,40%); fill-opacity: 0.3" transform="scale(0.5)"/></svg>')
    ),
    column(11,
           h4("Disclaimer"),
           h5("This software is in the public domain because it contains materials that originally came from the U.S. Geological Survey (USGS), an agency of the United States Department of Interior. For more information, see the official USGS copyright policy at https://www.usgs.gov/visual-id/credit_usgs.html#copyright
              Although this software program has been used by the USGS, no warranty, expressed or implied, is made by the USGS or the U.S. Government as to the accuracy and functioning of the program and related program material nor shall the fact of distribution constitute any such warranty, and no responsibility is assumed by the USGS in connection therewith.
              This software is provided 'AS IS.'"))

  )
)

ui <- tagList(
  tags$header(class = "navbar",
              tags$a(href = 'https://www.usgs.gov/', type = "button",
                        img(src = 'logo.png', target="_blank",
                            title = "USGS", height = "60px"),
                     style = "float: left; 
                              padding: 10px 50px 10px 50px;"),
              # tags$a(
              #   type = "button",
              #   style="color: #fff; 
              #         font-weight: 300;
              #         font-size: 16px;
              #         font-family: 'Helvetica Neue',Helvetica,Arial,sans-serif;
              #         float: left; 
              #         padding: 30px 30px 30px 10px;",
              #   href = "https://rconnect.usgs.gov/toxEval_docs/articles/Introduction.html", 
              #   target="_blank",
              #   "Introduction"
              # ),
              # tags$a(
              #   type = "button",
              #   # class = "btn action-button",
              #   style="color: #fff;
              #         font-weight: 300;
              #         font-size: 16px;
              #         font-family: 'Helvetica Neue',Helvetica,Arial,sans-serif;
              #         float: left; 
              #         padding: 30px 30px 20px 10px;",
              #   href = "https://rconnect.usgs.gov/toxEval_docs/articles/basicWorkflow.html",
              #   target="_blank",
              #   "Workflow"
              # ),
              # tags$a(
              #   style="color: #fff;
              #         font-weight: 300;
              #         font-size: 16px;
              #         font-family: 'Helvetica Neue',Helvetica,Arial,sans-serif;
              #         float: left; 
              #         padding: 30px 30px 20px 10px;",
              #   href = "https://rconnect.usgs.gov/toxEval_docs/articles/PrepareData.html",
              #   target="_blank",
              #   "Prepare Data"
              # ),
              # tags$a(
              #   style="color: #fff;
              #         font-weight: 300;
              #         font-size: 16px;
              #         vertical-align: middle;
              #         font-family: 'Helvetica Neue',Helvetica,Arial,sans-serif;
              #         float: left; 
              #         padding: 30px 30px 20px 10px;",
              #   href = "https://rconnect.usgs.gov/toxEval_docs/articles/shinyApp.html",
              #   target="_blank",
              #   "Using the Shiny App"
              # ),
              tags$button(
                id = 'close',
                type = "button",
                class = "btn action-button",
                style='color: #000000;
                       float: right;
                       margin-right:13px;margin-top:20px;margin-bottom:7px',
                onclick = "setTimeout(function(){window.close();},500);",  # close browser
                "Stop ToxEval"
              )),
  dashboardPage(header = dashboardHeader(disable = TRUE),
                sidebar = sidebar,
                body = body),
  # tags$footer(htmltools::includeHTML("www/footer.html") )
)
