Project: Example bumphunter.
This report is meant to help explore a set of genomic regions and was generated using the regionReport
(Collado-Torres, Jaffe, and Leek, 2016) package. While the report is rich, it is meant to just start the exploration of the results and exemplify some of the code used to do so. If you need a more in-depth analysis for your specific data set you might want to use the customCode
argument.
Most plots were made with using ggplot2
(Wickham, 2016).
## knitrBoostrap and device chunk options
library('knitr')
opts_chunk$set(bootstrap.show.code = FALSE, dev = device, crop = NULL)
if(!outputIsHTML) opts_chunk$set(bootstrap.show.code = FALSE, dev = device, echo = FALSE)
#### Libraries needed
## Bioconductor
library('bumphunter')
library('derfinder')
library('derfinderPlot')
library('GenomeInfoDb')
library('GenomicRanges')
library('ggbio')
## Transcription database to use by default
if(is.null(txdb)) {
library('TxDb.Hsapiens.UCSC.hg19.knownGene')
txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene::TxDb.Hsapiens.UCSC.hg19.knownGene
}
## CRAN
library('ggplot2')
if(!is.null(theme)) theme_set(theme)
library('grid')
library('gridExtra')
library('knitr')
library('RColorBrewer')
library('mgcv')
library('whisker')
library('DT')
library('sessioninfo')
#### Code setup
## For ggplot
tmp <- regions
names(tmp) <- seq_len(length(tmp))
regions.df <- as.data.frame(tmp)
regions.df$width <- width(tmp)
rm(tmp)
## Special subsets: need at least 3 points for a density plot
keepChr <- table(regions.df$seqnames) > 2
regions.df.plot <- subset(regions.df, seqnames %in% names(keepChr[keepChr]))
if(hasSignificant) {
## Keep only those sig
regions.df.sig <- regions.df[significantVar, ]
keepChr <- table(regions.df.sig$seqnames) > 2
regions.df.sig <- subset(regions.df.sig, seqnames %in% names(keepChr[keepChr]))
}
## Find which chrs are present in the data set
chrs <- levels(seqnames(regions))
## areaVar initialize
areaVar <- NULL
p2a <- ggplot(regions.df.plot, aes(x=log10(width), colour=seqnames)) +
geom_line(stat='density') + labs(title='Density of region lengths') +
xlab('Region width (log10)') + scale_colour_discrete(limits=chrs) +
theme(legend.title=element_blank())
p2a
This plot shows the density of the region lengths for all regions.
for(i in seq_len(length(densityVars))) {
densityVarName <- names(densityVars[i])
densityVarName <- ifelse(is.null(densityVarName), densityVars[i], densityVarName)
cat(knit_child(text = whisker.render(templateDensityInUse, list(varName = densityVars[i], densityVarName = densityVarName)), quiet = TRUE), sep = '\n')
}
p3aarea <- ggplot(regions.df.plot[is.finite(regions.df.plot[, 'area']), ], aes(x=area, colour=seqnames)) +
geom_line(stat='density') + labs(title='Density of Area') +
xlab('Area') + scale_colour_discrete(limits=chrs) +
theme(legend.title=element_blank())
p3aarea
This plot shows the density of the Area for all regions.
p3avalue <- ggplot(regions.df.plot[is.finite(regions.df.plot[, 'value']), ], aes(x=value, colour=seqnames)) +
geom_line(stat='density') + labs(title='Density of Value') +
xlab('Value') + scale_colour_discrete(limits=chrs) +
theme(legend.title=element_blank())
p3avalue
This plot shows the density of the Value for all regions.
p3aclusterL <- ggplot(regions.df.plot[is.finite(regions.df.plot[, 'clusterL']), ], aes(x=clusterL, colour=seqnames)) +
geom_line(stat='density') + labs(title='Density of Cluster Length') +
xlab('Cluster Length') + scale_colour_discrete(limits=chrs) +
theme(legend.title=element_blank())
p3aclusterL
This plot shows the density of the Cluster Length for all regions.
The following plots were made using ggbio
(Yin, Cook, and Lawrence, 2012) which in turn uses ggplot2
(Wickham, 2016). For more details check plotOverview
in derfinderPlot
(Collado-Torres, Jaffe, and Leek, 2017).
This plot shows the genomic locations of the regions found in the analysis. The significant regions are highlighted and the of the regions is shown on top of each chromosome (skipped because there was no applicable variable).
## Annotate regions with bumphunter
if(is.null(annotation)) {
genes <- annotateTranscripts(txdb = txdb)
annotation <- matchGenes(x = regions, subject = genes)
}
## Warning in library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package
## called 'org.Hs.eg.db'
## Make the plot
plotOverview(regions=regions, annotation=annotation, type='annotation', base_size=overviewParams$base_size, areaRel=overviewParams$areaRel, legend.position=c(0.97, 0.12))
This genomic overview plot shows the annotation region type for the regions as determined using bumphunter
(Jaffe, Murakami, Lee et al., 2012). Note that the regions are shown only if the annotation information is available. Below is a table of the actual number of results per annotation region type.
annoReg <- table(annotation$region, useNA='always')
annoReg.df <- data.frame(Region=names(annoReg), Count=as.vector(annoReg))
if(outputIsHTML) {
kable(annoReg.df, format = 'markdown', align=rep('c', 3))
} else {
kable(annoReg.df)
}
Region | Count |
---|---|
upstream | 10 |
promoter | 0 |
overlaps 5’ | 0 |
inside | 0 |
overlaps 3’ | 0 |
close to 3’ | 0 |
downstream | 5 |
covers | 0 |
NA | 0 |
This genomic overview plot shows the annotation region type for the statistically significant regions. Note that the regions are shown only if the annotation information is available. Plot skipped because there are no significant regions.
Below is a table summarizing the number of genomic states per region as determined using derfinder
(Collado-Torres, Nellore, Frazee et al., 2017).
## Construct genomic state object
genomicState <- makeGenomicState(txdb = txdb, chrs = chrs, verbose = FALSE)
## 'select()' returned 1:1 mapping between keys and columns
## Annotate regions by genomic state
annotatedRegions <- annotateRegions(regions, genomicState$fullGenome, verbose = FALSE)
## Genomic states table
info <- do.call(rbind, lapply(annotatedRegions$countTable, function(x) { data.frame(table(x)) }))
colnames(info) <- c('Number of Overlapping States', 'Frequency')
info$State <- gsub('\\..*', '', rownames(info))
rownames(info) <- NULL
if(outputIsHTML) {
kable(info, format = 'markdown', align=rep('c', 4))
} else {
kable(info)
}
Number of Overlapping States | Frequency | State |
---|---|---|
0 | 15 | exon |
1 | 15 | intergenic |
0 | 15 | intron |
The following is a venn diagram showing how many regions overlap known exons, introns, and intergenic segments, none of them, or multiple of these groups.
## Venn diagram for all regions
venn <- vennRegions(annotatedRegions, counts.col = 'blue',
main = 'Regions overlapping genomic states')
Below is an interactive table with the top 15 regions (out of 15) as ranked by p-value without ranking because no p-value information was provided. Inf and -Inf are shown as 1e100 and -1e100 respectively. Use the search function to find your region of interest or sort by one of the columns.
## Add annotation information
regions.df <- cbind(regions.df, annotation)
## Rank by p-value (first pvalue variable supplied)
if(hasPvalueVars){
topRegions <- head(regions.df[order(regions.df[, pvalueVars[1]],
decreasing = FALSE), ], nBestRegions)
topRegions <- cbind(data.frame('pvalueRank' = seq_len(nrow(topRegions))),
topRegions)
} else {
topRegions <- head(regions.df, nBestRegions)
}
## Clean up -Inf, Inf if present
## More details at https://github.com/ramnathv/rCharts/issues/259
replaceInf <- function(df, colsubset=seq_len(ncol(df))) {
for(i in colsubset) {
inf.idx <- !is.finite(df[, i])
if(any(inf.idx)) {
inf.sign <- sign(df[inf.idx, i])
df[inf.idx, i] <- inf.sign * 1e100
}
}
return(df)
}
topRegions <- replaceInf(topRegions, which(sapply(topRegions, function(x) {
class(x) %in% c('numeric', 'integer')})))
## Make the table
greptext <- 'value$|area$|mean|log2FoldChange'
greppval <- 'pvalues$|qvalues$|fwer$'
if(hasPvalueVars) {
greppval <- paste0(paste(pvalueVars, collapse = '$|'), '$|', greppval)
}
if(hasDensityVars) {
greptext <- paste0(paste(densityVars, collapse = '$|'), '$|', greptext)
}
for(i in which(grepl(greppval, colnames(topRegions)))) topRegions[, i] <- format(topRegions[, i], scientific = TRUE)
if(outputIsHTML) {
datatable(topRegions, options = list(pagingType='full_numbers', pageLength=10, scrollX='100%'), rownames = FALSE) %>% formatRound(which(grepl(greptext, colnames(topRegions))), digits)
} else {
## Only print the top part if your output is a PDF file
df_top <- head(topRegions, 20)
for(i in which(grepl(greptext, colnames(topRegions)))) df_top[, i] <- round(df_top[, i], digits)
kable(df_top)
}
This report was generated in path /__w/regionReport/regionReport/docs/reference using the following call to renderReport()
:
## renderReport(regions = regions, project = "Example bumphunter",
## pvalueVars = NULL, densityVars = c(Area = "area", Value = "value",
## `Cluster Length` = "clusterL"), significantVar = NULL,
## outdir = "bumphunter-example", output = "bumphunter-example",
## device = "png")
Date the report was generated.
## [1] "2024-12-12 21:49:23 UTC"
Wallclock time spent generating the report.
## Time difference of 40.021 secs
R
session information.
## [1m[36m─ Session info ───────────────────────────────────────────────────────────────────────────────────────────────────────[39m[22m
## [3m[90msetting [39m[23m [3m[90mvalue[39m[23m
## 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 C
## ctype en_US.UTF-8
## tz UTC
## date 2024-12-12
## pandoc 3.5 @ /usr/bin/ (via rmarkdown)
##
## [1m[36m─ Packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────[39m[22m
## [3m[90mpackage [39m[23m [3m[90m*[39m[23m [3m[90mversion [39m[23m [3m[90mdate (UTC)[39m[23m [3m[90mlib[39m[23m [3m[90msource[39m[23m
## abind 1.4-8 [90m2024-09-12[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## annotate 1.84.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## AnnotationDbi * 1.68.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## AnnotationFilter 1.30.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## backports 1.5.0 [90m2024-05-23[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## base64enc 0.1-3 [90m2015-07-28[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## bibtex 0.5.1 [90m2023-01-26[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## Biobase * 2.66.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## BiocFileCache 2.14.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## BiocGenerics * 0.52.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## BiocIO 1.16.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## BiocManager 1.30.25 [90m2024-08-28[39m [90m[2][39m [90mCRAN (R 4.4.2)[39m
## BiocParallel * 1.40.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## BiocStyle * 2.34.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## biomaRt 2.62.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## Biostrings 2.74.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## biovizBase * 1.54.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## bit 4.5.0.1 [90m2024-12-03[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## bit64 4.5.2 [90m2024-09-22[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## bitops 1.0-9 [90m2024-10-03[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## blob 1.2.4 [90m2023-03-17[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## bookdown 0.41 [90m2024-10-16[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## BSgenome 1.74.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## bslib 0.8.0 [90m2024-07-29[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## bumphunter * 1.48.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## cachem 1.1.0 [90m2024-05-16[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## checkmate 2.3.2 [90m2024-07-29[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## cli 3.6.3 [90m2024-06-21[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## cluster 2.1.8 [90m2024-12-11[39m [90m[3][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## codetools 0.2-20 [90m2024-03-31[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## colorspace 2.1-1 [90m2024-07-26[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## crayon 1.5.3 [90m2024-06-20[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## crosstalk 1.2.1 [90m2023-11-23[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## curl 6.0.1 [90m2024-11-14[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## data.table 1.16.4 [90m2024-12-06[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## DBI 1.2.3 [90m2024-06-02[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## dbplyr 2.5.0 [90m2024-03-19[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## DEFormats * 1.34.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## DelayedArray 0.32.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## derfinder * 1.40.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## derfinderHelper 1.40.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## derfinderPlot * 1.40.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## desc 1.4.3 [90m2023-12-10[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## DESeq2 * 1.46.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## devtools 2.4.5 [90m2022-10-11[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## DEXSeq * 1.52.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## dichromat 2.0-0.1 [90m2022-05-02[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## digest 0.6.37 [90m2024-08-19[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## doRNG 1.8.6 [90m2023-01-16[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## downlit 0.4.4 [90m2024-06-10[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## dplyr 1.1.4 [90m2023-11-17[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## DT * 0.33 [90m2024-04-04[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## edgeR * 4.4.1 [90m2024-12-02[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## ellipsis 0.3.2 [90m2021-04-29[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## ensembldb 2.30.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## evaluate 1.0.1 [90m2024-10-10[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## fansi 1.0.6 [90m2023-12-08[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## farver 2.1.2 [90m2024-05-13[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## fastmap 1.2.0 [90m2024-05-15[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## filelock 1.0.3 [90m2023-12-11[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## foreach * 1.5.2 [90m2022-02-02[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## foreign 0.8-87 [90m2024-06-26[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## Formula 1.2-5 [90m2023-02-24[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## fs 1.6.5 [90m2024-10-30[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## genefilter 1.88.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## geneplotter 1.84.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## generics 0.1.3 [90m2022-07-05[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## GenomeInfoDb * 1.42.1 [90m2024-11-28[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## GenomeInfoDbData 1.2.13 [90m2024-12-10[39m [90m[1][39m [90mBioconductor[39m
## GenomicAlignments 1.42.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## GenomicFeatures * 1.58.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## GenomicFiles 1.42.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## GenomicRanges * 1.58.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## GGally 2.2.1 [90m2024-02-14[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## ggbio * 1.54.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## ggplot2 * 3.5.1 [90m2024-04-23[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## ggstats 0.7.0 [90m2024-09-22[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## glue 1.8.0 [90m2024-09-30[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## graph 1.84.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## gridExtra * 2.3 [90m2017-09-09[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## gtable 0.3.6 [90m2024-10-25[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## Hmisc 5.2-1 [90m2024-12-02[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## hms 1.1.3 [90m2023-03-21[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## htmlTable 2.4.3 [90m2024-07-21[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## htmltools 0.5.8.1 [90m2024-04-04[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## htmlwidgets 1.6.4 [90m2023-12-06[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## httpuv 1.6.15 [90m2024-03-26[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## httr 1.4.7 [90m2023-08-15[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## httr2 1.0.7 [90m2024-11-26[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## hwriter 1.3.2.1 [90m2022-04-08[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## IRanges * 2.40.1 [90m2024-12-05[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## iterators * 1.0.14 [90m2022-02-05[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## jquerylib 0.1.4 [90m2021-04-26[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## jsonlite 1.8.9 [90m2024-09-20[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## KEGGREST 1.46.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## knitr * 1.49 [90m2024-11-08[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## knitrBootstrap 1.0.3 [90m2024-02-06[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## labeling 0.4.3 [90m2023-08-29[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## later 1.4.1 [90m2024-11-27[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## lattice 0.22-6 [90m2024-03-20[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## lazyeval 0.2.2 [90m2019-03-15[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## lifecycle 1.0.4 [90m2023-11-07[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## limma * 3.62.1 [90m2024-11-03[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## locfit * 1.5-9.10 [90m2024-06-24[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## lubridate 1.9.4 [90m2024-12-08[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## magrittr 2.0.3 [90m2022-03-30[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## markdown 1.13 [90m2024-06-04[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## Matrix 1.7-1 [90m2024-10-18[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## MatrixGenerics * 1.18.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## matrixStats * 1.4.1 [90m2024-09-08[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## memoise 2.0.1 [90m2021-11-26[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## mgcv * 1.9-1 [90m2023-12-21[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## mime 0.12 [90m2021-09-28[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## miniUI 0.1.1.1 [90m2018-05-18[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## munsell 0.5.1 [90m2024-04-01[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## nlme * 3.1-166 [90m2024-08-14[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## nnet 7.3-19 [90m2023-05-03[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## OrganismDbi 1.48.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## pasilla * 1.34.0 [90m2024-11-05[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## pheatmap * 1.0.12 [90m2019-01-04[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## pillar 1.9.0 [90m2023-03-22[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## pkgbuild 1.4.5 [90m2024-10-28[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## pkgconfig 2.0.3 [90m2019-09-22[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## pkgdown 2.1.1 [90m2024-09-17[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## pkgload 1.4.0 [90m2024-06-28[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## plyr 1.8.9 [90m2023-10-02[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## png 0.1-8 [90m2022-11-29[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## prettyunits 1.2.0 [90m2023-09-24[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## profvis 0.4.0 [90m2024-09-20[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## progress 1.2.3 [90m2023-12-06[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## promises 1.3.2 [90m2024-11-28[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## ProtGenerics 1.38.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## purrr 1.0.2 [90m2023-08-10[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## qvalue 2.38.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## R6 2.5.1 [90m2021-08-19[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## ragg 1.3.3 [90m2024-09-11[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rappdirs 0.3.3 [90m2021-01-31[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## RBGL 1.82.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## RColorBrewer * 1.1-3 [90m2022-04-03[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## Rcpp 1.0.13-1 [90m2024-11-02[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## RCurl 1.98-1.16 [90m2024-07-11[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## RefManageR 1.4.0 [90m2022-09-30[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## regionReport * 1.41.0 [90m2024-12-12[39m [90m[1][39m [90mBioconductor[39m
## remotes 2.5.0 [90m2024-03-17[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## reshape2 1.4.4 [90m2020-04-09[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## restfulr 0.0.15 [90m2022-06-16[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rjson 0.2.23 [90m2024-09-16[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rlang 1.1.4 [90m2024-06-04[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rmarkdown 2.29 [90m2024-11-04[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rngtools 1.5.2 [90m2021-09-20[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rpart 4.1.23 [90m2023-12-05[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## Rsamtools 2.22.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## RSQLite 2.3.9 [90m2024-12-03[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rstudioapi 0.17.1 [90m2024-10-22[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## rtracklayer 1.66.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## S4Arrays 1.6.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## S4Vectors * 0.44.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## sass 0.4.9 [90m2024-03-15[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## scales 1.3.0 [90m2023-11-28[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## sessioninfo * 1.2.2 [90m2021-12-06[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## shiny 1.9.1 [90m2024-08-01[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## SparseArray 1.6.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## statmod 1.5.0 [90m2023-01-06[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## stringi 1.8.4 [90m2024-05-06[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## stringr 1.5.1 [90m2023-11-14[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## SummarizedExperiment * 1.36.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## survival 3.7-0 [90m2024-06-05[39m [90m[3][39m [90mCRAN (R 4.4.2)[39m
## systemfonts 1.1.0 [90m2024-05-15[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## textshaping 0.4.1 [90m2024-12-06[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## tibble 3.2.1 [90m2023-03-20[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## tidyr 1.3.1 [90m2024-01-24[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## tidyselect 1.2.1 [90m2024-03-11[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## timechange 0.3.0 [90m2024-01-18[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## TxDb.Hsapiens.UCSC.hg19.knownGene * 3.2.2 [90m2024-06-24[39m [90m[1][39m [90mBioconductor[39m
## txdbmaker 1.2.1 [90m2024-11-25[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## UCSC.utils 1.2.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## urlchecker 1.0.1 [90m2021-11-30[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## usethis 3.1.0 [90m2024-11-26[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## utf8 1.2.4 [90m2023-10-22[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## VariantAnnotation 1.52.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## vctrs 0.6.5 [90m2023-12-01[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## whisker * 0.4.1 [90m2022-12-05[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## withr 3.0.2 [90m2024-10-28[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## xfun 0.49 [90m2024-10-31[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## XML 3.99-0.17 [90m2024-06-25[39m [90m[1][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## xml2 1.3.6 [90m2023-12-04[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## xtable 1.8-4 [90m2019-04-21[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## XVector 0.46.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
## yaml 2.3.10 [90m2024-07-26[39m [90m[2][39m [1m[35mRSPM (R 4.4.0)[39m[22m
## zlibbioc 1.52.0 [90m2024-10-29[39m [90m[1][39m [90mBioconductor 3.20 (R 4.4.2)[39m
##
## [90m [1] /__w/_temp/Library[39m
## [90m [2] /usr/local/lib/R/site-library[39m
## [90m [3] /usr/local/lib/R/library[39m
##
## [1m[36m──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────[39m[22m
Pandoc version used: 3.5.
This report was created with regionReport
(Collado-Torres, Jaffe, and Leek, 2016) using rmarkdown
(Allaire, Xie, Dervieux et al., 2024) while knitr
(Xie, 2014) and DT
(Xie, Cheng, and Tan, 2024) were running behind the scenes. whisker
(de Jonge, 2022) was used for creating templates for the pvalueVars
and densityVars
.
Citations made with RefManageR (McLean, 2017). The BibTeX file can be found here.
[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] 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.
[3] 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.
[4] 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.
[5] 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.
[6] 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.
[7] H. Wickham. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York, 2016. ISBN: 978-3-319-24277-4. URL: https://ggplot2.tidyverse.org.
[8] 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.
[9] 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.
[10] 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.
[11] E. de Jonge. whisker: mustache for R, Logicless Templating. R package version 0.4.1. 2022. URL: https://github.com/edwindj/whisker.