# Example report using bumphunter results

Leonardo Collado-Torres1,2*

1Lieber Institute for Brain Development, Johns Hopkins Medical Campus
2Center for Computational Biology, Johns Hopkins University

*leo.collado@libd.org

#### 20 January 2018

#### Package

regionReport 1.12.2

# 1 bumphunter example

The bumphunter package can be used for methylation analyses where you are interested in identifying differentially methylated regions. The vignette explains in greater detail the data set we are using in this example.

    ## Load bumphunter
    library('bumphunter')

    ## Loading required package: S4Vectors

    ## Loading required package: stats4

    ## Loading required package: BiocGenerics

    ## Loading required package: parallel

    ##
    ## Attaching package: 'BiocGenerics'

    ## The following objects are masked from 'package:parallel':
    ##
    ##     clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    ##     clusterExport, clusterMap, parApply, parCapply, parLapply,
    ##     parLapplyLB, parRapply, parSapply, parSapplyLB

    ## The following objects are masked from 'package:stats':
    ##
    ##     IQR, mad, sd, var, xtabs

    ## The following objects are masked from 'package:base':
    ##
    ##     Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    ##     as.data.frame, cbind, colMeans, colSums, colnames, do.call,
    ##     duplicated, eval, evalq, get, grep, grepl, intersect,
    ##     is.unsorted, lapply, lengths, mapply, match, mget, order, paste,
    ##     pmax, pmax.int, pmin, pmin.int, rank, rbind, rowMeans, rowSums,
    ##     rownames, sapply, setdiff, sort, table, tapply, union, unique,
    ##     unsplit, which, which.max, which.min

    ##
    ## Attaching package: 'S4Vectors'

    ## The following object is masked from 'package:base':
    ##
    ##     expand.grid

    ## Loading required package: IRanges

    ## Loading required package: GenomeInfoDb

    ## Loading required package: GenomicRanges

    ## Loading required package: foreach

    ## Loading required package: iterators

    ## Loading required package: locfit

    ## locfit 1.5-9.1    2013-03-22

    ## Create data from the vignette
    pos <- list(pos1=seq(1, 1000, 35),
                pos2=seq(2001, 3000, 35),
                pos3=seq(1, 1000, 50))
    chr <- rep(paste0('chr', c(1, 1, 2)), times = sapply(pos, length))
    pos <- unlist(pos, use.names = FALSE)

    ## Find clusters
    cl <- clusterMaker(chr, pos, maxGap = 300)

    ## Build simulated bumps
    Indexes <- split(seq_along(cl), cl)
    beta1 <- rep(0, length(pos))
    for(i in seq(along=Indexes)){
        ind <- Indexes[[i]]
        x <- pos[ind]
        z <- scale(x, median(x), max(x)/12)
        beta1[ind] <- i*(-1)^(i+1)*pmax(1-abs(z)^3,0)^3 ##multiply by i to vary size
    }

    ## Build data
    beta0 <- 3 * sin(2 * pi * pos / 720)
    X <- cbind(rep(1, 20), rep(c(0, 1), each = 10))
    set.seed(23852577)
    error <- matrix(rnorm(20 * length(beta1), 0, 1), ncol = 20)
    y <- t(X[, 1]) %x% beta0 + t(X[, 2]) %x% beta1 + error

    ## Perform bumphunting
    tab <- bumphunter(y, X, chr, pos, cl, cutoff=.5)

    ## [bumphunterEngine] Using a single core (backend: doSEQ, version: 1.4.4).

    ## [bumphunterEngine] Computing coefficients.

    ## [bumphunterEngine] Finding regions.

    ## [bumphunterEngine] Found 15 bumps.

    ## Explore data
    lapply(tab, head)

    ## $table
    ##     chr start  end      value       area cluster indexStart indexEnd  L
    ## 10 chr1  2316 2631 -1.5814747 15.8147473       2         39       48 10
    ## 7  chr2   451  551  1.5891293  4.7673878       3         68       70  3
    ## 2  chr1   456  526  1.0678828  3.2036485       1         14       16  3
    ## 5  chr1  2176 2211  0.7841794  1.5683589       2         35       36  2
    ## 6  chr1  2841 2841  1.2010184  1.2010184       2         54       54  1
    ## 4  chr1   771  771  0.7780902  0.7780902       1         23       23  1
    ##    clusterL
    ## 10       29
    ## 7        20
    ## 2        29
    ## 5        29
    ## 6        29
    ## 4        29
    ##
    ## $coef
    ##             [,1]
    ## [1,]  0.60960932
    ## [2,] -0.09052769
    ## [3,] -0.21482638
    ## [4,]  0.13053755
    ## [5,] -0.21723642
    ## [6,]  0.39934961
    ##
    ## $fitted
    ##             [,1]
    ## [1,]  0.60960932
    ## [2,] -0.09052769
    ## [3,] -0.21482638
    ## [4,]  0.13053755
    ## [5,] -0.21723642
    ## [6,]  0.39934961
    ##
    ## $pvaluesMarginal
    ## [1] NA

