установка реактивных переменных в Shiny App

Question

Я пытаюсь написать свое первое блестящее приложение.Я получил много помощи с этого сайта, но снова застрял и думаю, что это связано с тем, что мои реактивные переменные не инициализируются.

У меня проблемы с функцией создания таблицы - anova_tab2.Если я не включаю раздел if (compFacts == "Metarhizium) {}, то он работает нормально, но как только я добавляю операторы if, я получаю ошибку об" аргументе нулевой длины ", который, я думаю,означает, что if (compFacts ==) не работает, я предполагаю, потому что переменная compFacts пуста?

Я включил только часть своей серверной функции, вызывающую проблему. Спасибо

Вот мой частичный сценарий:

ui <- fluidPage(
  # Make a title to display in the app
  titlePanel(" Exploring the Effect of Metarhizium on the Soil and Root Microbiome "),
  # Make the Sidebar layout
  sidebarLayout(
    # Put in the sidebar all the input functions
    sidebarPanel(
      tabsetPanel(id="tabs",
        tabPanel("otu", selectInput('dataset', 'dataset', names(abundance_tables),selected=names(abundance_tables)[1]),
                 uiOutput("otu"), br(),
                 # Add comment
                 p("For details on OTU identification please refer to the original publications")),
        tabPanel("anova", sliderInput('pval','p-value for significance',
                                      value=0.1,min=0,max=0.5,step=0.00001),
                 selectInput('dataset', 'dataset', names(abundance_tables)),
                 selectInput('fact_ofInt','factor of interest',FactorsOfInt,selected="Metarhizium"))
        ) 
    ),
    # Put in the main panel of the layout the output functions 
      mainPanel(
        conditionalPanel(condition="input.tabs == 'otu'",
                         plotOutput('plot'),
                         tableOutput("table1"),
                         p("anova for Soil Samples"),
                         tableOutput("tableS"),
                         tableOutput("tableR")

        ),
        conditionalPanel(condition="input.tabs == 'anova'",
                         #plotOutput('plot2')
                        # verbatimTextOutput("anovaText")
                        tableOutput("anova_tab2")

        )
      )
  )
)

# ========================= SERVER ####
  server <- function(input, output){
    # Return the requested dataset ----
    datasetInput <- reactive({
      abundance_tables[[input$dataset]]
    })
    pvalInput<-reactive({
      input$pval
   })
    comparisonInput<-reactive({
      input$FactorsOfInt
    })
    #
    # output otus to choose basaed on dataset selection
   output$otu <- renderUI({
     selectInput(inputId = "otu", label = "otu",
                       choices = colnames(datasetInput()),selected="Metarhizium")
    })
   otuInput<-reactive({
     input$otu
   })

  output$anova_tab2 <- renderTable({
    pval<-pvalInput()
    df<-datasetInput()
    compFacts<-comparisonInput()
    print(paste("compFacts:",compFacts))
    df_annot<-merge(df,sample_metadata,by="row.names",all.x=T)
    rownames(df_annot)<-df_annot[,1]
    df_annot<-df_annot[,-1]
    df_annot<-subset(df_annot,df_annot$Bean =="Bean")
    sum_tab.sig<-data.frame(Df=as.numeric(),Sum.Sq=as.numeric(),Mean.Sq=as.numeric(),Fval=as.numeric(),Pval=as.numeric(),OTU=as.character())
      for (i in 1:ncol(df)) {
      #  if (compFacts=="Metarhizium"){
          aov.ex <- aov(df_annot[,i]~df_annot$Fungi)
       # } 
        #else if (compFacts=="Insect"){
        #  aov.ex <- aov(df_annot[,i]~df_annot$Insect)
        #} else if (compFacts=="Sample_Type"){
        #  aov.ex <- aov(df_annot[,i]~df_annot$Location)
        #}
        ## summary of the anova given as a list
        sum_tab<-summary(aov.ex)[[1]]
        temp.sig<-sum_tab[sum_tab[,"Pr(>F)"]<pval,]  
        temp.sig<-transform(temp.sig,OTU=colnames(df)[i])
           sum_tab.sig<-rbind(sum_tab.sig,temp.sig)
          sum_tab.sig<-sum_tab.sig[!is.na(sum_tab.sig[,1]),]
          if (dim(sum_tab.sig)[1]==0){
            sum_tab.sig<-c(rep(0,5),"no significant trends at this pvalue")
          }
      }
    sum_tab.sig<-sum_tab.sig[-1,c(6,4,5)]
    colnames(sum_tab.sig)<-c("OTU","F-value","P-value")
    sum_tab.sig$'P-value'<-format(sum_tab.sig$'P-value',digits=3,scientific=5)
  return(sum_tab.sig)
  })
   ### end of server
  }
  ##