compartments: Estimates A/B compartments from Illumina methylation arrays
In minfi: Analyze Illumina Infinium DNA methylation arrays

Description Usage Arguments Details Value Author(s) References Examples

Estimates A/B compartments as revealed by Hi-C by computing the first eigenvector on a binned probe correlation matrix.

1 2	compartments(object, resolution=100*1000, what = "OpenSea", chr="chr22", method = c("pearson", "spearman"), keep=TRUE)

`object`	An object of class (Genomic)MethylSet or (Genomic)RatioSet
`resolution`	An integer specifying the binning resolution
`what`	Which subset of probes should be used?
`chr`	The chromosome to be analyzed.
`method`	Method of correlation.
`keep`	Should the correlation matrix be stored or not?

This function extracts A/B compartments from Illumina methylation microarrays. Analysis of Hi-C data has shown that the genome can be divided into two compartments (A/B compartments) that are cell-type specific and are associated with open and closed chromatin respectively. The approximately 170,000 open sea probes on the 450k array can be used to estimate these compartments by computing the first eigenvector on a binned correlation matrix. The binning resolution can be specified by resolution, and by default is set to a 100 kb. We do not recommend higher resolutions because of the low-resolution probe design of the 450k array.

an object of class GRanges containing the correlation matrix, the compartment eigenvector and the compartment labels (A or B) as metadata.

Jean-Philippe Fortin jfortin@jhsph.edu, Kasper D. Hansen kasperdanielhansen@gmail.com

JP Fortin and KD Hansen. Reconstructing A/B compartments as revealed by Hi-C using long-range correlations in epigenetic data. bioRxiv (2015). doi:10.1101/019000.

if (require(minfiData)) {
  GMset <- mapToGenome(MsetEx)
  ## compartments at 1MB resolution; we recommend 100kb.
  comps <- compartments(GMset, res = 10^6)
}