inst/Factocatdesapp/server.R

function(input, output, session) {
  
  resultat <- reactive({
    return(catdes(donnee = my_data,
                  num.var = which(colnames(my_data) == input$select_categorical_var),
                  proba = 1
    ))
  })
  
  output$select_quali <- renderUI({
    selectInput(
      inputId = "select_quali_var_test",
      label = gettext("Qualitative variables",domain="R-Factoshiny"),
      multiple = TRUE,
      choices = colnames(my_data)[sapply(my_data,is.factor) & !(colnames(my_data) %in% input$select_categorical_var)],
      selected = choix_var_quali
    )
  })
  
  tableau_vtest <- reactive({
    
    if(!is.null(resultat()$quanti)){
      lvl <- length(resultat()$quanti)
      x <- resultat()$quanti[[1]][,"v.test", drop = FALSE]
      tabvtest <- x[sort(rownames(x)),,drop=FALSE]
      
      x <- resultat()$quanti[[1]][,"p.value", drop = FALSE]
      tabpvalue <- x[sort(rownames(x)),,drop=FALSE]
      
      x <- resultat()$quanti[[1]][,2, drop = FALSE]
      tabmean <- x[sort(rownames(x)),,drop=FALSE]
      
      for(i in 2:lvl){
        x <- resultat()$quanti[[i]][,"v.test", drop = FALSE]
        tabvtest <- cbind(tabvtest,as.data.frame(x[sort(rownames(x)),]))
        x <- resultat()$quanti[[i]][,"p.value", drop = FALSE]
        tabpvalue <- cbind(tabpvalue,as.data.frame(x[sort(rownames(x)),]))
        x <- resultat()$quanti[[i]][,2, drop = FALSE]
        tabmean <- cbind(tabmean,as.data.frame(x[sort(rownames(x)),]))
      }
      colnames(tabpvalue) <- colnames(tabvtest) <- colnames(tabmean) <- names(resultat()$quanti)
      tabmean$overall <- resultat()$quanti[[1]][sort(rownames(resultat()$quanti[[1]])),3]
      
      validate(
        need(as.numeric(input$select_proba_plot) > 0, paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),0))
      )
      if (input$select_proba_plot <= min(tabpvalue)) return(NULL)
      else {
	    sortie <- signif(t(as.matrix(tabvtest[apply(tabpvalue,1,min) <= input$select_proba_plot & rownames(tabvtest) %in% input$select_quanti_var,])),3)
        return(sortie)
      }
    }
  })
  
  tableau_quali <- reactive({
    
    if(!is.null(resultat()$category) & !(input$select_categorical_var %in% input$select_quali_var_test)){
      coldonnee <- colnames(my_data[input$select_quali_var_test])
      
      rows <- NULL
      for (i in 1:nrow(resultat()$category[[1]])) {
        if (strsplit(rownames(resultat()$category[[1]]),"=")[[i]][1]%in%coldonnee) rows <- c(rows,rownames(resultat()$category[[1]])[i])
      }
      lvl <- length(resultat()$category)
      x <- resultat()$category[[1]][rows,"v.test", drop = FALSE]
      tabvtest <- as.data.frame(x[sort(rownames(x)),,drop=FALSE])
      x <- resultat()$category[[1]][rows,"p.value", drop = FALSE]
      tabpvalue <- as.data.frame(x[sort(rownames(x)),,drop=FALSE])
      for(i in 2:lvl){
        x <- resultat()$category[[i]][rows,"v.test", drop = FALSE]
        tabvtest <- cbind(tabvtest,as.data.frame(x[sort(rownames(x)),,drop=FALSE]))
        x <- resultat()$category[[i]][rows,"p.value", drop = FALSE]
        tabpvalue <- cbind(tabpvalue,as.data.frame(x[sort(rownames(x)),,drop=FALSE]))
      }
      colnames(tabpvalue) <- colnames(tabvtest) <- names(resultat()$category)
      validate(
        need(as.numeric(input$select_proba_plot) > 0, paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),0))
      )
      # validate(
        # need(input$select_proba_plot > min(tabpvalue), paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),signif(min(tabpvalue),3)))
      # )
      if (input$select_proba_plot <= min(tabpvalue)) return(NULL)
	  else {
	    sortie <- signif(tabvtest[apply(tabpvalue,1,min) <= input$select_proba_plot,,drop = FALSE],3)
        return(sortie)
      }
    }
    
  })
  
  output$barplot <-  renderPlot({
    plot.catdes(x = resultat(), 
                level = input$select_proba_plot,
                col.upper = input$col_up,
                col.lower = input$col_low,
                barplot = TRUE
    )
  })
  
  tab_quanti <- reactive({
    if(!is.null(resultat()$quanti)){
      if(!is.null(input$select_quanti_var)){
        lvl <- length(resultat()$category)
        tabVtest <- tableau_vtest()
		mini <- 1
		if (!is.null(resultat()$quanti.var)) mini <- min(resultat()$quanti.var[,2])
		# if (!is.null(resultat()$test.chi2)) mini <- min(mini,resultat()$test.chi2[,1])
        validate(
          need(!is.null(tabVtest), paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),mini,"."))
        )
        quant <- seq(min(tabVtest, na.rm = T), max(tabVtest, na.rm = T),length.out = 100)
        color <- grDevices::colorRampPalette(c(input$col_low,"white",input$col_up))(length(quant)+1)
        
        a <- DT::formatStyle(
          DT::datatable(t(as.matrix(tabVtest)),
                        extensions = c('Buttons','FixedColumns','FixedHeader'),
                        options = list( pageLength = ncol(tabVtest),
                            dom = 'Bfrtip', buttons = c('csv'), fixedColumns = TRUE, fixedHeader = TRUE)
          ),
          rownames(tabVtest),
          backgroundColor = DT::styleInterval(quant, color)
        )
        return(a)
      }
    }
  })
  
  
  output$tableau_df_quanti <- DT::renderDataTable({
    if(!is.null(resultat()$quanti)) return(tab_quanti())
  })
  
  
  tab_quali <- reactive({
    if(!is.null(input$select_quali_var_test)){
      
      if(!(input$select_categorical_var %in% input$select_quali_var_test)){
        tabQuali <- tableau_quali()
        validate(
          need(!is.null(tabQuali), paste(gettext("The p-value is too small. You should increase the p-value.",domain="R-Factoshiny")))
        )
        quant <- seq(min(tabQuali,na.rm = TRUE), max(tabQuali, na.rm = TRUE), length.out = 100)
        color <- grDevices::colorRampPalette(c(input$col_low,"white",input$col_up))(length(quant)+1)
        
        a <- DT::formatStyle(
          DT::datatable(tabQuali, extensions = c('Buttons','FixedColumns','FixedHeader'),
            options = list(pageLength = nrow(tabQuali), dom = 'Bfrtip',
              buttons = c('csv'), fixedColumns = TRUE, fixedHeader = TRUE)
          ),
          colnames(tabQuali), backgroundColor = DT::styleInterval(quant, color)
        )  
        return(a)
      }
    }
    
  })
  
  output$tableau_df_quali <- DT::renderDataTable({
    if (!is.null(tab_quali())) return(tab_quali())
  })
  
  tab_both <- reactive({
    tab <- NULL
    # if(!is.null(input$select_quali_var_test) & !is.null(input$select_quanti_var)) tab <- rbind(tableau_quali(), t(tableau_vtest()))
    # if(!is.null(input$select_quali_var_test) & is.null(input$select_quanti_var)) tab <- rbind(tableau_quali())
    # if(is.null(input$select_quali_var_test) & !is.null(input$select_quanti_var)) tab <- rbind(t(tableau_vtest()))
    if(!is.null(input$select_quali_var_test)) tab <- rbind(tab,tableau_quali())
    if(!is.null(input$select_quanti_var) & !is.null(tableau_vtest())) tab <- rbind(tab,t(tableau_vtest()))
	mini <- 1
	if (!is.null(resultat()$quanti.var)) mini <- min(resultat()$quanti.var[input$select_quanti_var,2])
	if (!is.null(resultat()$category)) mini <- min(mini,resultat()$test.chi2[input$select_quali_var_test,1])
    validate(
      need(!is.null(tab), paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),signif(mini,3),"."))
    )
      # validate(
        # need(!is.null(tab), paste(gettext("The p-value is too small. You should increase the p-value.",domain="R-Factoshiny")))
      # )

    if (!is.null(tab)){
      quant <- seq(min(tab, na.rm = T), max(tab, na.rm = T), length.out = 100)
      color <- grDevices::colorRampPalette(c(input$col_low,"white",input$col_up))(length(quant)+1)
      
      a <- DT::formatStyle(
        DT::datatable(tab,extensions = c('Buttons','FixedColumns','FixedHeader'),
          options = list(pageLength = nrow(tab),dom = 'Bfrtip',
            buttons = c('csv'),fixedColumns = TRUE,fixedHeader = TRUE)
        ), colnames(tab),backgroundColor = DT::styleInterval(quant, color)
      )
      return(a)
    }
  })
  
  output$tableau_df_both <- DT::renderDataTable({
    if (!is.null(tab_both())) return(tab_both())
  })
  
  observe({ input$catdesMAJ
    output$resu_catdes <- renderPrint({ isolate(catdes(donnee = my_data[,c(input$select_categorical_var,input$select_quali_var_test,input$select_quanti_var)],
                                                       num.var = 1, proba = input$select_proba_plot))
    }) 
  })
  
  tableau_link_quanti <- reactive({
    if(!is.null(input$select_quanti_var)){
      validate(
        need(input$select_proba_plot > min(resultat()$quanti.var[,"P-value"]), paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),signif(min(resultat()$quanti.var[,"P-value"]),3)))
      )
      tab <- as.data.frame(resultat()$quanti.var[rownames(resultat()$quanti.var) %in% input$select_quanti_var  & resultat()$quanti.var[,"P-value"] <= input$select_proba_plot,c("Eta2","P-value"), drop = FALSE])
      return(tab)
    }
  })
  
  df_link_quanti <- reactive({
    if(!is.null(input$select_quanti_var)){
      quant <- seq(min(tableau_link_quanti()[,1]), max(tableau_link_quanti()[,1]), length.out = 100)
      color <- grDevices::colorRampPalette(c(input$col_low,"white",input$col_up))(length(quant)+1)
      a <- DT::formatStyle(
        DT::datatable(signif(tableau_link_quanti(),3),
                      extensions = c('Buttons','FixedColumns','FixedHeader'),
                      options = list(pageLength = nrow(tableau_link_quanti()),
                        dom = 'Bfrtip',buttons = c('csv'),fixedColumns = TRUE,fixedHeader = TRUE)
        ),
        columns = colnames(tableau_link_quanti()),
        valueColumns = 'P-value',
        backgroundColor = DT::styleInterval(quant, color)
      )
    }
    return(a)
  })
  
  output$table_link_quanti <- DT::renderDataTable({
    if (!is.null(df_link_quanti())) return(df_link_quanti())
  })
  
  tableau_link_chisquare <- reactive({
    if(!is.null(resultat()$test.chi2)){
      
      validate(
        need(input$select_proba_plot > min(resultat()$test.chi2[,"p.value"]), paste(gettext("The p-value should be greater than",domain="R-Factoshiny"),signif(min(resultat()$test.chi2[,"p.value"]),3)))
      )
      tab <- (resultat()$test.chi2[rownames(resultat()$test.chi2) %in% input$select_quali_var_test & resultat()$test.chi2[,"p.value"] <= input$select_proba_plot,"p.value", drop = FALSE])
      return(tab)
    }
  })
  
  tab_chisquare <- reactive({
    if(!is.null(resultat()$test.chi2)){
      if(nrow(tableau_link_chisquare()) > 0){
        
        quant <- seq(min(tableau_link_chisquare()[,"p.value"]), max(tableau_link_chisquare()[,"p.value"]), length.out = 100)
        color <- grDevices::colorRampPalette(c(input$col_up,"white",input$col_low))(length(quant)+1)
        
        a <- DT::formatStyle(
          DT::datatable(
            signif(tableau_link_chisquare(),3),
            extensions = c('Buttons','FixedColumns','FixedHeader'),
            options = list(pageLength = nrow(tableau_link_chisquare()),
                           dom = 'Bfrtip', buttons = c('csv'), fixedColumns = TRUE, fixedHeader = TRUE)
          ),
          
          columns = colnames(tableau_link_chisquare()),
          valueColumns = "p.value",
          backgroundColor = DT::styleInterval(quant, color)
        )
        return(a)
      }
    }
  })
  
  output$table_link_chisquare <- DT::renderDataTable({
    if(!is.null(tab_chisquare())) return(tab_chisquare())
  })
  
  output$quanti_quali_both1 <- renderUI({
    xx <- gettext("Quantitative",domain="R-Factoshiny")
    old.x <- gettext("Both",domain="R-Factoshiny")
    if(!is.null(input$quanti_quali_both)) old.x <- input$quanti_quali_both
    if(!is.null(input$select_quali_var_test) & !is.null(input$select_quanti_var)) xx <- c(gettext("Both",domain="R-Factoshiny"),gettext("Quantitative",domain="R-Factoshiny"),gettext("Qualitative",domain="R-Factoshiny"))
    if(is.null(input$select_quali_var_test) & !is.null(input$select_quanti_var)) xx <- c(gettext("Quantitative",domain="R-Factoshiny"))
    if(!is.null(input$select_quali_var_test) & is.null(input$select_quanti_var)) xx <- c(gettext("Qualitative",domain="R-Factoshiny"))
    
    if (old.x%in%xx){
      radioButtons(inputId = "quanti_quali_both",choices = xx,
        inline = TRUE,label = gettext("Describe by ... variables",domain="R-Factoshiny"),selected = old.x)
    } else {
      radioButtons(inputId = "quanti_quali_both", choices = xx,
        inline = TRUE,label = gettext("Describe by ... variables",domain="R-Factoshiny"))
    }
  })
  
  observeEvent(input$download_tabquanti,{
      NameFile <- tcltk::tclvalue(tcltk::tcl("tk_getSaveFile"))
      if (!(any(strsplit(NameFile,split="[.]")[[1]]=="html")) & !(any(strsplit(NameFile,split="[.]")[[1]]=="htm"))) NameFile <- paste0(NameFile,".html")
      htmlwidgets::saveWidget(widget =  df_link_quanti(), file = NameFile)
  })
  
  observeEvent(input$download_tabquali,{
      NameFile <- tcltk::tclvalue(tcltk::tcl("tk_getSaveFile"))
      if (!(any(strsplit(NameFile,split="[.]")[[1]]=="html")) & !(any(strsplit(NameFile,split="[.]")[[1]]=="htm"))) NameFile <- paste0(NameFile,".html")
      htmlwidgets::saveWidget(widget =  tab_chisquare(), file = NameFile)
  })
  
  observeEvent(input$download_cate_quanti,{
      NameFile <- tcltk::tclvalue(tcltk::tcl("tk_getSaveFile"))
      if (!(any(strsplit(NameFile,split="[.]")[[1]]=="html")) & !(any(strsplit(NameFile,split="[.]")[[1]]=="htm"))) NameFile <- paste0(NameFile,".html")
      htmlwidgets::saveWidget(widget =  tab_quanti(), file = NameFile)
  })
  
  observeEvent(input$download_cate_quali,{
      NameFile <- tcltk::tclvalue(tcltk::tcl("tk_getSaveFile"))
      if (!(any(strsplit(NameFile,split="[.]")[[1]]=="html")) & !(any(strsplit(NameFile,split="[.]")[[1]]=="htm"))) NameFile <- paste0(NameFile,".html")
      htmlwidgets::saveWidget(widget =  tab_quali(), file = NameFile)
  })
  
  observeEvent(input$download_cate_both,{
      NameFile <- tcltk::tclvalue(tcltk::tcl("tk_getSaveFile"))
      if (!(any(strsplit(NameFile,split="[.]")[[1]]=="html")) & !(any(strsplit(NameFile,split="[.]")[[1]]=="htm"))) NameFile <- paste0(NameFile,".html")
      htmlwidgets::saveWidget(widget =  tab_both(), file = NameFile)
    
  })
  
  liste_retourner <- reactive({
    retour <- list()
    retour$donnees = my_data
    retour$explain = input$select_categorical_var
    retour$proba = input$select_proba_plot
    retour$col_basse = input$col_low
    retour$col_haute = input$col_up
    retour$var_quanti = input$select_quanti_var
    retour$var_quali = input$select_quali_var_test
    class(retour) <- c("catdesshiny", "list")
    
    return(retour)
  })
  
  observeEvent(input$Quit,{
        stopApp(returnValue=liste_retourner())
  })
  
}
husson/Factoshiny documentation built on Sept. 19, 2023, 5:05 a.m.