Once we have the regions we can proceed to build the required `GRanges` object.

    library('GenomicRanges')

    ## Build GRanges with sequence lengths
    regions <- GRanges(seqnames = tab$table$chr,
        IRanges(start = tab$table$start, end = tab$table$end),
        strand = '*', value = tab$table$value, area = tab$table$area,
        cluster = tab$table$cluster, L = tab$table$L, clusterL = tab$table$clusterL)

    ## Assign chr lengths
    data(hg19Ideogram, package = 'biovizBase')
    seqlengths(regions) <- seqlengths(hg19Ideogram)[names(seqlengths(regions))]

    ## Explore the regions
    regions

    ## GRanges object with 15 ranges and 5 metadata columns:
    ##        seqnames       ranges strand |              value              area
    ##           <Rle>    <IRanges>  <Rle> |          <numeric>         <numeric>
    ##    [1]     chr1 [2316, 2631]      * |  -1.58147472668612  15.8147472668612
    ##    [2]     chr2 [ 451,  551]      * |   1.58912926544627  4.76738779633882
    ##    [3]     chr1 [ 456,  526]      * |   1.06788283103948  3.20364849311844
    ##    [4]     chr1 [2176, 2211]      * |  0.784179442360784  1.56835888472157
    ##    [5]     chr1 [2841, 2841]      * |   1.20101842543356  1.20101842543356
    ##    ...      ...          ...    ... .                ...               ...
    ##   [11]     chr1 [ 631,  631]      * |  0.618602646982714 0.618602646982714
    ##   [12]     chr1 [   1,    1]      * |  0.609609318450113 0.609609318450113
    ##   [13]     chr1 [2911, 2911]      * | -0.576423040796044 0.576423040796044
    ##   [14]     chr2 [ 251,  251]      * | -0.556159535339166 0.556159535339166
    ##   [15]     chr1 [2806, 2806]      * | -0.521605772372843 0.521605772372843
    ##          cluster         L  clusterL
    ##        <numeric> <numeric> <integer>
    ##    [1]         2        10        29
    ##    [2]         3         3        20
    ##    [3]         1         3        29
    ##    [4]         2         2        29
    ##    [5]         2         1        29
    ##    ...       ...       ...       ...
    ##   [11]         1         1        29
    ##   [12]         1         1        29
    ##   [13]         2         1        29
    ##   [14]         3         1        20
    ##   [15]         2         1        29
    ##   -------
    ##   seqinfo: 2 sequences from an unspecified genome

