recountmethylation 1.2.3
The recountmethylation package provides access to databases of DNA
methylation (DNAm) data from over 62,000 cumulative sample records with
IDATs in the Gene Expression Omnibus (GEO, available by November, 2020).
Samples were run using either of 2 Illumina BeadArray platform types,
either the older HM450K platform or the newer EPIC/HM850K platform. The
database compilation files include mined, mapped, and model-based sample
metadata, and DNAm data in the form of either raw/unnormalized red and
green signals, raw/unnormalized methylated and unmethylated signals, or
normalized DNAm fractions (a.k.a. “Beta-values”) [1].
Normalization was performed using the out-of-band signal correction (a.k.a.
“noob”) method, a type of within-sample normalization [2].
This User’s Guide shows how to use the recountmethylation package to obtain,
load, and query the DNAm databases with 2 small example files. Background about
DNAm arrays, DNAm measurement, SummarizedExperiment objects, database file
types, and samples metadata is also provided. Further analysis examples are
contained in the data_analyses vignette.
Databases accessed with recountmethylation contain data from GEO
(ncbi.nlm.nih.gov/geo/), a live public database where alterations to
online records can cause discrepancies with stored data over time.
We cannot guarantee the accuracy of stored data, and advise users
cross-check their findings with latest available records.
Database compilation file download and access is managed by the get_db
functions, where the DNAm array platform type using the platform argument
(see ?get_db for details). Both HM450K and EPIC/HM850K platforms are
currently supported (see below for platform details). Note you will need
between 50-180 Gb of disk space to store a single database file. Files pair
sample metadata and assay data in various formats, including HDF5-SummarizedExperiment database directories, and HDF5 database files
with the .h5 extension.
The databases are located at https://methylation.recount.bio/, and file details are viewable as follows:
sm <- as.data.frame(smfilt(get_servermatrix()))
if(is(sm, "data.frame")){knitr::kable(sm, align = "c")}| filename | date | time | size (bytes) |
|---|---|---|---|
| remethdb_h5-rg_epic_0-0-2_1589820348.h5 | 08-Jan-2021 | 09:46 | 66751358297 |
| remethdb_h5se-gm_epic_0-0-2_1589820348 | 06-Jan-2021 | 09:56 | assays.h5 = 56956363488;se.rds = 8475111 |
| remethdb_h5se-gr_epic_0-0-2_1607018051 | 06-Jan-2021 | 10:09 | assays.h5 = 82090895411;se.rds = 8475201 |
| remethdb_h5se-rg_epic_0-0-2_1589820348 | 06-Jan-2021 | 10:19 | assays.h5 = 68707689800;se.rds = 3059883 |
| remethdb_h5-rg_hm450k_0-0-2_1607018051.h5 | 07-Jan-2021 | 10:01 | 193342823766 |
| remethdb_h5se-gm_hm450k_0-0-2_1607018051 | 07-Jan-2021 | 10:23 | assays.h5 = 130935841655;se.rds = 5372091 |
| remethdb_h5se-gr_hm450k_0-0-2_1607018051 | 07-Jan-2021 | 10:50 | assays.h5 = 184355830172;se.rds = 5372008 |
| remethdb_h5se-rg_hm450k_0-0-2_1607018051 | 06-Jan-2021 | 10:44 | assays.h5 = 164788908310;se.rds = 3179962 |
| remethdb-h5se_gr-test_0-0-1_1590090412 | 29-May-2020 | 07:28 | assays.h5 = 132596;se.rds = 68522 |
| remethdb-h5_rg-test_0-0-1_1590090412.h5 | 31-May-2020 | 07:26 | 252757 |
The DNAm array database files are indexed on ExperimentHub, and are
viewable as follows.
hub = ExperimentHub::ExperimentHub() # connect to the hubs
rmdat <- AnnotationHub::query(hub, "recountmethylation") # query the hubs
rmdat## ExperimentHub with 6 records
## # snapshotDate(): 2021-05-18
## # $dataprovider: GEO/GDS
## # $species: Homo sapiens
## # $rdataclass: HDF5-SummarizedExperiment, HDF5Database
## # additional mcols(): taxonomyid, genome, description,
## # coordinate_1_based, maintainer, rdatadateadded, preparerclass, tags,
## # rdatapath, sourceurl, sourcetype
## # retrieve records with, e.g., 'object[["EH3773"]]'
##
## title
## EH3773 | remethdb-h5se_gr_0-0-1_1590090412
## EH3774 | remethdb-h5se_gm_0-0-1_1590090412
## EH3775 | remethdb_h5se-rg_0-0-1_1590090412
## EH3776 | remethdb-h5se_gr-test_0-0-1_159009041
## EH3777 | remethdb-h5_rg_0-0-1_1590090412.h5
## EH3778 | remethdb-h5_rg-test_0-0-1_1590090412.h5In addition to using the getdb functions, the HDF5 (“.h5”" extension)
files may alternatively be downloaded from the hubs, as follows.
eid <- "EH3778" # h5 test file id
fpath <- rmdat[[eid]] # download with default caching
rhdf5::h5ls(fpath)## group name otype dclass dim
## 0 / greensignal H5I_DATASET FLOAT 2 x 11162
## 1 / greensignal.colnames H5I_DATASET STRING 11162
## 2 / greensignal.rownames H5I_DATASET STRING 2
## 3 / mdpost H5I_DATASET STRING 2 x 19
## 4 / mdpost.colnames H5I_DATASET STRING 19
## 5 / redsignal H5I_DATASET FLOAT 2 x 11162
## 6 / redsignal.colnames H5I_DATASET STRING 11162
## 7 / redsignal.rownames H5I_DATASET STRING 2Note that whether downloads use the hubs or getdb functions, caching
is implemented to check for previously downloaded database files.
This section includes essential background about DNAm array platforms, assays and file types, and sample metadata.
Databases include human samples run on the Illumina Infinium HM450K BeadArray platform. HM450K is a popular 2-channel platform that probes over 480,000 CpG loci genome-wide, with enriched coverage at CG islands, genes, and enhancers [3]. The more recently released EPIC/HM850K platform contains an expanded probe set targeting over 850,000 CpGs, including more than 90% of the HM450K probes, with greater coverage of potential intergenic regulatory regions [4].
Array processing generates 2 intensity files (IDATs) per sample, one each for
the red and green color channels. These raw files also contain control signals
useful for quality evaluations [5]. The BeadArray probes use
either of 2 bead technologies, known as Type I and Type II, where the majority
(72%) of probes use the latter. For Type II probes, a single bead assay informs
a single probe, while Type I probes use 2 beads each. Practically, this means
the bead-specific matrices found in RGChannelSet objects are larger than the
probe-specific matrices found in derived object types (e.g. for HM450K samples,
622,399 assays for red/green signal matrices versus 485,512 assays for
methylated/unmethylated signal, DNAm fractions matrices, see below).
SummarizedExperiment object classesDNAm array sample IDATs can be read into an R session as an object of class
RGChannelSet, a type of SummarizedExperiment. These objects support
analyses of high-throughput genomics datasets, and they include slots for
assay matrices, sample metadata, and experiment metadata. During a typical
workflow, normalization and preprocessing convert RGChannelSet objects into
new types like MethylSet and RatioSet. While not all IDAT information is
accessible from every object type (e.g. only RGChannelSets can contain
control assays), derived objects like MethylSets and RatioSets may be
smaller and/or faster to access.
Three SummarizedExperiment databases are provided as
HDF5-SummarizedExperiment files, including an unnormalized RGChannelSet
(red/green signals), an unnormalized MethylSet (methylated/unmethylated
signals) and a normalized GenomicRatioSet (DNAm fractions). For the latter,
DNAm fractions (logit2 Beta-values, or M-values) were normalized using the
out-of-band signal or “noob” method, an effective within-sample normalization
that removes signal artifacts [2].
Database files are stored as either HDF5 or HDF5-SummarizedExperiment. For
most R users, the latter files will be most convenient to work with. HDF5, or
hierarchical data format 5, combines compression and chunking for convenient
handling of large datasets. HDF5-SummarizedExperiment files combine the
benefits of HDF5 and SummarizedExperiment entities using a
DelayedArray-powered backend. Once an HDF5-SummarizedExperiment file is
loaded, it can be treated similarly to a SummarizedExperiment object in
active memory. That is, summary and subset operations execute rapidly, and
realization of large data chunks in active memory is delayed until called for
by the script (see examples).
Sample metadata are included with DNAm assays in the database files. Currently,
metadata variables include GEO record IDs for samples (GSM) and studies (GSE),
sample record titles, learned labels for tissue and disease, sample type
predictions from the MetaSRA-pipeline, and DNAm model-based predictions for
age, sex, and blood cell types. Access sample metadata from
SummarizedExperiment objects using the pData minfi function (see examples).
Examples in the data_analyses vignette illustrate some ways to utilize the
provided sample metadata.
Provided metadata derives from the GSE-specific SOFT files, which contain
experiment, sample, and platform metadata. Considerable efforts were made to
learn, harmonize, and predict metadata labels. Certain types of info lacking
in the recountmethylation metadata may be available in the SOFT files,
especially if it is sample non-specific (e.g. methods text, PubMed ID, etc.)
or redundant with DNAm-derived metrics (e.g. DNAm summaries, predicted sex,
etc.).
It is good practice to validate the harmonized metadata with original metadata records, especially where labels are ambiguous or there is insufficient information for a given query. GEO GSM and GSE records can be viewed from a browser, or SOFT files may be downloaded directly. Packages like GEOmetadb and GEOquery are also useful to query and summarize GEO metadata.
HDF5-SummarizedExperiment exampleThis example shows basic handling for HDF5-SummarizedExperiment (a.k.a.
“h5se”) files. For these files, the getdb function returns the loaded file.
Thanks to a DelayedArray backend, even full-sized h5se databases can be
treated as if they were fully loaded into active memory.
The test h5se dataset includes sample metadata and noob-normalized
DNAm fractions (Beta-values) for chromosome 22 probes for 2 samples.
Datasets can be downloaded using the getdb series of functions
(see ?getdb for details), where the dfp argument specifies the
download destination. The test h5se file is included in the package
“inst” directory, and can be loaded as follows.
dn <- "remethdb-h5se_gr-test_0-0-1_1590090412"
path <- system.file("extdata", dn, package = "recountmethylation")
h5se.test <- HDF5Array::loadHDF5SummarizedExperiment(path)Common characterization functions can be used on the dataset after it has been
loaded. These include functions for SummarizedExperiment-like objects, such
as the getBeta, pData, and getAnnotation minfi functions. First, inspect
the dataset using standard functions like class, dim, and summary as
follows.
class(h5se.test)## [1] "GenomicRatioSet"
## attr(,"package")
## [1] "minfi"dim(h5se.test)## [1] 8552 2summary(h5se.test)## [1] "GenomicRatioSet object of length 8552 with 0 metadata columns"Access the sample metadata for the 2 available samples using pData.
h5se.md <- minfi::pData(h5se.test)
dim(h5se.md)## [1] 2 19colnames(h5se.md)## [1] "gsm" "gsm_title" "gseid" "disease"
## [5] "tissue" "sampletype" "arrayid_full" "basename"
## [9] "age" "predage" "sex" "predsex"
## [13] "predcell.CD8T" "predcell.CD4T" "predcell.NK" "predcell.Bcell"
## [17] "predcell.Mono" "predcell.Gran" "storage"Next get CpG probe-specific DNAm fractions, or “Beta-values”, with getBeta
(rows are probes, columns are samples).
h5se.bm <- minfi::getBeta(h5se.test)
dim(h5se.bm)## [1] 8552 2colnames(h5se.bm) <- h5se.test$gsm
knitr::kable(head(h5se.bm), align = "c")| GSM1038308 | GSM1038309 | |
|---|---|---|
| cg00017461 | 0.9807283 | 0.9746836 |
| cg00077299 | 0.3476970 | 0.3456837 |
| cg00079563 | 0.8744652 | 0.9168005 |
| cg00087182 | 0.9763206 | 0.9760947 |
| cg00093544 | 0.0225112 | 0.0265087 |
| cg00101350 | 0.9736359 | 0.9789818 |
Access manifest information for probes with getAnnotation. This includes the
bead addresses, probe type, and genome coordinates and regions. For full details
about the probe annotations, consult the minfi and Illumina platform documentation.
an <- minfi::getAnnotation(h5se.test)
dim(an)## [1] 8552 33colnames(an)## [1] "chr" "pos"
## [3] "strand" "Name"
## [5] "AddressA" "AddressB"
## [7] "ProbeSeqA" "ProbeSeqB"
## [9] "Type" "NextBase"
## [11] "Color" "Probe_rs"
## [13] "Probe_maf" "CpG_rs"
## [15] "CpG_maf" "SBE_rs"
## [17] "SBE_maf" "Islands_Name"
## [19] "Relation_to_Island" "Forward_Sequence"
## [21] "SourceSeq" "Random_Loci"
## [23] "Methyl27_Loci" "UCSC_RefGene_Name"
## [25] "UCSC_RefGene_Accession" "UCSC_RefGene_Group"
## [27] "Phantom" "DMR"
## [29] "Enhancer" "HMM_Island"
## [31] "Regulatory_Feature_Name" "Regulatory_Feature_Group"
## [33] "DHS"ant <- as.matrix(t(an[c(1:4), c(1:3, 5:6, 9, 19, 24, 26)]))
knitr::kable(ant, align = "c")| cg00017461 | cg00077299 | cg00079563 | cg00087182 | |
|---|---|---|---|---|
| chr | chr22 | chr22 | chr22 | chr22 |
| pos | 30663316 | 18632618 | 43253521 | 24302043 |
| strand | - | + | + | + |
| AddressA | 31616369 | 13618325 | 65630302 | 37797387 |
| AddressB | 70798487 | 37626331 | 55610348 | 20767312 |
| Type | I | I | I | I |
| Relation_to_Island | OpenSea | N_Shore | Island | N_Shore |
| UCSC_RefGene_Name | OSM | USP18 | ARFGAP3;ARFGAP3 | GSTT2B;GSTT2 |
| UCSC_RefGene_Group | TSS1500 | TSS200 | TSS200;TSS200 | Body;Body |
HDF5 database and exampleTo provide more workflow options, bead-specific red and green signal data have
been provided with sample metadata in an HDF5/h5 file. This example shows
how to handle objects of this type with recountmethylation.
The test h5 file includes metadata and bead-specific signals from
chromosome 22 for the same 2 samples as in the h5se test file.
Note getdb functions for h5 files simply return the database path.
Since the test h5 file has also been included in the package “inst” folder,
get the path to load the file as follows.
dn <- "remethdb-h5_rg-test_0-0-1_1590090412.h5"
h5.test <- system.file("extdata", "h5test", dn,
package = "recountmethylation")Use the file path to read data into an RGChannelSet with the getrg
function. Setting all.gsm = TRUE obtains data for all samples in the
database files, while passing a vector of GSM IDs to gsmv argument
will query a subset of available samples. Signals from all available
probes are retrieved by default, and probe subsets can be obtained by
passing a vector of valid bead addresses to the cgv argument.
h5.rg <- getrg(dbn = h5.test, all.gsm = TRUE)To avoid exhausting active memory with the full-sized h5 dataset, provide
either gsmv or cgv to getrg, and set either all.cg or all.gsm to
FALSE (see ?getrg for details).
As in the previous example, use pData and getAnnotation to get sample
metadata and array manifest information, respectively. Access the green and
red signal matrices in the RGChannelSet with the getRed and getGreen
minfi functions.
h5.red <- minfi::getRed(h5.rg)
h5.green <- minfi::getGreen(h5.rg)
dim(h5.red)## [1] 11162 2knitr::kable(head(h5.red), align = "c")| GSM1038308 | GSM1038309 | |
|---|---|---|
| 10601475 | 1234 | 1603 |
| 10603366 | 342 | 344 |
| 10603418 | 768 | 963 |
| 10605304 | 2368 | 2407 |
| 10605460 | 3003 | 3322 |
| 10608343 | 357 | 399 |
knitr::kable(head(h5.green), align = "c")| GSM1038308 | GSM1038309 | |
|---|---|---|
| 10601475 | 6732 | 8119 |
| 10603366 | 288 | 356 |
| 10603418 | 267 | 452 |
| 10605304 | 4136 | 4395 |
| 10605460 | 1395 | 1762 |
| 10608343 | 840 | 1269 |
identical(rownames(h5.red), rownames(h5.green))## [1] TRUERows in these signal matrices map to bead addresses rather than probe IDs.
These matrices have more rows than the h5se test Beta-value matrix because
any type I probes use data from 2 beads each.
This section demonstrates validation using the test databases. Full code to reproduce this section is provided but not evaluated, as it involves a download from the GEO servers. As the disclaimer notes, it is good practice to validate data against the latest available GEO files. This step may be most useful for newer samples published close to the end compilation date (through November 7, 2020 for current version), which may be more prone to revisions at initial publication.
Use the gds_idat2rg function to download IDATs for the 2 test samples
and load these into a new RGChannelSet object. Do this by passing a vector
of GSM IDs to gsmv and the download destination to dfp. (note, chunks in
this section are fully executable, but not evaluated for this vignette).
# download from GEO
dlpath <- tempdir()
gsmv <- c("GSM1038308", "GSM1038309")
geo.rg <- gds_idat2rg(gsmv, dfp = dlpath)
colnames(geo.rg) <- gsub("\\_.*", "", colnames(geo.rg))Extract the red and green signal matrices from geo.rg.
geo.red <- minfi::getRed(geo.rg)
geo.green <- minfi::getGreen(geo.rg)Match indices and labels between the GEO and h5 test signal matrices.
int.addr <- intersect(rownames(geo.red), rownames(h5.red))
geo.red <- geo.red[int.addr,]
geo.green <- geo.green[int.addr,]
geo.red <- geo.red[order(match(rownames(geo.red), rownames(h5.red))),]
geo.green <- geo.green[order(match(rownames(geo.green), rownames(h5.green))),]
identical(rownames(geo.red), rownames(h5.red))
identical(rownames(geo.green), rownames(h5.green))
class(h5.red) <- "integer"
class(h5.green) <- "integer"Finally, compare the signal matrix data.
identical(geo.red, h5.red)identical(geo.green, h5.green)Before comparing the GEO-downloaded data to data from the h5se.test database,
normalize the data using the same out-of-band or “noob” normalization technique
that was used to generate data in the h5se database.
geo.gr <- minfi::preprocessNoob(geo.rg)Next, extract the Beta-values.
geo.bm <- as.matrix(minfi::getBeta(geo.gr))Now match row and column labels and indices.
h5se.bm <- as.matrix(h5se.bm)
int.cg <- intersect(rownames(geo.bm), rownames(h5se.bm))
geo.bm <- geo.bm[int.cg,]
geo.bm <- geo.bm[order(match(rownames(geo.bm), rownames(h5se.bm))),]Finally, compare the two datasets.
identical(summary(geo.bm), summary(h5se.bm))identical(rownames(geo.bm), rownames(h5se.bm))This section describes how to address potential issues with accessing the
database files or working with the DelayedArray based objects locally.
If repeated attempts to download the database compilation files fail, you may try the following:
First ensure your internet connection is stable and there is sufficient space at the download destination for the database file.
Second, try increasing your timeout duration beyond the default before
repeating the download attempt with getdb. Check the current timeout
for an R session with getOptions('timeout'), then manually increase
the timeout duration with options(timeout = new.time).
Finally, you may attempt to download a server file using command line
calls to your system terminal or console. For instance, on a Mac you
might try wget -r <file_url>. If this doesn’t work, you can again
attempt to increase the timeout duration and repeat the download attempt.
DelayedArray inputsUnexpected function behaviors may arise when using DelayedArray-based inputs.
These essentially arise from lacking interoperativity between normal matrices
and the DelayedArray-based matrices. Known examples include:
minfi::detectionP():Throws error for specific subsets of data, such as for queries of exactly 50 samples.
detectionP(rg[,1:50])
"Error in .local(Red, Green, locusNames, controlIdx, TypeI.Red, TypeI.Green, dim(Red_grid) == dim(detP_sink_grid) are not all TRUE"minfi::preprocessFunnorm():Throws error when called for an RGChannelSet of type HDF5-SummarizedExperiment.
preprocessFunnorm(rg)
"Error: 'preprocessFunnorm()' only supports matrix-backed minfi objects.""These and other related errors may be addressed by instantiating the data query,
or the data chunk, as a new non-DelayedArray object. For example, remake a
subset of the full h5se dataset, rg, as follows.
rg.h5se <- loadHDF5SummarizedExperiment(rg.path) # full h5se RGChannelSet
rg.sub <- rg.h5se[,c(1:20)] # subset samples of interest
# re-make as non-DA object
rg.new <- RGChannelSet(Red = getRed(rg.sub),
Green = getGreen(rg.sub),
annotation = annotation(rg.sub))
gr <- preprocessFunnorm(rg.new) # repeat preprocessingAlternatively, non-DelayedArray RGChannelSet objects can be readily generated from
the full h5 RGChannelSet database with the provided function getrg().
Consult the Data Analyses vignette and main manuscript for analysis examples and details about data compilations.
sessionInfo()## R version 4.1.1 (2021-08-10)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.13-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.13-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] stats4 parallel stats graphics grDevices utils datasets
## [8] methods base
##
## other attached packages:
## [1] IlluminaHumanMethylation450kanno.ilmn12.hg19_0.6.0
## [2] ExperimentHub_2.0.0
## [3] AnnotationHub_3.0.1
## [4] BiocFileCache_2.0.0
## [5] dbplyr_2.1.1
## [6] HDF5Array_1.20.0
## [7] DelayedArray_0.18.0
## [8] Matrix_1.3-4
## [9] limma_3.48.3
## [10] gridExtra_2.3
## [11] ggplot2_3.3.5
## [12] knitr_1.36
## [13] recountmethylation_1.2.3
## [14] minfi_1.38.0
## [15] bumphunter_1.34.0
## [16] locfit_1.5-9.4
## [17] iterators_1.0.13
## [18] foreach_1.5.1
## [19] Biostrings_2.60.2
## [20] XVector_0.32.0
## [21] SummarizedExperiment_1.22.0
## [22] Biobase_2.52.0
## [23] MatrixGenerics_1.4.3
## [24] matrixStats_0.61.0
## [25] GenomicRanges_1.44.0
## [26] GenomeInfoDb_1.28.4
## [27] IRanges_2.26.0
## [28] S4Vectors_0.30.2
## [29] BiocGenerics_0.38.0
## [30] rhdf5_2.36.0
## [31] BiocStyle_2.20.2
##
## loaded via a namespace (and not attached):
## [1] plyr_1.8.6 splines_4.1.1
## [3] BiocParallel_1.26.2 digest_0.6.28
## [5] htmltools_0.5.2 magick_2.7.3
## [7] fansi_0.5.0 magrittr_2.0.1
## [9] memoise_2.0.0 tzdb_0.1.2
## [11] readr_2.0.2 annotate_1.70.0
## [13] askpass_1.1 siggenes_1.66.0
## [15] prettyunits_1.1.1 colorspace_2.0-2
## [17] blob_1.2.2 rappdirs_0.3.3
## [19] xfun_0.26 dplyr_1.0.7
## [21] crayon_1.4.1 RCurl_1.98-1.5
## [23] jsonlite_1.7.2 genefilter_1.74.1
## [25] GEOquery_2.60.0 survival_3.2-13
## [27] glue_1.4.2 gtable_0.3.0
## [29] zlibbioc_1.38.0 Rhdf5lib_1.14.2
## [31] scales_1.1.1 DBI_1.1.1
## [33] rngtools_1.5.2 Rcpp_1.0.7
## [35] xtable_1.8-4 progress_1.2.2
## [37] bit_4.0.4 mclust_5.4.7
## [39] preprocessCore_1.54.0 httr_1.4.2
## [41] RColorBrewer_1.1-2 ellipsis_0.3.2
## [43] farver_2.1.0 pkgconfig_2.0.3
## [45] reshape_0.8.8 XML_3.99-0.8
## [47] sass_0.4.0 utf8_1.2.2
## [49] later_1.3.0 labeling_0.4.2
## [51] tidyselect_1.1.1 rlang_0.4.11
## [53] AnnotationDbi_1.54.1 BiocVersion_3.13.1
## [55] munsell_0.5.0 tools_4.1.1
## [57] cachem_1.0.6 generics_0.1.0
## [59] RSQLite_2.2.8 evaluate_0.14
## [61] stringr_1.4.0 fastmap_1.1.0
## [63] yaml_2.2.1 bit64_4.0.5
## [65] beanplot_1.2 scrime_1.3.5
## [67] purrr_0.3.4 KEGGREST_1.32.0
## [69] nlme_3.1-153 doRNG_1.8.2
## [71] sparseMatrixStats_1.4.2 mime_0.12
## [73] nor1mix_1.3-0 xml2_1.3.2
## [75] biomaRt_2.48.3 compiler_4.1.1
## [77] interactiveDisplayBase_1.30.0 filelock_1.0.2
## [79] curl_4.3.2 png_0.1-7
## [81] tibble_3.1.5 bslib_0.3.1
## [83] stringi_1.7.5 highr_0.9
## [85] GenomicFeatures_1.44.2 lattice_0.20-45
## [87] multtest_2.48.0 vctrs_0.3.8
## [89] pillar_1.6.3 lifecycle_1.0.1
## [91] rhdf5filters_1.4.0 BiocManager_1.30.16
## [93] jquerylib_0.1.4 data.table_1.14.2
## [95] bitops_1.0-7 httpuv_1.6.3
## [97] rtracklayer_1.52.1 R6_2.5.1
## [99] BiocIO_1.2.0 promises_1.2.0.1
## [101] bookdown_0.24 codetools_0.2-18
## [103] MASS_7.3-54 assertthat_0.2.1
## [105] openssl_1.4.5 rjson_0.2.20
## [107] withr_2.4.2 GenomicAlignments_1.28.0
## [109] Rsamtools_2.8.0 GenomeInfoDbData_1.2.6
## [111] mgcv_1.8-38 hms_1.1.1
## [113] quadprog_1.5-8 grid_4.1.1
## [115] tidyr_1.1.4 base64_2.0
## [117] rmarkdown_2.11 DelayedMatrixStats_1.14.3
## [119] illuminaio_0.34.0 shiny_1.7.1
## [121] restfulr_0.0.131. Maden, S. K., Thompson, R. F., Hansen, K. D., and Nellore, A. (2020). Human methylome variation across Infinium 450K data on the Gene Expression Omnibus. bioRxiv. Available at: https://www.biorxiv.org/content/10.1101/2020.11.17.387548v2.
2. Triche, T. J., Weisenberger, D. J., Van Den Berg, D., Laird, P. W., and Siegmund, K. D. (2013). Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucleic Acids Research 41, e90.
3. Sandoval, J., Heyn, H. A., Moran, S., Serra-Musach, J., Pujana, M. A., Bibikova, M., and Esteller, M. (2011). Validation of a DNA methylation microarry for 450,000 CpG sites in the human genome. Epigenetics 6, 692–702. Available at: http://www.landesbioscience.com/journals/epigenetics/article/16196/?nocache=1384341162.
4. Pidsley, R., Zotenko, E., Peters, T. J., Lawrence, M. G., Risbridger, G. P., Molloy, P., Van Djik, S., Muhlhausler, B., Stirzaker, C., and Clark, S. J. (2016). Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biology 17. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5055731/ [Accessed April 19, 2019].