derfinder quick start guide
30 April 2018
Package
derfinder 1.14.0
1 Basics
1.1 Install derfinder
R is an open-source statistical environment which can be easily modified to enhance its functionality via packages. derfinder 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 derfinder by using the following commands in your R session:
## try http:// if https:// URLs are not supported
source("https://bioconductor.org/biocLite.R")
biocLite("derfinder")
## Check that you have a valid Bioconductor installation
biocValid()1.2 Required knowledge
derfinder 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 Rsamtools, GenomicAlignments and rtracklayer that allow you to import the data. A derfinder user is not expected to deal with those packages directly but will need to be familiar with GenomicRanges to understand the results derfinder generates. It might also prove to be highly beneficial to check the BiocParallel package for performing parallel computations.
If you are asking yourself the question “Where do I start using Bioconductor?” you might be interested in this blog post.
1.3 Asking for help
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 derfinder 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 would like to highlight the derfinder user Jessica Hekman. She has used derfinder with non-human data, and in the process of doing so discovered some small bugs or sections of the documentation that were not clear.
1.4 Citing derfinder
We hope that derfinder will be useful for your research. Please use the following information to cite the package and the overall approach. Thank you!
## Citation info
citation('derfinder')##
## 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 (URL: http://doi.org/10.1093/nar/gkw852), <URL:
## http://nar.oxfordjournals.org/content/early/2016/09/29/nar.gkw852>.
##
## Frazee AC, Sabunciyan S, Hansen KD, Irizarry RA, Leek JT (2014).
## "Differential expression analysis of RNA-seq data at single-base
## resolution." _Biostatistics_, *15 (3)*, 413-426. doi:
## 10.1093/biostatistics/kxt053 (URL:
## http://doi.org/10.1093/biostatistics/kxt053), <URL:
## http://biostatistics.oxfordjournals.org/content/15/3/413.long>.
##
## Collado-Torres L, Jaffe AE, Leek JT (2017). _derfinder:
## Annotation-agnostic differential expression analysis of RNA-seq data
## at base-pair resolution via the DER Finder approach_. doi:
## 10.18129/B9.bioc.derfinder (URL:
## http://doi.org/10.18129/B9.bioc.derfinder),
## https://github.com/lcolladotor/derfinder - R package version 1.14.0,
## <URL: http://www.bioconductor.org/packages/derfinder>.
##
## To see these entries in BibTeX format, use 'print(<citation>,
## bibtex=TRUE)', 'toBibtex(.)', or set
## 'options(citation.bibtex.max=999)'.2 Quick start to using to derfinder
Here is a very quick example of a DER Finder analysis. This analysis is explained in more detail later on in this document.
## Load libraries
library('derfinder')
library('derfinderData')
library('GenomicRanges')
## Determine the files to use and fix the names
files <- rawFiles(system.file('extdata', 'AMY', package = 'derfinderData'),
samplepatt = 'bw', fileterm = NULL)
names(files) <- gsub('.bw', '', names(files))
## Load the data from disk -- On Windows you have to load data from Bam files
fullCov <- fullCoverage(files = files, chrs = 'chr21', verbose = FALSE)
## Get the region matrix of Expressed Regions (ERs)
regionMat <- regionMatrix(fullCov, cutoff = 30, L = 76, verbose = FALSE)
## Get pheno table
pheno <- subset(brainspanPheno, structure_acronym == 'AMY')
## Identify which ERs are differentially expressed, that is, find the DERs
library('DESeq2')
## Round matrix
counts <- round(regionMat$chr21$coverageMatrix)
## Round matrix and specify design
dse <- DESeqDataSetFromMatrix(counts, pheno, ~ group + gender)
## Perform DE analysis
dse <- DESeq(dse, test = 'LRT', reduced = ~ gender, fitType = 'local')
## Extract results
mcols(regionMat$chr21$regions) <- c(mcols(regionMat$chr21$regions), results(dse))
## Save info in an object with a shorter name
ers <- regionMat$chr21$regions
ers3 Introduction
derfinder is an R package that implements the DER Finder approach (Frazee, Sabunciyan, Hansen, Irizarry, et al., 2014) for RNA-seq data. Briefly, this approach is an alternative to feature-counting and transcript assembly. The basic idea is to identify contiguous base-pairs in the genome that present differential expression signal. These base-pairs are grouped into _d_ifferentially _e_xpressed _r_regions (DERs). This approach is annotation-agnostic which is a feature you might be interested in. In particular, derfinder contains functions that allow you to identify DERs via two alternative methods. You can find more details on the full derfinder users guide.
4 Sample DER Finder analysis
This is a brief overview of what a DER Finder analysis looks like. In particular, here we will be identifying expressed regions (ERs) without relying on annotation. Next, we’ll identify candidate differentially expressed regions (DERs). Finally, we’ll compare the DERs with known annotation features.
We first load the required packages.
## Load libraries
library('derfinder')
library('derfinderData')
library('GenomicRanges')Next, we need to locate the chromosome 21 coverage files for a set of 12 samples. These samples are a small subset from the BrainSpan Atlas of the Human Brain (BrainSpan, 2011) publicly available data. The function rawFiles() helps us in locating these files.
## Determine the files to use and fix the names
files <- rawFiles(system.file('extdata', 'AMY', package = 'derfinderData'),
samplepatt = 'bw', fileterm = NULL)
names(files) <- gsub('.bw', '', names(files))Next, we can load the full coverage data into memory using the fullCoverage() function. Note that the BrainSpan data is already normalized by the total number of mapped reads in each sample. However, that won’t be the case with most data sets in which case you might want to use the totalMapped and targetSize arguments. The function getTotalMapped() will be helpful to get this information.
## Load the data from disk
fullCov <- fullCoverage(files = files, chrs = 'chr21', verbose = FALSE,
totalMapped = rep(1, length(files)), targetSize = 1)Now that we have the data, we can identify expressed regions (ERs) by using a cutoff of 30 on the base-level mean coverage from these 12 samples. Once the regions have been identified, we can calculate a coverage matrix with one row per ER and one column per sample (12 in this case). For doing this calculation we need to know the length of the sequence reads, which in this study were 76 bp long.
## Get the region matrix of Expressed Regions (ERs)
regionMat <- regionMatrix(fullCov, cutoff = 30, L = 76, verbose = FALSE)regionMatrix() returns a list of elements, each with three pieces of output. The actual ERs are arranged in a GRanges object named regions.
## regions output
regionMat$chr21$regions## GRanges object with 45 ranges and 6 metadata columns:
## seqnames ranges strand | value area
## <Rle> <IRanges> <Rle> | <numeric> <numeric>
## 1 chr21 9827018-9827582 * | 313.671741938907 177224.534195483
## 2 chr21 15457301-15457438 * | 215.084607424943 29681.6758246422
## 3 chr21 20230140-20230192 * | 38.8325471608144 2058.12499952316
## 4 chr21 20230445-20230505 * | 41.3245082031834 2520.79500039419
## 5 chr21 27253318-27253543 * | 34.9131305372469 7890.3675014178
## .. ... ... ... . ... ...
## 41 chr21 33039644-33039688 * | 34.4705371008979 1551.17416954041
## 42 chr21 33040784-33040798 * | 32.134222178989 482.013332684835
## 43 chr21 33040890 * | 30.0925002098083 30.0925002098083
## 44 chr21 33040900-33040901 * | 30.1208333174388 60.2416666348775
## 45 chr21 48019401-48019414 * | 31.1489284609755 436.084998453657
## indexStart indexEnd cluster clusterL
## <integer> <integer> <Rle> <Rle>
## 1 1 565 1 565
## 2 566 703 2 138
## 3 704 756 3 366
## 4 757 817 3 366
## 5 818 1043 4 765
## .. ... ... ... ...
## 41 2180 2224 17 45
## 42 2225 2239 18 118
## 43 2240 2240 18 118
## 44 2241 2242 18 118
## 45 2243 2256 19 14
## -------
## seqinfo: 1 sequence from an unspecified genomebpCoverage is the base-level coverage list which can then be used for plotting.
## Base-level coverage matrices for each of the regions
## Useful for plotting
lapply(regionMat$chr21$bpCoverage[1:2], head, n = 2)## $`1`
## HSB97 HSB92 HSB178 HSB159 HSB153 HSB113 HSB130 HSB136 HSB126 HSB145 HSB123
## 1 36.54 47.40 15.71 3.52 16.71 93.20 5.62 98.34 28.22 5.88 3.32
## 2 51.42 65.04 29.35 5.83 30.53 124.76 11.32 199.28 63.68 10.39 7.25
## HSB135
## 1 185.17
## 2 374.85
##
## $`2`
## HSB97 HSB92 HSB178 HSB159 HSB153 HSB113 HSB130 HSB136 HSB126 HSB145
## 566 22.21 13.95 23.2 4.97 38.61 45.59 34.75 104.21 15.92 19.87
## 567 22.21 13.95 23.2 4.97 38.65 45.59 35.15 104.30 15.92 19.87
## HSB123 HSB135
## 566 7.94 141.61
## 567 7.94 141.64## Check dimensions. First region is 565 long, second one is 138 bp long.
## The columns match the number of samples (12 in this case).
lapply(regionMat$chr21$bpCoverage[1:2], dim)## $`1`
## [1] 565 12
##
## $`2`
## [1] 138 12The end result of the coverage matrix is shown below. Note that the coverage has been adjusted for read length. Because reads might not fully align inside a given region, the numbers are generally not integers but can be rounded if needed.
## Dimensions of the coverage matrix
dim(regionMat$chr21$coverageMatrix)## [1] 45 12## Coverage for each region. This matrix can then be used with limma or other pkgs
head(regionMat$chr21$coverageMatrix)## HSB97 HSB92 HSB178 HSB159 HSB153 HSB113
## 1 1325.9622371 1168.519079 1407.288552 297.939342 1461.2986829 3653.1093346
## 2 200.7834228 127.543553 193.695657 67.940131 313.8513139 333.3740816
## 3 20.4623685 23.758816 33.127368 28.145395 52.8756585 35.3828948
## 4 28.2011836 29.546183 33.576974 31.244343 51.6146040 42.3398681
## 5 35.5428946 100.085790 97.151316 78.069605 36.3046051 77.7402631
## 6 0.3956579 1.136447 1.139079 1.000132 0.4221053 0.7988158
## HSB130 HSB136 HSB126 HSB145 HSB123 HSB135
## 1 1106.722895 5570.221054 1397.068687 1330.158818 277.072106 8987.460401
## 2 267.354342 1162.313418 463.909476 257.114210 99.987237 1198.713552
## 3 23.483947 17.168947 30.725394 22.895921 20.153553 16.786842
## 4 24.724736 19.309606 41.094474 33.802632 29.931579 32.634080
## 5 171.861974 93.503158 115.059342 90.950526 168.939342 180.638684
## 6 2.231053 1.007895 1.473421 1.171316 1.770263 1.697368We can then use the coverage matrix and packages such as limma, DESeq2 or edgeR to identify which ERs are differentially expressed. Here we’ll use DESeq2 for which we need some phenotype data.
## Get pheno table
pheno <- subset(brainspanPheno, structure_acronym == 'AMY')Now we can identify the DERs using a rounded version of the coverage matrix.
## Identify which ERs are differentially expressed, that is, find the DERs
library('DESeq2')## Loading required package: SummarizedExperiment## Loading required package: Biobase## Welcome to Bioconductor
##
## Vignettes contain introductory material; view with
## 'browseVignettes()'. To cite Bioconductor, see
## 'citation("Biobase")', and for packages 'citation("pkgname")'.## Loading required package: DelayedArray## Loading required package: matrixStats##
## Attaching package: 'matrixStats'## The following objects are masked from 'package:Biobase':
##
## anyMissing, rowMedians## Loading required package: BiocParallel##
## Attaching package: 'DelayedArray'## The following objects are masked from 'package:matrixStats':
##
## colMaxs, colMins, colRanges, rowMaxs, rowMins, rowRanges## The following objects are masked from 'package:base':
##
## aperm, apply## Round matrix
counts <- round(regionMat$chr21$coverageMatrix)
## Round matrix and specify design
dse <- DESeqDataSetFromMatrix(counts, pheno, ~ group + gender)## converting counts to integer mode## Perform DE analysis
dse <- DESeq(dse, test = 'LRT', reduced = ~ gender, fitType = 'local')## estimating size factors## estimating dispersions## gene-wise dispersion estimates## mean-dispersion relationship## final dispersion estimates## fitting model and testing## Extract results
mcols(regionMat$chr21$regions) <- c(mcols(regionMat$chr21$regions),
results(dse))
## Save info in an object with a shorter name
ers <- regionMat$chr21$regions
ers## GRanges object with 45 ranges and 12 metadata columns:
## seqnames ranges strand | value area
## <Rle> <IRanges> <Rle> | <numeric> <numeric>
## 1 chr21 9827018-9827582 * | 313.671741938907 177224.534195483
## 2 chr21 15457301-15457438 * | 215.084607424943 29681.6758246422
## 3 chr21 20230140-20230192 * | 38.8325471608144 2058.12499952316
## 4 chr21 20230445-20230505 * | 41.3245082031834 2520.79500039419
## 5 chr21 27253318-27253543 * | 34.9131305372469 7890.3675014178
## .. ... ... ... . ... ...
## 41 chr21 33039644-33039688 * | 34.4705371008979 1551.17416954041
## 42 chr21 33040784-33040798 * | 32.134222178989 482.013332684835
## 43 chr21 33040890 * | 30.0925002098083 30.0925002098083
## 44 chr21 33040900-33040901 * | 30.1208333174388 60.2416666348775
## 45 chr21 48019401-48019414 * | 31.1489284609755 436.084998453657
## indexStart indexEnd cluster clusterL baseMean
## <integer> <integer> <Rle> <Rle> <numeric>
## 1 1 565 1 565 2846.28716155756
## 2 566 703 2 138 451.519643002767
## 3 704 756 3 366 29.5780895864663
## 4 757 817 3 366 36.0602688380417
## 5 818 1043 4 765 101.646763062388
## .. ... ... ... ... ...
## 41 2180 2224 17 45 20.7820346002786
## 42 2225 2239 18 118 6.41054227728292
## 43 2240 2240 18 118 0.129716561731455
## 44 2241 2242 18 118 0.702291488002413
## 45 2243 2256 19 14 5.29329263119445
## log2FoldChange lfcSE stat
## <numeric> <numeric> <numeric>
## 1 -1.30120370440296 0.900113799302452 1.18751549559329
## 2 -1.063513200629 0.787523066873229 2.66980711856959
## 3 -0.536464039142821 0.307664617466879 16.3262317270667
## 4 -0.625613734632087 0.294301541893243 8.96934291317559
## 5 -0.0887821239498447 0.344720516444579 2.15123765308093
## .. ... ... ...
## 41 -0.654516821124693 0.32647613471869 1.20103944253287
## 42 -0.570136770711043 0.475008341718673 0.454717599746111
## 43 -0.745913545383179 3.11653965011299 -0.0333220146504161
## 44 -0.646216179154895 3.04557743767828 0.260500694468277
## 45 -0.320976518752516 0.744646534294096 29.1670784564985
## pvalue padj
## <numeric> <numeric>
## 1 0.27583128369531 1
## 2 0.102268416966693 0.920415752700237
## 3 5.33207203205748e-05 0.00119971620721293
## 4 0.00274547317006963 0.0411820975510445
## 5 0.14245501436384 0.947562263665391
## .. ... ...
## 41 0.273114026168067 1
## 42 0.500103113665161 1
## 43 1 1
## 44 0.609776379540419 1
## 45 6.63978743599099e-08 2.98790434619595e-06
## -------
## seqinfo: 1 sequence from an unspecified genomeWe can then compare the DERs against known annotation to see which DERs overlap known exons, introns, or intergenic regions. A way to visualize this information is via a Venn diagram which we can create using vennRegions() from the derfinderPlot package as shown in Figure 1.
## Find overlaps between regions and summarized genomic annotation
annoRegs <- annotateRegions(ers, genomicState$fullGenome, verbose = FALSE)
library('derfinderPlot')
venn <- vennRegions(annoRegs, counts.col = 'blue',
main = 'Venn diagram using TxDb.Hsapiens.UCSC.hg19.knownGene annotation')
Figure 1: Venn diagram showing ERs by annotation class
We can also identify the nearest annotated feature. In this case, we’ll look for the nearest known gene from the UCSC hg19 annotation.
## Load database of interest
library('TxDb.Hsapiens.UCSC.hg19.knownGene')## Loading required package: GenomicFeatures## Loading required package: AnnotationDbitxdb <- keepSeqlevels(TxDb.Hsapiens.UCSC.hg19.knownGene, 'chr21')
## Find nearest feature
library('bumphunter')## Loading required package: foreach## Loading required package: iterators## Loading required package: locfit## locfit 1.5-9.1 2013-03-22genes <- annotateTranscripts(txdb)## No annotationPackage supplied. Trying org.Hs.eg.db.## Loading required package: org.Hs.eg.db## ## Getting TSS and TSE.## Getting CSS and CSE.## Getting exons.## Annotating genes.annotation <- matchGenes(ers, subject = genes)
## View annotation results
head(annotation)## name
## 1 <NA>
## 2 LIPI
## 3 TMPRSS15
## 4 TMPRSS15
## 5 APP
## 6 APP
## annotation
## 1 <NA>
## 2 NM_001302998 NM_001302999 NM_001303000 NM_001303001 NM_145317 NM_198996 NP_001289927 NP_001289928 NP_001289929 NP_001289930 NP_945347 XM_006723965 XP_006724028
## 3 NM_002772 NP_002763 XM_011529654 XM_011529655 XM_011529656 XM_011529657 XM_011529658 XM_011529659 XP_011527956 XP_011527957 XP_011527958 XP_011527959 XP_011527960 XP_011527961
## 4 NM_002772 NP_002763 XM_011529654 XM_011529655 XM_011529656 XM_011529657 XM_011529658 XM_011529659 XP_011527956 XP_011527957 XP_011527958 XP_011527959 XP_011527960 XP_011527961
## 5 NM_000484 NM_001136016 NM_001136129 NM_001136130 NM_001136131 NM_001204301 NM_001204302 NM_001204303 NM_201413 NM_201414 NP_000475 NP_001129488 NP_001129601 NP_001129602 NP_001129603 NP_001191230 NP_001191231 NP_001191232 NP_958816 NP_958817 XM_024452075 XP_024307843
## 6 NM_000484 NM_001136016 NM_001136129 NM_001136130 NM_001136131 NM_001204301 NM_001204302 NM_001204303 NM_201413 NM_201414 NP_000475 NP_001129488 NP_001129601 NP_001129602 NP_001129603 NP_001191230 NP_001191231 NP_001191232 NP_958816 NP_958817 XM_024452075 XP_024307843
## description region distance subregion insideDistance exonnumber
## 1 close to 3' close to 3' 815 <NA> NA NA
## 2 downstream downstream 125774 <NA> NA NA
## 3 upstream upstream 454170 <NA> NA NA
## 4 upstream upstream 454475 <NA> NA NA
## 5 inside exon inside 289903 inside exon 0 16
## 6 inside exon inside 289899 inside exon 0 16
## nexons UTR strand geneL codingL Geneid subjectHits
## 1 1 <NA> + 60 NA 100500815 15
## 2 8 <NA> - 102077 101886 149998 149
## 3 25 <NA> - 134537 133653 5651 579
## 4 25 <NA> - 134537 133653 5651 579
## 5 16 3'UTR - 290585 230434 351 354
## 6 16 3'UTR - 290585 230434 351 354## You can use derfinderPlot::plotOverview() to visualize this informationWe can check the base-level coverage information for some of our DERs. In this example we do so for the first 5 ERs (Figures ??, ??, ??, ??, ??).
## Extract the region coverage
regionCov <- regionMat$chr21$bpCoverage
plotRegionCoverage(regions = ers, regionCoverage = regionCov,
groupInfo = pheno$group, nearestAnnotation = annotation,
annotatedRegions = annoRegs, whichRegions = seq_len(5), txdb = txdb,
scalefac = 1, ask = FALSE, verbose = FALSE)
Figure 2: Base-pair resolution plot of differentially expressed region 1
Figure 2: Base-pair resolution plot of differentially expressed region 2
Figure 2: Base-pair resolution plot of differentially expressed region 3
Figure 2: Base-pair resolution plot of differentially expressed region 4
Figure 2: Base-pair resolution plot of differentially expressed region 5
You can then use the regionReport package to generate interactive HTML reports exploring the results.
If you are interested in using derfinder we recommend checking the derfinder users guide and good luck with your analyses!
5 Reproducibility
This package was made possible thanks to:
- R (R Core Team, 2018)
- AnnotationDbi (Pagès, Carlson, Falcon, and Li, 2017)
- BiocParallel (Morgan, Obenchain, Lang, Thompson, et al., 2018)
- bumphunter (Jaffe, Murakami, Lee, Leek, et al., 2012) and (Jaffe, Murakami, Lee, Leek, et al., 2012)
- derfinderHelper
- GenomeInfoDb (Arora, Morgan, Carlson, and Pagès, 2017)
- GenomicAlignments (Lawrence, Huber, Pagès, Aboyoun, et al., 2013)
- GenomicFeatures (Lawrence, Huber, Pagès, Aboyoun, et al., 2013)
- GenomicFiles (Bioconductor Package Maintainer, Obenchain, Love, Shepherd, et al., 2018)
- GenomicRanges (Lawrence, Huber, Pagès, Aboyoun, et al., 2013)
- Hmisc (Harrell Jr, Dupont, and others., 2018)
- IRanges (Lawrence, Huber, Pagès, Aboyoun, et al., 2013)
- qvalue (Bass, Dabney, and Robinson, 2015)
- Rsamtools (Morgan, Pagès, Obenchain, and Hayden, 2018)
- rtracklayer (Lawrence, Gentleman, and Carey, 2009)
- S4Vectors (Pagès, Lawrence, and Aboyoun, 2017)
- biovizBase (Yin, Lawrence, and Cook, 2018)
- derfinderData (Collado-Torres, Jaffe, and Leek, 2015)
- devtools (Wickham, Hester, and Chang, 2018)
- ggplot2 (Wickham, 2009)
- knitcitations (Boettiger, 2017)
- knitr (Xie, 2014)
- BiocStyle (Oles, Morgan, and Huber, 2018)
- rmarkdown (Allaire, Xie, McPherson, Luraschi, et al., 2018)
- testthat (Wickham, 2011)
- TxDb.Hsapiens.UCSC.hg19.knownGene (Carlson and Maintainer, 2015)
Code for creating the vignette
## Create the vignette
library('rmarkdown')
system.time(render('derfinder-quickstart.Rmd', 'BiocStyle::html_document'))
## Extract the R code
library('knitr')
knit('derfinder-quickstart.Rmd', tangle = TRUE)## Clean up
file.remove('quickstartRef.bib')## [1] TRUEDate the vignette was generated.
## [1] "2018-04-30 22:42:47 EDT"Wallclock time spent generating the vignette.
## Time difference of 1.262 minsR session information.
## Session info ----------------------------------------------------------------------------------------------------------## setting value
## version R version 3.5.0 (2018-04-23)
## system x86_64, mingw32
## ui RTerm
## language (EN)
## collate C
## tz America/New_York
## date 2018-04-30## Packages --------------------------------------------------------------------------------------------------------------## package * version date source
## acepack 1.4.1 2016-10-29 CRAN (R 3.5.0)
## annotate 1.58.0 2018-04-30 Bioconductor
## AnnotationDbi * 1.42.0 2018-04-30 Bioconductor
## AnnotationFilter 1.4.0 2018-04-30 Bioconductor
## assertthat 0.2.0 2017-04-11 CRAN (R 3.5.0)
## backports 1.1.2 2017-12-13 CRAN (R 3.5.0)
## base * 3.5.0 2018-04-23 local
## base64enc 0.1-3 2015-07-28 CRAN (R 3.5.0)
## bibtex 0.4.2 2017-06-30 CRAN (R 3.5.0)
## Biobase * 2.40.0 2018-04-30 Bioconductor
## BiocGenerics * 0.26.0 2018-04-30 Bioconductor
## BiocInstaller 1.30.0 2018-04-30 Bioconductor
## BiocParallel * 1.14.0 2018-04-30 Bioconductor
## BiocStyle * 2.8.0 2018-04-30 Bioconductor
## biomaRt 2.36.0 2018-04-30 Bioconductor
## Biostrings 2.48.0 2018-04-30 Bioconductor
## biovizBase 1.28.0 2018-05-01 Bioconductor
## bit 1.1-12 2014-04-09 CRAN (R 3.5.0)
## bit64 0.9-7 2017-05-08 CRAN (R 3.5.0)
## bitops 1.0-6 2013-08-17 CRAN (R 3.5.0)
## blob 1.1.1 2018-03-25 CRAN (R 3.5.0)
## bookdown 0.7 2018-02-18 CRAN (R 3.5.0)
## BSgenome 1.48.0 2018-05-01 Bioconductor
## bumphunter * 1.22.0 2018-05-01 Bioconductor
## checkmate 1.8.5 2017-10-24 CRAN (R 3.5.0)
## cluster 2.0.7-1 2018-04-13 CRAN (R 3.5.0)
## codetools 0.2-15 2016-10-05 CRAN (R 3.5.0)
## colorspace 1.3-2 2016-12-14 CRAN (R 3.5.0)
## compiler 3.5.0 2018-04-23 local
## curl 3.2 2018-03-28 CRAN (R 3.5.0)
## data.table 1.10.4-3 2017-10-27 CRAN (R 3.5.0)
## datasets * 3.5.0 2018-04-23 local
## DBI 0.8 2018-03-02 CRAN (R 3.5.0)
## DelayedArray * 0.6.0 2018-04-30 Bioconductor
## derfinder * 1.14.0 2018-05-01 Bioconductor
## derfinderData * 0.113.0 2018-04-27 Bioconductor
## derfinderHelper 1.14.0 2018-04-30 Bioconductor
## derfinderPlot * 1.14.0 2018-05-01 Bioconductor
## DESeq2 * 1.20.0 2018-05-01 Bioconductor
## devtools * 1.13.5 2018-02-18 CRAN (R 3.5.0)
## dichromat 2.0-0 2013-01-24 CRAN (R 3.5.0)
## digest 0.6.15 2018-01-28 CRAN (R 3.5.0)
## doRNG 1.6.6 2017-04-10 CRAN (R 3.5.0)
## ensembldb 2.4.0 2018-05-01 Bioconductor
## evaluate 0.10.1 2017-06-24 CRAN (R 3.5.0)
## foreach * 1.4.4 2017-12-12 CRAN (R 3.5.0)
## foreign 0.8-70 2017-11-28 CRAN (R 3.5.0)
## Formula 1.2-2 2017-07-10 CRAN (R 3.5.0)
## genefilter 1.62.0 2018-04-30 Bioconductor
## geneplotter 1.58.0 2018-04-30 Bioconductor
## GenomeInfoDb * 1.16.0 2018-04-30 Bioconductor
## GenomeInfoDbData 1.1.0 2018-04-27 Bioconductor
## GenomicAlignments 1.16.0 2018-04-30 Bioconductor
## GenomicFeatures * 1.32.0 2018-05-01 Bioconductor
## GenomicFiles 1.16.0 2018-05-01 Bioconductor
## GenomicRanges * 1.32.0 2018-04-30 Bioconductor
## GGally 1.3.2 2017-08-02 CRAN (R 3.5.0)
## ggbio 1.28.0 2018-05-01 Bioconductor
## ggplot2 2.2.1 2016-12-30 CRAN (R 3.5.0)
## graph 1.58.0 2018-04-30 Bioconductor
## graphics * 3.5.0 2018-04-23 local
## grDevices * 3.5.0 2018-04-23 local
## grid 3.5.0 2018-04-23 local
## gridExtra 2.3 2017-09-09 CRAN (R 3.5.0)
## gtable 0.2.0 2016-02-26 CRAN (R 3.5.0)
## highr 0.6 2016-05-09 CRAN (R 3.5.0)
## Hmisc 4.1-1 2018-01-03 CRAN (R 3.5.0)
## htmlTable 1.11.2 2018-01-20 CRAN (R 3.5.0)
## htmltools 0.3.6 2017-04-28 CRAN (R 3.5.0)
## htmlwidgets 1.2 2018-04-19 CRAN (R 3.5.0)
## httr 1.3.1 2017-08-20 CRAN (R 3.5.0)
## IRanges * 2.14.0 2018-04-30 Bioconductor
## iterators * 1.0.9 2017-12-12 CRAN (R 3.5.0)
## jsonlite 1.5 2017-06-01 CRAN (R 3.5.0)
## knitcitations * 1.0.8 2017-07-04 CRAN (R 3.5.0)
## knitr 1.20 2018-02-20 CRAN (R 3.5.0)
## lattice 0.20-35 2017-03-25 CRAN (R 3.5.0)
## latticeExtra 0.6-28 2016-02-09 CRAN (R 3.5.0)
## lazyeval 0.2.1 2017-10-29 CRAN (R 3.5.0)
## limma 3.36.0 2018-04-30 Bioconductor
## locfit * 1.5-9.1 2013-04-20 CRAN (R 3.5.0)
## lubridate 1.7.4 2018-04-11 CRAN (R 3.5.0)
## magrittr 1.5 2014-11-22 CRAN (R 3.5.0)
## Matrix 1.2-14 2018-04-13 CRAN (R 3.5.0)
## matrixStats * 0.53.1 2018-02-11 CRAN (R 3.5.0)
## memoise 1.1.0 2017-04-21 CRAN (R 3.5.0)
## methods * 3.5.0 2018-04-23 local
## munsell 0.4.3 2016-02-13 CRAN (R 3.5.0)
## nnet 7.3-12 2016-02-02 CRAN (R 3.5.0)
## org.Hs.eg.db * 3.6.0 2018-04-27 Bioconductor
## OrganismDbi 1.22.0 2018-05-01 Bioconductor
## parallel * 3.5.0 2018-04-23 local
## pillar 1.2.2 2018-04-26 CRAN (R 3.5.0)
## pkgconfig 2.0.1 2017-03-21 CRAN (R 3.5.0)
## pkgmaker 0.22 2014-05-14 CRAN (R 3.5.0)
## plyr 1.8.4 2016-06-08 CRAN (R 3.5.0)
## prettyunits 1.0.2 2015-07-13 CRAN (R 3.5.0)
## progress 1.1.2 2016-12-14 CRAN (R 3.5.0)
## ProtGenerics 1.12.0 2018-04-30 Bioconductor
## qvalue 2.12.0 2018-04-30 Bioconductor
## R6 2.2.2 2017-06-17 CRAN (R 3.5.0)
## RBGL 1.56.0 2018-04-30 Bioconductor
## RColorBrewer 1.1-2 2014-12-07 CRAN (R 3.5.0)
## Rcpp 0.12.16 2018-03-13 CRAN (R 3.5.0)
## RCurl 1.95-4.10 2018-01-04 CRAN (R 3.5.0)
## RefManageR 1.2.0 2018-04-25 CRAN (R 3.5.0)
## registry 0.5 2017-12-03 CRAN (R 3.5.0)
## reshape 0.8.7 2017-08-06 CRAN (R 3.5.0)
## reshape2 1.4.3 2017-12-11 CRAN (R 3.5.0)
## rlang 0.2.0 2018-02-20 CRAN (R 3.5.0)
## rmarkdown 1.9 2018-03-01 CRAN (R 3.5.0)
## rngtools 1.2.4 2014-03-06 CRAN (R 3.5.0)
## rpart 4.1-13 2018-02-23 CRAN (R 3.5.0)
## rprojroot 1.3-2 2018-01-03 CRAN (R 3.5.0)
## Rsamtools 1.32.0 2018-04-30 Bioconductor
## RSQLite 2.1.0 2018-03-29 CRAN (R 3.5.0)
## rstudioapi 0.7 2017-09-07 CRAN (R 3.5.0)
## rtracklayer 1.40.0 2018-05-01 Bioconductor
## S4Vectors * 0.18.0 2018-04-30 Bioconductor
## scales 0.5.0 2017-08-24 CRAN (R 3.5.0)
## splines 3.5.0 2018-04-23 local
## stats * 3.5.0 2018-04-23 local
## stats4 * 3.5.0 2018-04-23 local
## stringi 1.1.7 2018-03-12 CRAN (R 3.5.0)
## stringr 1.3.0 2018-02-19 CRAN (R 3.5.0)
## SummarizedExperiment * 1.10.0 2018-04-30 Bioconductor
## survival 2.42-3 2018-04-16 CRAN (R 3.5.0)
## tibble 1.4.2 2018-01-22 CRAN (R 3.5.0)
## tools 3.5.0 2018-04-23 local
## TxDb.Hsapiens.UCSC.hg19.knownGene * 3.2.2 2018-04-28 Bioconductor
## utils * 3.5.0 2018-04-23 local
## VariantAnnotation 1.26.0 2018-05-01 Bioconductor
## withr 2.1.2 2018-03-15 CRAN (R 3.5.0)
## xfun 0.1 2018-01-22 CRAN (R 3.5.0)
## XML 3.98-1.11 2018-04-16 CRAN (R 3.5.0)
## xml2 1.2.0 2018-01-24 CRAN (R 3.5.0)
## xtable 1.8-2 2016-02-05 CRAN (R 3.5.0)
## XVector 0.20.0 2018-04-30 Bioconductor
## yaml 2.1.18 2018-03-08 CRAN (R 3.5.0)
## zlibbioc 1.26.0 2018-04-30 Bioconductor6 Bibliography
This vignette was generated using BiocStyle (Oles, Morgan, and Huber, 2018) with knitr (Xie, 2014) and rmarkdown (Allaire, Xie, McPherson, Luraschi, et al., 2018) running behind the scenes.
Citations made with knitcitations (Boettiger, 2017).
[1] J. Allaire, Y. Xie, J. McPherson, J. Luraschi, et al. rmarkdown: Dynamic Documents for R. R package version 1.9. 2018. URL: https://CRAN.R-project.org/package=rmarkdown.
[1] S. Arora, M. Morgan, M. Carlson and H. Pagès. GenomeInfoDb: Utilities for manipulating chromosome and other ‘seqname’ identifiers. 2017. DOI: 10.18129/B9.bioc.GenomeInfoDb.
[1] J. D. S. w. c. f. A. J. Bass, A. Dabney and D. Robinson. qvalue: Q-value estimation for false discovery rate control. R package version 2.12.0. 2015. URL: http://github.com/jdstorey/qvalue.
[1] Bioconductor Package Maintainer, V. Obenchain, M. Love, L. Shepherd, et al. GenomicFiles: Distributed computing by file or by range. R package version 1.16.0. 2018.
[1] C. Boettiger. knitcitations: Citations for ‘Knitr’ Markdown Files. R package version 1.0.8. 2017. URL: https://CRAN.R-project.org/package=knitcitations.
[1] BrainSpan. “Atlas of the Developing Human Brain [Internet]. Funded by ARRA Awards 1RC2MH089921-01, 1RC2MH090047-01, and 1RC2MH089929-01.” 2011. URL: http://www.brainspan.org/.
[1] M. Carlson and B. P. Maintainer. TxDb.Hsapiens.UCSC.hg19.knownGene: Annotation package for TxDb object(s). R package version 3.2.2. 2015.
[1] L. Collado-Torres, A. Jaffe and J. Leek. derfinderData: Processed BigWigs from BrainSpan for examples. R package version 0.113.0. 2015. URL: https://github.com/lcolladotor/derfinderData.
[1] F. E. Harrell Jr, w. c. f. C. Dupont and m. others. Hmisc: Harrell Miscellaneous. R package version 4.1-1. 2018. URL: https://CRAN.R-project.org/package=Hmisc.
[1] A. E. Jaffe, P. Murakami, H. Lee, J. T. Leek, et al. “Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies”. In: International Journal of Epidemiology (2012).
[1] M. Morgan, V. Obenchain, M. Lang, R. Thompson, et al. BiocParallel: Bioconductor facilities for parallel evaluation. R package version 1.14.0. 2018. URL: https://github.com/Bioconductor/BiocParallel.
[1] M. Morgan, H. Pagès, V. Obenchain and N. Hayden. Rsamtools: Binary alignment (BAM), FASTA, variant call (BCF), and tabix file import. R package version 1.32.0. 2018. URL: http://bioconductor.org/packages/release/bioc/html/Rsamtools.html.
[1] A. Oles, M. Morgan and W. Huber. BiocStyle: Standard styles for vignettes and other Bioconductor documents. R package version 2.8.0. 2018. URL: https://github.com/Bioconductor/BiocStyle.
[1] H. Pagès, M. Carlson, S. Falcon and N. Li. AnnotationDbi: Annotation Database Interface. 2017. DOI: 10.18129/B9.bioc.AnnotationDbi.
[1] H. Pagès, M. Lawrence and P. Aboyoun. S4Vectors: S4 implementation of vector-like and list-like objects. 2017. DOI: 10.18129/B9.bioc.S4Vectors.
[1] R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria, 2018. URL: https://www.R-project.org/.
[1] H. Wickham. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York, 2009. ISBN: 978-0-387-98140-6. URL: http://ggplot2.org.
[1] H. Wickham. “testthat: Get Started with Testing”. In: The R Journal 3 (2011), pp. 5–10. URL: https://journal.r-project.org/archive/2011-1/RJournal_2011-1_Wickham.pdf.
[1] H. Wickham, J. Hester and W. Chang. devtools: Tools to Make Developing R Packages Easier. R package version 1.13.5. 2018. URL: https://CRAN.R-project.org/package=devtools.
[1] 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. URL: http://www.crcpress.com/product/isbn/9781466561595.
[1] T. Yin, M. Lawrence and D. Cook. biovizBase: Basic graphic utilities for visualization of genomic data. R package version 1.28.0. 2018.