Now that we have identified a set of differentially methylated regions we can proceed to creating the HTML report. Note that this report has less information than the DiffBind example because we don’t have a p-value variable.

    ## Load regionReport
    library('regionReport')

    ## Make it so that the report will be available as a vignette
    original <- readLines(system.file('regionExploration', 'regionExploration.Rmd',
        package = 'regionReport'))
    vignetteInfo <- c(
        'vignette: >',
        '  %\\VignetteEngine{knitr::rmarkdown}',
        '  %\\VignetteIndexEntry{Basic genomic regions exploration}',
        '  %\\VignetteEncoding{UTF-8}'
    )
    new <- c(original[1:16], vignetteInfo, original[17:length(original)])
    writeLines(new, 'regionReportBumphunter.Rmd')

    ## Now create the report
    report <- renderReport(regions, 'Example bumphunter', pvalueVars = NULL,
        densityVars = c('Area' = 'area', 'Value' = 'value',
        'Cluster Length' = 'clusterL'), significantVar = NULL,
        output = 'bumphunterExampleOutput', outdir = '.',
        template = 'regionReportBumphunter.Rmd', device = 'png')

    ## Writing 11 Bibtex entries ... OK
    ## Results written to file './bumphunterExampleOutput.bib'
    ## processing file: bumphunterExampleOutput.Rmd

    ## Warning: running command 'C:\Windows\system32\cmd.exe /c convert
    ## "bumphunterExampleOutput_files/figure-html/regLen2-1.png" -trim
    ## "bumphunterExampleOutput_files/figure-html/regLen2-1.png"' had status 4

    ## Warning in shell(paste(c(cmd, args), collapse = " ")): 'convert
    ## "bumphunterExampleOutput_files/figure-html/regLen2-1.png" -trim
    ## "bumphunterExampleOutput_files/figure-html/regLen2-1.png"' execution failed
    ## with error code 4

    ## 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.
    ## Scale for 'x' is already present. Adding another scale for 'x', which will
    ## replace the existing scale.

    ## Warning: running command 'C:\Windows\system32\cmd.exe /c convert
    ## "bumphunterExampleOutput_files/figure-html/genomeOverview2-1.png" -trim
    ## "bumphunterExampleOutput_files/figure-html/genomeOverview2-1.png"' had status
    ## 4

    ## Warning in shell(paste(c(cmd, args), collapse = " ")): 'convert
    ## "bumphunterExampleOutput_files/figure-html/genomeOverview2-1.png" -trim
    ## "bumphunterExampleOutput_files/figure-html/genomeOverview2-1.png"' execution
    ## failed with error code 4

    ## Warning: running command 'C:\Windows\system32\cmd.exe /c convert
    ## "bumphunterExampleOutput_files/figure-html/vennDiagram-1.png" -trim
    ## "bumphunterExampleOutput_files/figure-html/vennDiagram-1.png"' had status 4

    ## Warning in shell(paste(c(cmd, args), collapse = " ")): 'convert
    ## "bumphunterExampleOutput_files/figure-html/vennDiagram-1.png" -trim
    ## "bumphunterExampleOutput_files/figure-html/vennDiagram-1.png"' execution
    ## failed with error code 4

    ## output file: bumphunterExampleOutput.knit.md

    ## "C:/PROGRA~2/Pandoc/pandoc" +RTS -K512m -RTS bumphunterExampleOutput.utf8.md --to html --from markdown+autolink_bare_uris+ascii_identifiers+tex_math_single_backslash+smart --output bumphunterExampleOutput.html --email-obfuscation none --self-contained --standalone --section-divs --table-of-contents --toc-depth 3 --variable toc_float=1 --variable toc_selectors=h1,h2,h3 --variable toc_collapsed=1 --variable toc_smooth_scroll=1 --variable toc_print=1 --template "C:\Users\biocbuild\bbs-3.6-bioc\tmpdir\RtmpK8ganD/BiocStyle/template.html" --no-highlight --variable highlightjs=1 --number-sections --css "C:\Users\biocbuild\bbs-3.6-bioc\R\library\BiocStyle\resources\html\bioconductor2.css" --variable "theme:bootstrap" --include-in-header "C:\Users\biocbuild\bbs-3.6-bioc\tmpdir\RtmpK8ganD\rmarkdown-str11f49439e.html" --mathjax --variable "mathjax-url:https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" --variable code_folding=hide --variable code_menu=1

    ##
    ## Output created: bumphunterExampleOutput.html

    ## Clean up
    file.remove('regionReportBumphunter.Rmd')

    ## [1] TRUE

You can view the final report here.

In case the link does not work, a pre-compiled version of this document and its corresponding report are available at leekgroup.github.io/regionReportSupp/.

