Contents

1 Description

This vignettes will guide you throught a tipycal usage of the easyreporting package, while performing a simplified bioinformatics analysis workflow.

2 Requirements

For the usage you just need to load the easyreporting package, which will load the rmarkdown and tools packages.

if(!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("easyreporting")
library("easyreporting")

3 easyreporting instance creation

For simplicity we setup a project directory path starting from the working directory for our report, but you can just enter any path. The filenamepath and the title parameters are mandatory, while the author paramenter is optional.

Once created the easyreporting class instance, we can use it into our further code to make other operations. It stores some variables for us, in order to not be called again during next opreations. For example the name and the path of the report, the type of report and the general rmarkdown options of the document.

proj.path <- file.path(tempdir(), "bioinfo_report")
bioEr <- easyreporting(filenamePath=proj.path, title="bioinfo_report",
                        author=c("Dario Righelli"))

4 Loading Data

For this vignette we previously downloaded an xlsx table from GEO with accession number GSE134118.

For importing the xls file, we prepared an ad-hoc function called importData stored in the importFunctions.R file.

mkdTitle(bioEr, title="Loading Counts Data")
mkdCodeChunkComplete(object=bioEr, message=paste0("geneCounts <- importData(sy",
        "stem.file('extdata/GSE134118_Table_S3.xlsx',",
        "package='easyreporting'))"), 
        sourceFilesList=system.file("script/importFunctions.R", 
                                            package="easyreporting"))

5 Counts exploration

In order to explore the counts we can perform a PCA with plotPCA function within the EDASeq package. To trace this step we again use a mkdCodeChunkComplete function, inserting the function call as the message to track.

mkdTitle(bioEr, title="Plot PCA on count data", level=2)
mkdCodeChunkComplete(bioEr, 
                        message="EDASeq::plotPCA(as.matrix(geneCounts))")

6 Differential expression

Let’s suppose we stored a function for applying edgeR on our counts in an R file. We called the file geneFunctions.R and the function applyEdgeRExample.

It will be easy for us to source the file and to call the function with a single call of easyreporting package. By aid of the mkdCodeChunkCommented function we can add a comment preceeding the code chunk.

mkdTitle(bioEr, "Differential Expression Analysis")
mkdCodeChunkCommented(bioEr, 
    codeMsg="degList <- applyEdgeRExample(counts=geneCounts, 
                    samples=colnames(geneCounts), contrast='Pleura - Broth')", 
    commentMsg=paste0("As we saw from the PCA, the groups are well separated",
        ", so we can perform a Differential Expression analysis with edgeR."),
    sourceFilesList=system.file("script/geneFunctions.R", 
                                package="easyreporting"))

6.1 Inspecting DEGs

Usually to have a graphical representation of the obtained results we can plot an MD plot by using the plotMD function of the limma package.

Also in this case, we create a new title and a new chunk within the function to call.

mkdTitle(bioEr, "MD Plot of DEGs", level=2)
mkdCodeChunkComplete(bioEr, message="limma::plotMD(degList$test)")

7 Compiling the report

Once finished our analysis it is possible to compile the produced rmarkdown report simply by using the compile method. The compile method appends a sessionInfo() to the report to trace all the packages and versions used for the analysis.

compile(bioEr)
#> processing file: bioinfo_report.Rmd
#> output file: bioinfo_report.knit.md
#> /usr/bin/pandoc +RTS -K512m -RTS bioinfo_report.utf8.md --to html4 --from markdown+autolink_bare_uris+tex_math_single_backslash+smart --output bioinfo_report.html --email-obfuscation none --self-contained --standalone --section-divs --template /home/biocbuild/bbs-3.11-bioc/R/library/rmarkdown/rmd/h/default.html --no-highlight --variable highlightjs=1 --variable 'theme:bootstrap' --include-in-header /tmp/Rtmp8ahaoS/rmarkdown-str15f81a2ddaa9.html --mathjax --variable 'mathjax-url:https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML' --lua-filter /home/biocbuild/bbs-3.11-bioc/R/library/rmarkdown/rmd/lua/pagebreak.lua --lua-filter /home/biocbuild/bbs-3.11-bioc/R/library/rmarkdown/rmd/lua/latex-div.lua
#> 
#> Output created: bioinfo_report.html

8 Session Info

