vignettes/regionReport.Rmd
regionReport.RmdR is an open-source statistical environment which can be
easily modified to enhance its functionality via packages. regionReport
is a R package available via the Bioconductor
repository for packages. R can be installed on any
operating system from CRAN
after which you can install regionReport
by using the following commands in your R session:
if (!requireNamespace("BiocManager", quietly = TRUE)) {
install.packages("BiocManager")
}
BiocManager::install("regionReport")
## Check that you have a valid Bioconductor installation
BiocManager::valid()regionReport is based on many other packages and in particular in those that have implemented the infrastructure needed for dealing with RNA-seq data. That is, packages like GenomicFeatures that allow you to import the data and DESeq2 for generating differential expression results. A regionReport user is not expected to deal with those packages directly.
If you are asking yourself the question “Where do I start using Bioconductor?” you might be interested in this blog post.
As package developers, we try to explain clearly how to use our
packages and in which order to use the functions. But R and
Bioconductor have a steep learning curve so it is critical
to learn where to ask for help. The blog post quoted above mentions some
but we would like to highlight the Bioconductor support site
as the main resource for getting help: remember to use the
regionReport tag and check the older
posts. Other alternatives are available such as creating GitHub
issues and tweeting. However, please note that if you want to receive
help you should adhere to the posting
guidelines. It is particularly critical that you provide a small
reproducible example and your session information so package developers
can track down the source of the error.
We hope that regionReport will be useful for your research. Please use the following information to cite the package and the overall approach. Thank you!
## Citation info
citation("regionReport")## To cite package 'regionReport' in publications use:
##
## Collado-Torres L, Jaffe AE, Leek JT (2016). "regionReport:
## Interactive reports for region-level and feature-level genomic
## analyses [version2; referees: 2 approved, 1 approved with
## reservations]." _F1000Research_, *4*, 105.
## doi:10.12688/f1000research.6379.2
## <https://doi.org/10.12688/f1000research.6379.2>,
## <http://f1000research.com/articles/4-105/v2>.
##
## Collado-Torres L, Nellore A, Frazee AC, Wilks C, Love MI, Langmead B,
## Irizarry RA, Leek JT, Jaffe AE (2017). "Flexible expressed region
## analysis for RNA-seq with derfinder." _Nucl. Acids Res._.
## doi:10.1093/nar/gkw852 <https://doi.org/10.1093/nar/gkw852>,
## <http://nar.oxfordjournals.org/content/early/2016/09/29/nar.gkw852>.
##
## Collado-Torres L, Jaffe AE, Leek JT (2017). _regionReport: Generate
## HTML or PDF reports for a set of genomic regions or DESeq2/edgeR
## results_. doi:10.18129/B9.bioc.regionReport
## <https://doi.org/10.18129/B9.bioc.regionReport>,
## https://github.com/leekgroup/regionReport - R package version 1.41.0,
## <http://www.bioconductor.org/packages/regionReport>.
##
## To see these entries in BibTeX format, use 'print(<citation>,
## bibtex=TRUE)', 'toBibtex(.)', or set
## 'options(citation.bibtex.max=999)'.regionReport (Collado-Torres, Jaffe, and Leek, 2016) creates HTML or PDF reports for a set of genomic regions such as derfinder (Collado-Torres, Nellore, Frazee, Wilks, Love, Langmead, Irizarry, Leek, and Jaffe, 2017) results or for feature-level analyses performed with DESeq2 (Love, Huber, and Anders, 2014) or edgeR (Chen, Chen, Lun, Baldoni, and Smyth, 2024; Chen, Lun, and Smyth, 2016; Chen, Lun, and Smyth, 2014). The HTML reports are styled with rmarkdown (Allaire, Xie, Dervieux, McPherson, Luraschi, Ushey, Atkins, Wickham, Cheng, Chang, and Iannone, 2024) by default but can optionally be styled with knitrBootstrap (Hester, 2024).
This package includes a basic exploration for a general set of
genomic regions which can be easily customized to include the
appropriate conclusions and/or further exploration of the results. Such
a report can be generated using renderReport(). regionReport
has a separate template for running a basic exploration analysis of
derfinder
results by using derfinderReport().
derfinderReport() is specific to single base-level approach
derfinder
results. A third template is included for exploring DESeq2 or
edgeR
differential expression results.
All reports are written in R
Markdown format and include all the code for making the plots and
explorations in the report itself. For all templates, regionReport
relies on knitr (Xie,
2014), rmarkdown
(Allaire, Xie, Dervieux et al., 2024), DT (Xie, Cheng, and
Tan, 2024) and optionally knitrBootstrap
(Hester, 2024) for generating the report. The reports can be either in
HTML or PDF format and can be easily customized.
The plots in regionReport for exploring DESeq2 are powered by ggplot2 (Wickham, 2016) and pheatmap (Kolde, 2019).
The regionReport supplementary website regionReportSupp has examples of using regionReport with DESeq2 results. In particular, please look at DESeq2.html which has the code for generating some DESeq2 results based on the DESeq2 vignette. Then it uses those results to create HTML and PDF versions of the same report. The resulting reports are available in the following locations:
Note that in both examples we changed the ggplot2 theme
to theme_bw(). Also, in the PDF version we used the option
device = 'pdf' instead of the default
device = 'png' in DESeq2Report() since PDF
figures are more appropriate for PDF reports: they look better than PNG
figures.
If you want to create a similar HTML report as the one linked in this
section, simply run
example('DESeq2Report', 'regionReport', ask=FALSE). The
only difference will be the ggplot2 theme
for the plots.
regionReport
has the edgeReport() function that takes as input a
DGEList object and the results from the differential
expression analysis using edgeR.
edgeReport() internally uses DEFormats
to convert the results to DESeq2’s
format and then uses DESeqReport() to create the final
report. The report looks nearly the same whether you performed the
differential expression analysis with DESeq2 or
edgeR in
order to make more homogenous the exploratory data analysis step.
The regionReport supplementary website regionReportSupp has examples of using regionReport with edgeR results. In particular, please look at edgeR.html which has the code for generating some random data with DEFormats and performing the differential expression analysis with edgeR. Then it uses those results to create HTML and PDF versions of the same report. The resulting reports are available in the following locations:
Note that in both examples we changed the ggplot2 theme
to theme_linedraw(). Also, in the PDF version we used the
option device = 'pdf' instead of the default
device = 'png' in edgeReport() since PDF
figures are more appropriate for PDF reports: they look better than PNG
figures.
If you want to create a similar HTML report as the one linked in this
section, simply run
example('edgeReport', 'regionReport', ask=FALSE). The only
difference will be the ggplot2 theme
for the plots and the amount of data simulated with DEFormats.
The plots in regionReport for region reports are powered by derfinderPlot (Collado-Torres, Jaffe, and Leek, 2017), ggbio (Yin, Cook, and Lawrence, 2012), and ggplot2 (Wickham, 2016).
The regionReport
supplementary website regionReportSupp
has examples of using regionReport
with results from DiffBind
and derfinder.
Included as a vignette, this package also has an example using a small
data set derived from bumphunter.
These represent different uses of regionReport
for results from ChIP-seq, methylation, and RNA-seq data. In particular,
the DiffBind
example illustrates how to expand a basic report created with
renderReport().
For a general use case, you first have to identify a set of genomic
regions of interest and store it as a GRanges object. In a
typical workflow you will have some variables measured for each of the
regions, such as p-values and scores. renderReport() uses
the set of regions and three main arguments:
pvalueVars: this is a character vector (named
optionally) with the names of the variables that are bound between 0 and
1, such as p-values. For each of these variables,
renderReport() explores the distribution by chromosome, the
overall distribution, and makes a table with commonly used cutoffs.densityVars: is another character vector (named
optionally) with another set of variables you wish to explore by making
density graphs. This is commonly used for scores and other similar
numerical variables.significantVar: is a logical vector separating the
regions into by whether they are statistically significant. For example,
this information is used to explore the width of all the regions and
compare it the significant ones.Other parameters control the name of the report, where it’ll be located, the transcripts database used to annotate the nearest genes, graphical parameters, etc.
Here is a short example of how to use renderReport().
Note that we are using regions produced by derfinder
just for convenience sake.
## Load derfinder
library("derfinder")
regions <- genomeRegions$regions
## Assign chr length
library("GenomicRanges")
seqlengths(regions) <- c("chr21" = 48129895)
## The output will be saved in the 'derfinderReport-example' directory
dir.create("renderReport-example", showWarnings = FALSE, recursive = TRUE)
## Generate the HTML report
report <- renderReport(regions, "Example run",
pvalueVars = c(
"Q-values" = "qvalues", "P-values" = "pvalues"
), densityVars = c(
"Area" = "area", "Mean coverage" = "meanCoverage"
),
significantVar = regions$qvalues <= 0.05, nBestRegions = 20,
outdir = "renderReport-example"
)See the report created by this example here.
For derfinder
results created via the expressed regions-level approach you can use
renderReport() to explore the results. If you use DESeq2 to
perform the differential expression analysis of the expressed regions
you can then use DESeq2Report().
Prior to using regionReport::derfinderReport() you must
use derfinder
to analyze a specific data set. While there are many ways to do so, we
recommend using analyzeChr() with the same
prefix argument. Then merging the results with
mergeResults(). This is the recommended pipeline for
the single base-level approach.
Below, we run derfinder for the example data included in the package. The steps are:
## Load derfinder
library("derfinder")
## The output will be saved in the "derfinderReport-example" directory
dir.create("derfinderReport-example", showWarnings = FALSE, recursive = TRUE)The following code runs derfinder.
## Save the current path
initialPath <- getwd()
setwd(file.path(initialPath, "derfinderReport-example"))
## Generate output from derfinder
## Collapse the coverage information
collapsedFull <- collapseFullCoverage(list(genomeData$coverage),
verbose = TRUE
)
## Calculate library size adjustments
sampleDepths <- sampleDepth(collapsedFull,
probs = c(0.5), nonzero = TRUE,
verbose = TRUE
)
## Build the models
group <- genomeInfo$pop
adjustvars <- data.frame(genomeInfo$gender)
models <- makeModels(sampleDepths, testvars = group, adjustvars = adjustvars)
## Analyze chromosome 21
analysis <- analyzeChr(
chr = "21", coverageInfo = genomeData, models = models,
cutoffFstat = 1, cutoffType = "manual", seeds = 20140330, groupInfo = group,
mc.cores = 1, writeOutput = TRUE, returnOutput = TRUE
)
## Save the stats options for later
optionsStats <- analysis$optionsStats
## Change the directory back to the original one
setwd(initialPath)For convenience, we have included the derfinder results as part of regionReport. Note that the above functions are routinely checked as part of derfinder.
## Copy previous results
file.copy(system.file(file.path("extdata", "chr21"),
package = "derfinder",
mustWork = TRUE
), "derfinderReport-example", recursive = TRUE)## [1] TRUENext, proceed to merging the results.
## Merge the results from the different chromosomes. In this case, there's
## only one: chr21
mergeResults(
chrs = "chr21", prefix = "derfinderReport-example",
genomicState = genomicState$fullGenome
)## 2024-12-12 21:49:50.872773 mergeResults: Saving options used## 2024-12-12 21:49:50.873816 Loading chromosome chr21## Neither 'cutoffFstatUsed' nor 'optionsStats' were supplied, so the FWER calculation step will be skipped.## 2024-12-12 21:49:50.938771 mergeResults: Saving fullNullSummary## 2024-12-12 21:49:50.939364 mergeResults: Re-calculating the p-values## 2024-12-12 21:49:51.019439 mergeResults: Saving fullRegions## 2024-12-12 21:49:51.020447 mergeResults: assigning genomic states## 2024-12-12 21:49:51.115923 annotateRegions: counting## 2024-12-12 21:49:51.18319 annotateRegions: annotating## 2024-12-12 21:49:51.206474 mergeResults: Saving fullAnnotatedRegions## 2024-12-12 21:49:51.20737 mergeResults: Saving fullFstats## 2024-12-12 21:49:51.208124 mergeResults: Saving fullTime
## Load optionsStats
load(file.path("derfinderReport-example", "chr21", "optionsStats.Rdata"), verbose = TRUE)## Loading objects:
## optionsStatsOnce the derfinder output has been generated and merged, use derfinderReport() to create the HTML report.
## Load derfindeReport
library("regionReport")
## Generate the HTML report
report <- derfinderReport(
prefix = "derfinderReport-example", browse = FALSE,
nBestRegions = 15, makeBestClusters = TRUE,
fullCov = list("21" = genomeDataRaw$coverage), optionsStats = optionsStats
)Once the output is generated, you can browse the report from
R using browseURL() as shown below.
## Browse the report
browseURL(report)You can view a pre-compiled version of this report here.
Note that the reports require an active Internet connection to render correctly.
The report is self-explanatory and will change some of the text depending on the input options.
If the report is taking too long to compile (say more than 3 hours),
you might want to consider setting nBestCluters to a small
number or even set makeBestClusters to FALSE.
This package was made possible thanks to:
Code for creating the vignette
## Create the vignette
library("rmarkdown")
system.time(render("regionReport.Rmd", "BiocStyle::html_document"))
## Extract the R code
library("knitr")
knit("regionReport.Rmd", tangle = TRUE)
## Clean up
unlink("derfinderReport-example", recursive = TRUE)Date the vignette was generated.
## [1] "2024-12-12 21:49:52 UTC"Wallclock time spent generating the vignette.
## Time difference of 11.344 secsR session information.
## ─ Session info ───────────────────────────────────────────────────────────────────────────────────────────────────────
## setting value
## version R version 4.4.2 (2024-10-31)
## os Ubuntu 24.04.1 LTS
## system x86_64, linux-gnu
## ui X11
## language en
## collate en_US.UTF-8
## ctype en_US.UTF-8
## tz UTC
## date 2024-12-12
## pandoc 3.5 @ /usr/bin/ (via rmarkdown)
##
## ─ Packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────
## package * version date (UTC) lib source
## abind 1.4-8 2024-09-12 [1] RSPM (R 4.4.0)
## AnnotationDbi 1.68.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## backports 1.5.0 2024-05-23 [1] RSPM (R 4.4.0)
## base64enc 0.1-3 2015-07-28 [2] RSPM (R 4.4.0)
## bibtex 0.5.1 2023-01-26 [1] RSPM (R 4.4.0)
## Biobase 2.66.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## BiocGenerics 0.52.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## BiocIO 1.16.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## BiocManager 1.30.25 2024-08-28 [2] CRAN (R 4.4.2)
## BiocParallel 1.40.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## BiocStyle * 2.34.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## Biostrings 2.74.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## bit 4.5.0.1 2024-12-03 [1] RSPM (R 4.4.0)
## bit64 4.5.2 2024-09-22 [1] RSPM (R 4.4.0)
## bitops 1.0-9 2024-10-03 [1] RSPM (R 4.4.0)
## blob 1.2.4 2023-03-17 [1] RSPM (R 4.4.0)
## bookdown 0.41 2024-10-16 [1] RSPM (R 4.4.0)
## BSgenome 1.74.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## bslib 0.8.0 2024-07-29 [2] RSPM (R 4.4.0)
## bumphunter 1.48.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## cachem 1.1.0 2024-05-16 [2] RSPM (R 4.4.0)
## checkmate 2.3.2 2024-07-29 [1] RSPM (R 4.4.0)
## cli 3.6.3 2024-06-21 [2] RSPM (R 4.4.0)
## cluster 2.1.8 2024-12-11 [3] RSPM (R 4.4.0)
## codetools 0.2-20 2024-03-31 [3] CRAN (R 4.4.2)
## colorspace 2.1-1 2024-07-26 [1] RSPM (R 4.4.0)
## crayon 1.5.3 2024-06-20 [2] RSPM (R 4.4.0)
## curl 6.0.1 2024-11-14 [2] RSPM (R 4.4.0)
## data.table 1.16.4 2024-12-06 [1] RSPM (R 4.4.0)
## DBI 1.2.3 2024-06-02 [1] RSPM (R 4.4.0)
## DEFormats 1.34.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## DelayedArray 0.32.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## derfinder * 1.40.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## derfinderHelper 1.40.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## desc 1.4.3 2023-12-10 [2] RSPM (R 4.4.0)
## DESeq2 1.46.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## digest 0.6.37 2024-08-19 [2] RSPM (R 4.4.0)
## doRNG 1.8.6 2023-01-16 [1] RSPM (R 4.4.0)
## dplyr 1.1.4 2023-11-17 [1] RSPM (R 4.4.0)
## edgeR 4.4.1 2024-12-02 [1] Bioconductor 3.20 (R 4.4.2)
## evaluate 1.0.1 2024-10-10 [2] RSPM (R 4.4.0)
## fansi 1.0.6 2023-12-08 [2] RSPM (R 4.4.0)
## fastmap 1.2.0 2024-05-15 [2] RSPM (R 4.4.0)
## foreach 1.5.2 2022-02-02 [1] RSPM (R 4.4.0)
## foreign 0.8-87 2024-06-26 [3] CRAN (R 4.4.2)
## Formula 1.2-5 2023-02-24 [1] RSPM (R 4.4.0)
## fs 1.6.5 2024-10-30 [2] RSPM (R 4.4.0)
## generics 0.1.3 2022-07-05 [1] RSPM (R 4.4.0)
## GenomeInfoDb 1.42.1 2024-11-28 [1] Bioconductor 3.20 (R 4.4.2)
## GenomeInfoDbData 1.2.13 2024-12-10 [1] Bioconductor
## GenomicAlignments 1.42.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## GenomicFeatures 1.58.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## GenomicFiles 1.42.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## GenomicRanges 1.58.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## ggplot2 3.5.1 2024-04-23 [1] RSPM (R 4.4.0)
## glue 1.8.0 2024-09-30 [2] RSPM (R 4.4.0)
## gridExtra 2.3 2017-09-09 [1] RSPM (R 4.4.0)
## gtable 0.3.6 2024-10-25 [1] RSPM (R 4.4.0)
## Hmisc 5.2-1 2024-12-02 [1] RSPM (R 4.4.0)
## htmlTable 2.4.3 2024-07-21 [1] RSPM (R 4.4.0)
## htmltools 0.5.8.1 2024-04-04 [2] RSPM (R 4.4.0)
## htmlwidgets 1.6.4 2023-12-06 [2] RSPM (R 4.4.0)
## httr 1.4.7 2023-08-15 [1] RSPM (R 4.4.0)
## IRanges 2.40.1 2024-12-05 [1] Bioconductor 3.20 (R 4.4.2)
## iterators 1.0.14 2022-02-05 [1] RSPM (R 4.4.0)
## jquerylib 0.1.4 2021-04-26 [2] RSPM (R 4.4.0)
## jsonlite 1.8.9 2024-09-20 [2] RSPM (R 4.4.0)
## KEGGREST 1.46.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## knitr 1.49 2024-11-08 [2] RSPM (R 4.4.0)
## knitrBootstrap 1.0.3 2024-02-06 [1] RSPM (R 4.4.0)
## lattice 0.22-6 2024-03-20 [3] CRAN (R 4.4.2)
## lifecycle 1.0.4 2023-11-07 [2] RSPM (R 4.4.0)
## limma 3.62.1 2024-11-03 [1] Bioconductor 3.20 (R 4.4.2)
## locfit 1.5-9.10 2024-06-24 [1] RSPM (R 4.4.0)
## lubridate 1.9.4 2024-12-08 [1] RSPM (R 4.4.0)
## magrittr 2.0.3 2022-03-30 [2] RSPM (R 4.4.0)
## markdown 1.13 2024-06-04 [1] RSPM (R 4.4.0)
## Matrix 1.7-1 2024-10-18 [3] CRAN (R 4.4.2)
## MatrixGenerics 1.18.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## matrixStats 1.4.1 2024-09-08 [1] RSPM (R 4.4.0)
## memoise 2.0.1 2021-11-26 [2] RSPM (R 4.4.0)
## munsell 0.5.1 2024-04-01 [1] RSPM (R 4.4.0)
## nnet 7.3-19 2023-05-03 [3] CRAN (R 4.4.2)
## pillar 1.9.0 2023-03-22 [2] RSPM (R 4.4.0)
## pkgconfig 2.0.3 2019-09-22 [2] RSPM (R 4.4.0)
## pkgdown 2.1.1 2024-09-17 [2] RSPM (R 4.4.0)
## plyr 1.8.9 2023-10-02 [1] RSPM (R 4.4.0)
## png 0.1-8 2022-11-29 [1] RSPM (R 4.4.0)
## qvalue 2.38.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## R6 2.5.1 2021-08-19 [2] RSPM (R 4.4.0)
## ragg 1.3.3 2024-09-11 [2] RSPM (R 4.4.0)
## Rcpp 1.0.13-1 2024-11-02 [2] RSPM (R 4.4.0)
## RCurl 1.98-1.16 2024-07-11 [1] RSPM (R 4.4.0)
## RefManageR * 1.4.0 2022-09-30 [1] RSPM (R 4.4.0)
## regionReport * 1.41.0 2024-12-12 [1] Bioconductor
## reshape2 1.4.4 2020-04-09 [1] RSPM (R 4.4.0)
## restfulr 0.0.15 2022-06-16 [1] RSPM (R 4.4.0)
## rjson 0.2.23 2024-09-16 [1] RSPM (R 4.4.0)
## rlang 1.1.4 2024-06-04 [2] RSPM (R 4.4.0)
## rmarkdown 2.29 2024-11-04 [2] RSPM (R 4.4.0)
## rngtools 1.5.2 2021-09-20 [1] RSPM (R 4.4.0)
## rpart 4.1.23 2023-12-05 [3] CRAN (R 4.4.2)
## Rsamtools 2.22.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## RSQLite 2.3.9 2024-12-03 [1] RSPM (R 4.4.0)
## rstudioapi 0.17.1 2024-10-22 [2] RSPM (R 4.4.0)
## rtracklayer 1.66.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## S4Arrays 1.6.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## S4Vectors 0.44.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## sass 0.4.9 2024-03-15 [2] RSPM (R 4.4.0)
## scales 1.3.0 2023-11-28 [1] RSPM (R 4.4.0)
## sessioninfo * 1.2.2 2021-12-06 [2] RSPM (R 4.4.0)
## SparseArray 1.6.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## statmod 1.5.0 2023-01-06 [1] RSPM (R 4.4.0)
## stringi 1.8.4 2024-05-06 [2] RSPM (R 4.4.0)
## stringr 1.5.1 2023-11-14 [2] RSPM (R 4.4.0)
## SummarizedExperiment 1.36.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## systemfonts 1.1.0 2024-05-15 [2] RSPM (R 4.4.0)
## textshaping 0.4.1 2024-12-06 [2] RSPM (R 4.4.0)
## tibble 3.2.1 2023-03-20 [2] RSPM (R 4.4.0)
## tidyselect 1.2.1 2024-03-11 [1] RSPM (R 4.4.0)
## timechange 0.3.0 2024-01-18 [1] RSPM (R 4.4.0)
## UCSC.utils 1.2.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## utf8 1.2.4 2023-10-22 [2] RSPM (R 4.4.0)
## VariantAnnotation 1.52.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## vctrs 0.6.5 2023-12-01 [2] RSPM (R 4.4.0)
## xfun 0.49 2024-10-31 [2] RSPM (R 4.4.0)
## XML 3.99-0.17 2024-06-25 [1] RSPM (R 4.4.0)
## xml2 1.3.6 2023-12-04 [2] RSPM (R 4.4.0)
## XVector 0.46.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
## yaml 2.3.10 2024-07-26 [2] RSPM (R 4.4.0)
## zlibbioc 1.52.0 2024-10-29 [1] Bioconductor 3.20 (R 4.4.2)
##
## [1] /__w/_temp/Library
## [2] /usr/local/lib/R/site-library
## [3] /usr/local/lib/R/library
##
## ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────This vignette was generated using BiocStyle (Oleś, 2024) with knitr (Xie, 2014) and rmarkdown (Allaire, Xie, Dervieux et al., 2024) running behind the scenes.
Citations made with RefManageR (McLean, 2017).
[1] J. Allaire, Y. Xie, C. Dervieux, et al. rmarkdown: Dynamic Documents for R. R package version 2.29. 2024. URL: https://github.com/rstudio/rmarkdown.
[2] S. Arora, M. Morgan, M. Carlson, et al. GenomeInfoDb: Utilities for manipulating chromosome and other ‘seqname’ identifiers. 2017. DOI: 10.18129/B9.bioc.GenomeInfoDb.
[3] B. Auguie. gridExtra: Miscellaneous Functions for "Grid" Graphics. R package version 2.3. 2017.
[4] M. Carlson and B. P. Maintainer. TxDb.Hsapiens.UCSC.hg19.knownGene: Annotation package for TxDb object(s). R package version 3.2.2. 2015.
[5] Y. Chen, L. Chen, A. T. L. Lun, et al. “edgeR 4.0: powerful differential analysis of sequencing data with expanded functionality and improved support for small counts and larger datasets”. In: bioRxiv (2024). DOI: 10.1101/2024.01.21.576131.
[6] Y. Chen, A. T. L. Lun, and G. K. Smyth. “Differential expression analysis of complex RNA-seq experiments using edgeR”. In: Statistical Analysis of Next Generation Sequencing Data. Ed. by S. Datta and D. Nettleton. New York: Springer, 2014, pp. 51-74.
[7] Y. Chen, A. T. L. Lun, and G. K. Smyth. “From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline”. In: F1000Research 5 (2016), p. 1438. DOI: 10.12688/f1000research.8987.2.
[8] L. Collado-Torres, A. E. Jaffe, and J. T. Leek. derfinderPlot: Plotting functions for derfinder. https://github.com/leekgroup/derfinderPlot - R package version 1.40.0. 2017. DOI: 10.18129/B9.bioc.derfinderPlot. URL: http://www.bioconductor.org/packages/derfinderPlot.
[9] L. Collado-Torres, A. E. Jaffe, and J. T. Leek. “regionReport: Interactive reports for region-level and feature-level genomic analyses [version2; referees: 2 approved, 1 approved with reservations]”. In: F1000Research 4 (2016), p. 105. DOI: 10.12688/f1000research.6379.2. URL: http://f1000research.com/articles/4-105/v2.
[10] L. Collado-Torres, A. Nellore, A. C. Frazee, et al. “Flexible expressed region analysis for RNA-seq with derfinder”. In: Nucl. Acids Res. (2017). DOI: 10.1093/nar/gkw852. URL: http://nar.oxfordjournals.org/content/early/2016/09/29/nar.gkw852.
[11] E. de Jonge. whisker: mustache for R, Logicless Templating. R package version 0.4.1. 2022. URL: https://github.com/edwindj/whisker.
[12] J. Hester. knitrBootstrap: ‘knitr’ Bootstrap Framework. R package version 1.0.3. 2024. URL: https://github.com/jimhester/knitrBootstrap#readme.
[13] A. E. Jaffe, P. Murakami, H. Lee, et al. “Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies”. In: International journal of epidemiology 41.1 (2012), pp. 200–209. DOI: 10.1093/ije/dyr238.
[14] R. Kolde. pheatmap: Pretty Heatmaps. R package version 1.0.12. 2019.
[15] M. Lawrence, W. Huber, H. Pagès, et al. “Software for Computing and Annotating Genomic Ranges”. In: PLoS Computational Biology 9 (8 2013). DOI: 10.1371/journal.pcbi.1003118. URL: http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1003118}.
[16] M. I. Love, W. Huber, and S. Anders. “Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2”. In: Genome Biology 15 (12 2014), p. 550. DOI: 10.1186/s13059-014-0550-8.
[17] M. W. McLean. “RefManageR: Import and Manage BibTeX and BibLaTeX References in R”. In: The Journal of Open Source Software (2017). DOI: 10.21105/joss.00338.
[18] E. Neuwirth. RColorBrewer: ColorBrewer Palettes. R package version 1.1-3. 2022.
[19] A. Oleś. BiocStyle: Standard styles for vignettes and other Bioconductor documents. R package version 2.34.0. 2024. DOI: 10.18129/B9.bioc.BiocStyle. URL: https://bioconductor.org/packages/BiocStyle.
[20] A. Oleś. DEFormats: Differential gene expression data formats converter. R package version 1.34.0. 2024. DOI: 10.18129/B9.bioc.DEFormats. URL: https://bioconductor.org/packages/DEFormats.
[21] R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria, 2024. URL: https://www.R-project.org/.
[22] H. Wickham. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York, 2016. ISBN: 978-3-319-24277-4. URL: https://ggplot2.tidyverse.org.
[23] H. Wickham, W. Chang, R. Flight, et al. sessioninfo: R Session Information. R package version 1.2.2, https://r-lib.github.io/sessioninfo/. 2021. URL: https://github.com/r-lib/sessioninfo#readme.
[24] S. N. Wood. “Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models”. In: Journal of the Royal Statistical Society (B) 73.1 (2011), pp. 3-36.
[25] Y. Xie. “knitr: A Comprehensive Tool for Reproducible Research in R”. In: Implementing Reproducible Computational Research. Ed. by V. Stodden, F. Leisch and R. D. Peng. ISBN 978-1466561595. Chapman and Hall/CRC, 2014.
[26] Y. Xie, J. Cheng, and X. Tan. DT: A Wrapper of the JavaScript Library ‘DataTables’. R package version 0.33. 2024. URL: https://github.com/rstudio/DT.
[27] T. Yin, D. Cook, and M. Lawrence. “ggbio: an R package for extending the grammar of graphics for genomic data”. In: Genome Biology 13.8 (2012), p. R77.
[28] T. Yin, M. Lawrence, and D. Cook. biovizBase: Basic graphic utilities for visualization of genomic data. R package version 1.54.0. 2024. DOI: 10.18129/B9.bioc.biovizBase. URL: https://bioconductor.org/packages/biovizBase.