Add functionality to view trait data for a given variable

Author: Nick Heyek

experiment-trait-data-visualizer/app.R

@@ -1,6 +1,7 @@
 library(shiny)
 library(traits)
 library(ggplot2)
+library(lubridate)
 
 knitr::opts_chunk$set(echo = FALSE, cache = TRUE)
 
@@ -9,14 +10,14 @@ options(betydb_key = readLines('~/.betykey', warn = FALSE),
         betydb_api_version = 'beta')
 
 experiments <- betydb_query(table='experiments')
-#varIdsObserved <- as.numeric(unique(traitData$variable_id))
+#varIdsObserved <- as.numeric()
 
 ui <- fluidPage(
   titlePanel("BETYdb Trait Data"),
   sidebarLayout (
     sidebarPanel(
-      selectInput('selectedExp', 'Experiment', experiments$name)
-      #selectInput('selectedVariable','Variable', varIdsObserved)
+      selectInput('selectedExp', 'Experiment', experiments$name),
+      uiOutput('selectVariable')
     ),
     mainPanel(
       plotOutput('traitPlot')
@@ -28,20 +29,26 @@ server <- function(input, output) {
 
   output$traitPlot <- renderPlot({
 
-    fullTraitData = data.frame()
+    fullTraitData <- data.frame()
 
     selectedExpRow <- subset(experiments, name==input$selectedExp)
     experimentStartDate <- as.Date(selectedExpRow$start_date)
     experimentEndDate <- as.Date(selectedExpRow$end_date)
 
     currDate <- experimentStartDate
     while (experimentEndDate - currDate != 0) {
-      currTraitData <- betydb_query(table='traits', date=paste0('~', currDate), limit='5')
+      currTraitData <- betydb_query(table='traits', date=paste0('~', currDate), limit='1')
       fullTraitData <- rbind(fullTraitData, currTraitData)
       currDate <- currDate + days(1)
     }
 
-    qplot(as.Date(fullTraitData$date), fullTraitData$mean, main="[Variable] for [Experiment]", 
+    output$selectedVariable <- renderUI({
+      variableIds <- as.numeric(unique(fullTraitData$variable_id))
+      selectInput <- selectInput('selectedVariable', 'Variable ID', variableIds)
+    })
+    variableTraitData <- subset(fullTraitData, variable_id=input$selectedVariable)
+    
+    qplot(as.Date(variableTraitData$date), variableTraitData$mean, main="[Variable] for [Experiment]", 
           xlab="Date", ylab="Unit", geom="auto")
   })
 }