if (require(minfiData) & require(GenomicRanges)){
set <- prepareMethylationSet(getBeta(MsetEx), pheno = pData(MsetEx))
range <- GenomicRanges::GRanges(seqnames=Rle("chrY"),
ranges = IRanges(3000000, end=12300000))
rangeNoSNPs <- DARegionAnalysis(set, variable_names = "sex", range = range)
plotRegion(rangeNoSNPs)
}
Run the code above in your browser using DataLab