In bioinformatist/LncPipe-Reporter: Automatically Aggregating and Summarizing lncRNA Analysis Results for Interactive Report

knitr::opts_chunk$set(
    echo = FALSE,
    message = FALSE,
    warning = FALSE,
    cache = FALSE,
    fig.height = 6,
    fig.width = 8
)

local({
  hook_plot = knitr::knit_hooks$get('plot')
  knitr::knit_hooks$set(plot = function(x, options) {
    x = paste(x, collapse = '.')
    if (!grepl('\\.svg', x)) return(hook_plot(x, options))
    # read the content of the svg image and write it out without <?xml ... ?>
    paste(readLines(x)[-1], collapse = '\n')
  })
})

options(DT.options = list(autoFill = TRUE, 
                          searchHighlight = TRUE, 
                          search = list(regex = TRUE, caseInsensitive = FALSE),
                          dom = 'Bfrtip', 
                          buttons = c('copy', 'csv', 'excel', 'pdf', 'print')))
# rm(list = setdiff(ls(), "params"))
# gc()
# The latest plotly should be used, for a lot of bugs fixed
# if (!require("devtools")) install.packages("devtools")
# install.packages(c("curl", "httr"))
# devtools::install_github('ramnathv/htmlwidgets')
# devtools::install_github("ropensci/plotly")
# source("https://bioconductor.org/biocLite.R")
# biocLite("edgeR")
# library(devtools)
# install_github("vqv/ggbiplot")
# require(plotly)
# require(data.table)
# require(cowplot)
# require(stringr)
# require(knitr)
# require(edgeR)
# require(ggbiplot)
# require(pheatmap)
# require(heatmaply)
# require(ggsci)

Sys.setlocale('LC_ALL','C')
type.list <- search_then_determine(path = input)
theme_set(cowplot::theme_cowplot())

Home

Column

Overview {.no-padding}

knitr::include_graphics('Figure_1_A.svg')

Column

Introduction

Thanks for using LncPipe! LncPipe is a nextflow-based pipeline which collected several popular tools and in-house scripts so as to identifiy and analyze novel lncRNAs from raw RNA sequencing data. As we known, identification and profiling of lncRNAs is a fundamental step to advance our knowledge on their function and regulatory mechanisms. We thus present LncPipe to facilitate related researches on LncRNA discovery and annotation. This report summarized all output files produced by LncPipe and shown them via fancy interactive plots. Briefly, LncPipeReporter's output contains sections as follows:

1. Reads alignment summary. Since the intial mapping step was carried out by STAR,HISAT or Tophat, we summarized the basic information that can help assess the library quanlity of each samples.
2. Novel lncRNAs basic charateristics such as length distribution, exon number, coding potiential and classifications.
3. Differential expression analysis if design file was feed. Abandunce of each transcript was assessed and compared between conditions. We also provide exploration analysis like PCA and clustering to help check the experimental design.
4. Comparison analysis was also performed between coding gene, known lncRNA and novel lncRNAs

Quick start

To run in R session:

library(LncPipeReporter)
run_reporter(input = system.file(file.path("extdata", "demo_results"),package = "LncPipeReporter"),
             output = 'reporter.html',
             theme = 'npg',
             cdf.percent = 10,
             max.lncrna.len = 10000,
             min.expressed.sample = 50,
             ask = FALSE)

To call in command line:

```{bash cmd_example, eval = FALSE, echo = TRUE}

Replace ... with arguments

$ Rscript -e "library(LncPipeReporter); run_reporter(input = '.', ...)"

> `...` stands for other arguments. You should use **single-quotes** here.

Column
-------------------------------------

### Authors 

[Qi Zhao](https://github.com/likelet) from Sun Yat-sen University Cancer Center, Guangdong, China.

[Yu Sun](http://icannotendure.space/) ([GitHub](https://github.com/bioinformatist)) from Nankai University. Tianjin, China.

### Release notes

This is the first version of LncPipeReporter.

### Citations

Coming soon!

### Lisence

[GPL v3 License](../LICENSE)

<!-- Quality control -->
<!-- ===================================== -->
<!-- Column -->
<!-- ------------------------------------- -->

<!-- ### Quality control -->
<!-- This sections presents the fastQC output. Samples were listed in the following table.  -->

<!-- ```r -->
<!-- fastqc.dirs <- list.files(path = params$input, pattern = '_fastqc$') -->
<!-- fastqc.table <- data.table::data.table( -->
<!-- Sample = substr(fastqc.dirs, 7, nchar(fastqc.dirs) - 7), -->
<!-- HTML = paste0( -->
<!-- '[**Click** for HTML](', -->
<!-- paste0(params$input, .Platform$file.sep, fastqc.dirs, .Platform$file.sep, "fastqc_report.html"), -->
<!-- ')' -->
<!-- ), -->
<!-- Summary = paste0( -->
<!-- '[**Click** for summmary](', -->
<!-- paste0(params$input, .Platform$file.sep, fastqc.dirs, .Platform$file.sep, "summary.txt"), -->
<!-- ')' -->
<!-- ), -->
<!-- Details = paste0( -->
<!-- '[**Click** for details](', -->
<!-- paste0(params$input, .Platform$file.sep, fastqc.dirs, .Platform$file.sep, "fastqc_data.txt"), -->
<!-- ')' -->
<!-- ) -->
<!-- ) -->
<!-- knitr::kable(fastqc.table, caption = "FastQC results") -->
<!-- rm(fastqc.dirs, fastqc.table) -->
<!-- invisible(gc()) -->
<!-- ``` -->


<!-- ```r -->
<!-- ``` -->

<!-- ```r -->
<!-- ``` -->

<!-- ```r -->
<!-- ``` -->

```r

bioinformatist/LncPipe-Reporter documentation built on May 6, 2019, 2:30 p.m.