Get the mean methylation difference across a specified region.

Using a chromDictMeth() output object, quickly calculate the mean methylation difference across a user-specified region.

MethylDiff(chromDictMeth, chrom, start, end)

Arguments

chromDictMeth

chromDictMeth() output object

chrom

character, chromosome

start

integer, region start

end

integer, region end

Value

Returns a list of data.tables for each chromosome, for faster analysis.

Examples

chromDictMethylDiff <- chromDictMeth(eigen, c("trt", "ctl"))
#> Removing low depth sites...
#> Calculating mean methylation differences...
#> Splitting by chromosome...
MethylDiff(chromDictMethylDiff, "chr3", 4920450, 4923267)
#> [1] -0.5493768