R/mod_star_fusion.R

In VincentAlcazer/hemRNA: A package for quick RNAseq data analysis

#' star_fusion UI Function
#'
#' @description A shiny Module.
#'
#' @param id,input,output,session Internal parameters for {shiny}.
#'
#' @noRd
#'
#' @importFrom shiny NS tagList
mod_star_fusion_ui <- function(id){
  ns <- NS(id)
  tagList(

    fluidPage(
      tabsetPanel(
        id = "fusion", type = "tabs",
        tabPanel("Overall graph",

                 column(8, shinycssloaders::withSpinner(plotOutput(ns("tileplot")),type=6)
                 ),
                 column(1)


        ),#tabsetpanel
        tabPanel("Result table",
                 column(8,
                        downloadButton(ns("download_table"), "Download table (.tsv)"),
                        shinycssloaders::withSpinner(DT::DTOutput(ns("result_table")),type=6)
                 ),
                 column(1)
        )
      ),##tabsetpanel

      column(2,
             absolutePanel(
               width = 200, right = 20, draggable = T,
               style = "opacity: 0.85",
               wellPanel(
                 numericInput(ns("FFPM"), label = "Min. Fusion Fragment per million (FFPM)",
                              min = 0, value = 0.5),
                 numericInput(ns("JunctionReadCount"), label = "Min. JunctionReadCount",
                              min = 0, value = 1),
                 numericInput(ns("SpanningFragCount"), label = "Min. SpanningFragCount",
                              min = 0, value = 1),
                 sliderInput(ns("y_size"), label = "y-axis font size",
                             min = 1, max = 30, value = 12, step= 1),
                 sliderInput(ns("x_size"), label = "x-axis font size",
                             min = 1, max = 30, value = 14, step= 1),

                 selectInput(ns("legend_ext"),
                             label = ("External legend"),
                             choices = c(
                               "No" = "none",
                               "Top" = "top",
                               "Right" = "right",
                               "Left" = "left",
                               "Bottom" = "bottom"
                             ),
                             multiple = F, selected = "right"
                 )
               )
             ) # Absolutepanel
      ) # Column

    ) # Fluidpage
  )
}

#' star_fusion Server Functions
#'
#' @noRd
mod_star_fusion_server <- function(id, r){
  moduleServer( id, function(input, output, session){
    ns <- session$ns


    df_fusion <- reactive({r$test$df_starfusion})

    plot <- reactive({

      req(df_fusion())

      plot <- df_fusion() %>%
        filter(FFPM >= input$FFPM,
               JunctionReadCount >= input$JunctionReadCount,
               SpanningFragCount >= input$SpanningFragCount) %>%
        ggplot(aes(x = sample_id, y = fusion, fill = FFPM)) +
        geom_tile() +
        scale_fill_viridis_c() +
        labs(x="", y="")

      return(plot)

    })

    height <- reactive({
      if(length(unique(df_fusion()[,"fusion"]) < 20)){
        height = 800
      } else if(length(unique(df_fusion()[,"fusion"]) < 40)){
        height = 1200
      } else if(length(unique(df_fusion()[,"fusion"]) < 60)){
        height = 1600
      } else {
        height = 2000
      }


    })


    output$tileplot <- renderPlot(height = height,{
      plot() +
        default_theme +
        theme(
          legend.position = input$legend_ext,
          axis.text.y = element_text(size = input$y_size, face = "bold"),
          axis.text.x = element_text(size = input$x_size, angle=90, vjust = 0.5, hjust = 1, face = "bold"))
    })

    output$result_table <- DT::renderDT(

      df_fusion() %>%
        select(sample_id, gene1, gene2, JunctionReadCount,SpanningFragCount, FFPM,  type, reading_frame, everything()) %>%
        filter(FFPM >= input$FFPM,
               JunctionReadCount >= input$JunctionReadCount,
               SpanningFragCount >= input$SpanningFragCount), # data
      class = "display nowrap compact", # style
      filter = "top", # location of column filters
      server = T,
      rownames = FALSE,
      options = list(
        scrollX = TRUE,
        "pagelength" = 20,
        lengthChange = TRUE,
        columnDefs = list(list(className = "dt-left", targets = "_all"))
      )
    )


    output$download_table <- downloadHandler(
      filename = function() {
        paste("Fusioncatcher.tsv")
      },
      content = function(file) {
        write.table(df_fusion() %>%
                      select(sample_id, gene1, gene2, FFPM, type, reading_frame, everything()),
                    file, row.names = FALSE, sep = "\t", quote = F)
      }
    )


  })
}

## To be copied in the UI
# mod_star_fusion_ui("star_fusion_ui_1")

## To be copied in the server
# mod_star_fusion_server("star_fusion_ui_1")