sessionInfo()
#> R version 4.0.0 (2020-04-24)
#> Platform: x86_64-pc-linux-gnu (64-bit)
#> Running under: Ubuntu 18.04.4 LTS
#> 
#> Matrix products: default
#> BLAS:   /home/biocbuild/bbs-3.11-bioc/R/lib/libRblas.so
#> LAPACK: /home/biocbuild/bbs-3.11-bioc/R/lib/libRlapack.so
#> 
#> locale:
#>  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
#>  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=C              
#>  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
#>  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
#>  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
#> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       
#> 
#> attached base packages:
#> [1] stats     graphics  grDevices utils     datasets  methods   base     
#> 
#> other attached packages:
#> [1] easyreporting_1.0.0 BiocStyle_2.16.0   
#> 
#> loaded via a namespace (and not attached):
#>  [1] Biobase_2.48.0              httr_1.4.1                 
#>  [3] edgeR_3.30.0                bit64_0.9-7                
#>  [5] splines_4.0.0               R.utils_2.9.2              
#>  [7] assertthat_0.2.1            statmod_1.4.34             
#>  [9] askpass_1.1                 BiocManager_1.30.10        
#> [11] stats4_4.0.0                aroma.light_3.18.0         
#> [13] BiocFileCache_1.12.0        latticeExtra_0.6-29        
#> [15] blob_1.2.1                  GenomeInfoDbData_1.2.3     
#> [17] cellranger_1.1.0            Rsamtools_2.4.0            
#> [19] yaml_2.2.1                  progress_1.2.2             
#> [21] pillar_1.4.3                RSQLite_2.2.0              
#> [23] lattice_0.20-41             limma_3.44.0               
#> [25] glue_1.4.0                  digest_0.6.25              
#> [27] GenomicRanges_1.40.0        RColorBrewer_1.1-2         
#> [29] XVector_0.28.0              htmltools_0.4.0            
#> [31] Matrix_1.2-18               R.oo_1.23.0                
#> [33] XML_3.99-0.3                pkgconfig_2.0.3            
#> [35] ShortRead_1.46.0            biomaRt_2.44.0             
#> [37] magick_2.3                  genefilter_1.70.0          
#> [39] bookdown_0.18               zlibbioc_1.34.0            
#> [41] purrr_0.3.4                 xtable_1.8-4               
#> [43] jpeg_0.1-8.1                BiocParallel_1.22.0        
#> [45] tibble_3.0.1                openssl_1.4.1              
#> [47] annotate_1.66.0             IRanges_2.22.0             
#> [49] ellipsis_0.3.0              SummarizedExperiment_1.18.0
#> [51] GenomicFeatures_1.40.0      BiocGenerics_0.34.0        
#> [53] survival_3.1-12             magrittr_1.5               
#> [55] crayon_1.3.4                readxl_1.3.1               
#> [57] memoise_1.1.0               evaluate_0.14              
#> [59] R.methodsS3_1.8.0           hwriter_1.3.2              
#> [61] tools_4.0.0                 prettyunits_1.1.1          
#> [63] hms_0.5.3                   lifecycle_0.2.0            
#> [65] matrixStats_0.56.0          stringr_1.4.0              
#> [67] S4Vectors_0.26.0            locfit_1.5-9.4             
#> [69] DelayedArray_0.14.0         AnnotationDbi_1.50.0       
#> [71] Biostrings_2.56.0           compiler_4.0.0             
#> [73] DESeq_1.40.0                GenomeInfoDb_1.24.0        
#> [75] rlang_0.4.5                 grid_4.0.0                 
#> [77] RCurl_1.98-1.2              rappdirs_0.3.1             
#> [79] bitops_1.0-6                rmarkdown_2.1              
#> [81] codetools_0.2-16            curl_4.3                   
#> [83] DBI_1.1.0                   EDASeq_2.22.0              
#> [85] R6_2.4.1                    GenomicAlignments_1.24.0   
#> [87] rtracklayer_1.48.0          knitr_1.28                 
#> [89] dplyr_0.8.5                 bit_1.1-15.2               
#> [91] stringi_1.4.6               parallel_4.0.0             
#> [93] Rcpp_1.0.4.6                png_0.1-7                  
#> [95] vctrs_0.2.4                 geneplotter_1.66.0         
#> [97] tidyselect_1.0.0            dbplyr_1.4.3               
#> [99] xfun_0.13