# 2 Reproducibility

    ## Date generated:
    Sys.time()

    ## [1] "2018-01-20 20:53:03 EST"

    ## Time spent making this page:
    proc.time()

    ##    user  system elapsed
    ##  163.26    4.68  171.09

    ## R and packages info:
    options(width = 120)
    devtools::session_info()

    ## Session info ----------------------------------------------------------------------------------------------------------

    ##  setting  value
    ##  version  R version 3.4.3 (2017-11-30)
    ##  system   x86_64, mingw32
    ##  ui       RTerm
    ##  language (EN)
    ##  collate  C
    ##  tz       America/New_York
    ##  date     2018-01-20

    ## Packages --------------------------------------------------------------------------------------------------------------

    ##  package                           * version   date       source
    ##  acepack                             1.4.1     2016-10-29 CRAN (R 3.4.3)
    ##  annotate                            1.56.1    2018-01-20 Bioconductor
    ##  AnnotationDbi                     * 1.40.0    2018-01-20 Bioconductor
    ##  AnnotationFilter                    1.2.0     2018-01-20 Bioconductor
    ##  AnnotationHub                       2.10.1    2018-01-20 Bioconductor
    ##  assertthat                          0.2.0     2017-04-11 CRAN (R 3.4.3)
    ##  backports                           1.1.2     2017-12-13 CRAN (R 3.4.3)
    ##  base                              * 3.4.3     2017-11-30 local
    ##  base64enc                           0.1-3     2015-07-28 CRAN (R 3.4.1)
    ##  bibtex                              0.4.2     2017-06-30 CRAN (R 3.4.3)
    ##  Biobase                           * 2.38.0    2018-01-20 Bioconductor
    ##  BiocGenerics                      * 0.24.0    2018-01-20 Bioconductor
    ##  BiocInstaller                       1.28.0    2018-01-20 Bioconductor
    ##  BiocParallel                        1.12.0    2018-01-20 Bioconductor
    ##  BiocStyle                         * 2.6.1     2018-01-20 Bioconductor
    ##  biomaRt                             2.34.2    2018-01-20 Bioconductor
    ##  Biostrings                          2.46.0    2018-01-20 Bioconductor
    ##  biovizBase                          1.26.0    2018-01-21 Bioconductor
    ##  bit                                 1.1-12    2014-04-09 CRAN (R 3.4.1)
    ##  bit64                               0.9-7     2017-05-08 CRAN (R 3.4.1)
    ##  bitops                              1.0-6     2013-08-17 CRAN (R 3.4.1)
    ##  blob                                1.1.0     2017-06-17 CRAN (R 3.4.3)
    ##  bookdown                            0.5       2017-08-20 CRAN (R 3.4.3)
    ##  BSgenome                            1.46.0    2018-01-20 Bioconductor
    ##  bumphunter                        * 1.20.0    2018-01-21 Bioconductor
    ##  checkmate                           1.8.5     2017-10-24 CRAN (R 3.4.3)
    ##  cluster                             2.0.6     2017-03-10 CRAN (R 3.4.3)
    ##  codetools                           0.2-15    2016-10-05 CRAN (R 3.4.3)
    ##  colorspace                          1.3-2     2016-12-14 CRAN (R 3.4.3)
    ##  compiler                            3.4.3     2017-11-30 local
    ##  curl                                3.1       2017-12-12 CRAN (R 3.4.3)
    ##  data.table                          1.10.4-3  2017-10-27 CRAN (R 3.4.3)
    ##  datasets                          * 3.4.3     2017-11-30 local
    ##  DBI                                 0.7       2017-06-18 CRAN (R 3.4.3)
    ##  DEFormats                           1.6.0     2018-01-20 Bioconductor
    ##  DelayedArray                        0.4.1     2018-01-20 Bioconductor
    ##  derfinder                         * 1.12.6    2018-01-21 Bioconductor
    ##  derfinderHelper                     1.12.0    2018-01-20 Bioconductor
    ##  derfinderPlot                     * 1.12.2    2018-01-21 Bioconductor
    ##  DESeq2                              1.18.1    2018-01-20 Bioconductor
    ##  devtools                          * 1.13.4    2017-11-09 CRAN (R 3.4.3)
    ##  dichromat                           2.0-0     2013-01-24 CRAN (R 3.4.1)
    ##  digest                              0.6.14    2018-01-14 CRAN (R 3.4.3)
    ##  doRNG                               1.6.6     2017-04-10 CRAN (R 3.4.3)
    ##  DT                                * 0.2       2016-08-09 CRAN (R 3.4.3)
    ##  edgeR                               3.20.7    2018-01-20 Bioconductor
    ##  ensembldb                           2.2.0     2018-01-21 Bioconductor
    ##  evaluate                            0.10.1    2017-06-24 CRAN (R 3.4.3)
    ##  foreach                           * 1.4.4     2017-12-12 CRAN (R 3.4.3)
    ##  foreign                             0.8-69    2017-06-22 CRAN (R 3.4.3)
    ##  Formula                             1.2-2     2017-07-10 CRAN (R 3.4.1)
    ##  genefilter                          1.60.0    2018-01-20 Bioconductor
    ##  geneplotter                         1.56.0    2018-01-20 Bioconductor
    ##  GenomeInfoDb                      * 1.14.0    2018-01-20 Bioconductor
    ##  GenomeInfoDbData                    1.0.0     2017-12-14 Bioconductor
    ##  GenomicAlignments                   1.14.1    2018-01-20 Bioconductor
    ##  GenomicFeatures                   * 1.30.0    2018-01-20 Bioconductor
    ##  GenomicFiles                        1.14.0    2018-01-21 Bioconductor
    ##  GenomicRanges                     * 1.30.1    2018-01-20 Bioconductor
    ##  GGally                              1.3.2     2017-08-02 CRAN (R 3.4.3)
    ##  ggbio                             * 1.26.0    2018-01-21 Bioconductor
    ##  ggplot2                           * 2.2.1     2016-12-30 CRAN (R 3.4.3)
    ##  graph                               1.56.0    2018-01-20 Bioconductor
    ##  graphics                          * 3.4.3     2017-11-30 local
    ##  grDevices                         * 3.4.3     2017-11-30 local
    ##  grid                              * 3.4.3     2017-11-30 local
    ##  gridExtra                         * 2.3       2017-09-09 CRAN (R 3.4.3)
    ##  gtable                              0.2.0     2016-02-26 CRAN (R 3.4.3)
    ##  highr                               0.6       2016-05-09 CRAN (R 3.4.3)
    ##  Hmisc                               4.1-1     2018-01-03 CRAN (R 3.4.3)
    ##  htmlTable                           1.11.2    2018-01-20 CRAN (R 3.4.3)
    ##  htmltools                           0.3.6     2017-04-28 CRAN (R 3.4.3)
    ##  htmlwidgets                         1.0       2018-01-20 CRAN (R 3.4.3)
    ##  httpuv                              1.3.5     2017-07-04 CRAN (R 3.4.3)
    ##  httr                                1.3.1     2017-08-20 CRAN (R 3.4.3)
    ##  interactiveDisplayBase              1.16.0    2018-01-20 Bioconductor
    ##  IRanges                           * 2.12.0    2018-01-20 Bioconductor
    ##  iterators                         * 1.0.9     2017-12-12 CRAN (R 3.4.3)
    ##  jsonlite                            1.5       2017-06-01 CRAN (R 3.4.3)
    ##  knitcitations                       1.0.8     2017-07-04 CRAN (R 3.4.3)
    ##  knitr                             * 1.18      2017-12-27 CRAN (R 3.4.3)
    ##  knitrBootstrap                      1.0.1     2017-07-19 CRAN (R 3.4.3)
    ##  labeling                            0.3       2014-08-23 CRAN (R 3.4.1)
    ##  lattice                             0.20-35   2017-03-25 CRAN (R 3.4.3)
    ##  latticeExtra                        0.6-28    2016-02-09 CRAN (R 3.4.3)
    ##  lazyeval                            0.2.1     2017-10-29 CRAN (R 3.4.3)
    ##  limma                               3.34.5    2018-01-20 Bioconductor
    ##  locfit                            * 1.5-9.1   2013-04-20 CRAN (R 3.4.3)
    ##  lubridate                           1.7.1     2017-11-03 CRAN (R 3.4.3)
    ##  magrittr                            1.5       2014-11-22 CRAN (R 3.4.3)
    ##  markdown                            0.8       2017-04-20 CRAN (R 3.4.3)
    ##  Matrix                              1.2-12    2017-11-20 CRAN (R 3.4.3)
    ##  matrixStats                         0.52.2    2017-04-14 CRAN (R 3.4.3)
    ##  memoise                             1.1.0     2017-04-21 CRAN (R 3.4.3)
    ##  methods                           * 3.4.3     2017-11-30 local
    ##  mgcv                              * 1.8-23    2018-01-15 CRAN (R 3.4.3)
    ##  mime                                0.5       2016-07-07 CRAN (R 3.4.1)
    ##  munsell                             0.4.3     2016-02-13 CRAN (R 3.4.3)
    ##  nlme                              * 3.1-131   2017-02-06 CRAN (R 3.4.3)
    ##  nnet                                7.3-12    2016-02-02 CRAN (R 3.4.3)
    ##  org.Hs.eg.db                      * 3.5.0     2017-12-04 Bioconductor
    ##  OrganismDbi                         1.20.0    2018-01-21 Bioconductor
    ##  parallel                          * 3.4.3     2017-11-30 local
    ##  pillar                              1.1.0     2018-01-14 CRAN (R 3.4.3)
    ##  pkgconfig                           2.0.1     2017-03-21 CRAN (R 3.4.3)
    ##  pkgmaker                            0.22      2014-05-14 CRAN (R 3.4.3)
    ##  plyr                                1.8.4     2016-06-08 CRAN (R 3.4.3)
    ##  prettyunits                         1.0.2     2015-07-13 CRAN (R 3.4.3)
    ##  progress                            1.1.2     2016-12-14 CRAN (R 3.4.3)
    ##  ProtGenerics                        1.10.0    2018-01-20 Bioconductor
    ##  qvalue                              2.10.0    2018-01-20 Bioconductor
    ##  R6                                  2.2.2     2017-06-17 CRAN (R 3.4.3)
    ##  RBGL                                1.54.0    2018-01-20 Bioconductor
    ##  RColorBrewer                      * 1.1-2     2014-12-07 CRAN (R 3.4.1)
    ##  Rcpp                                0.12.15   2018-01-20 CRAN (R 3.4.3)
    ##  RCurl                               1.95-4.10 2018-01-04 CRAN (R 3.4.3)
    ##  RefManageR                          0.14.20   2017-08-17 CRAN (R 3.4.3)
    ##  regionReport                      * 1.12.2    2018-01-21 Bioconductor
    ##  registry                            0.5       2017-12-03 CRAN (R 3.4.3)
    ##  reshape                             0.8.7     2017-08-06 CRAN (R 3.4.3)
    ##  reshape2                            1.4.3     2017-12-11 CRAN (R 3.4.3)
    ##  rlang                               0.1.6     2017-12-21 CRAN (R 3.4.3)
    ##  rmarkdown                           1.8       2017-11-17 CRAN (R 3.4.3)
    ##  RMySQL                              0.10.13   2017-08-14 CRAN (R 3.4.3)
    ##  rngtools                            1.2.4     2014-03-06 CRAN (R 3.4.3)
    ##  rpart                               4.1-12    2018-01-12 CRAN (R 3.4.3)
    ##  rprojroot                           1.3-2     2018-01-03 CRAN (R 3.4.3)
    ##  Rsamtools                           1.30.0    2018-01-20 Bioconductor
    ##  RSQLite                             2.0       2017-06-19 CRAN (R 3.4.1)
    ##  rstudioapi                          0.7       2017-09-07 CRAN (R 3.4.3)
    ##  rtracklayer                         1.38.2    2018-01-20 Bioconductor
    ##  S4Vectors                         * 0.16.0    2018-01-20 Bioconductor
    ##  scales                              0.5.0     2017-08-24 CRAN (R 3.4.3)
    ##  shiny                               1.0.5     2017-08-23 CRAN (R 3.4.3)
    ##  splines                             3.4.3     2017-11-30 local
    ##  stats                             * 3.4.3     2017-11-30 local
    ##  stats4                            * 3.4.3     2017-11-30 local
    ##  stringi                             1.1.6     2017-11-17 CRAN (R 3.4.2)
    ##  stringr                             1.2.0     2017-02-18 CRAN (R 3.4.3)
    ##  SummarizedExperiment                1.8.1     2018-01-20 Bioconductor
    ##  survival                            2.41-3    2017-04-04 CRAN (R 3.4.3)
    ##  tibble                              1.4.1     2017-12-25 CRAN (R 3.4.3)
    ##  tools                               3.4.3     2017-11-30 local
    ##  TxDb.Hsapiens.UCSC.hg19.knownGene * 3.2.2     2017-12-05 Bioconductor
    ##  utils                             * 3.4.3     2017-11-30 local
    ##  VariantAnnotation                   1.24.5    2018-01-21 Bioconductor
    ##  whisker                           * 0.3-2     2013-04-28 CRAN (R 3.4.3)
    ##  withr                               2.1.1     2017-12-19 CRAN (R 3.4.3)
    ##  XML                                 3.98-1.9  2017-06-19 CRAN (R 3.4.1)
    ##  xml2                                1.1.1     2017-01-24 CRAN (R 3.4.3)
    ##  xtable                              1.8-2     2016-02-05 CRAN (R 3.4.3)
    ##  XVector                             0.18.0    2018-01-20 Bioconductor
    ##  yaml                                2.1.16    2017-12-12 CRAN (R 3.4.3)
    ##  zlibbioc                            1.24.0    2018-01-20 Bioconductor