Question: How to do DMR analysis with BiSeq (for RRBS methylation data)?
gravatar for ahmad.moousavi
16 days ago by
Royan Institute, Tehran, Iran
ahmad.moousavi40 wrote:


I have 20 RRBS samples which I wanted to do do deferentially methylated region (DMR) analysis with BiSeq or other R packages, however the manual resource is not enough. I want to do DMR analysis on all chromosomes not specific regions.

Is there anybody who has experience with BiSeq:

> readBismark(files, colData)
> metadata <- list(Sequencer = "Instrument", Year = "2013")
> rowRanges <- GRanges(seqnames = "chr1",ranges = IRanges(start = c(1,2,3), end = c(1,2,3)))
> colData <- DataFrame(group = c("cancer", "control"), row.names = c("sample_1", "sample_2"))
> totalReads <- matrix(c(rep(10L, 3), rep(5L, 3)), ncol = 2)
> methReads <- matrix(c(rep(5L, 3), rep(5L, 3)), ncol = 2)
> BSraw(metadata = metadata,rowRanges = rowRanges,colData = colData,totalReads = totalReads,methReads = methReads)
sequencing dmr methylation • 131 views
ADD COMMENTlink modified 16 days ago • written 16 days ago by ahmad.moousavi40

This does not answer your BiSeq question directly, but I've been very happy using metilene for DMR calling.

ADD REPLYlink written 16 days ago by Chris Miller19k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 854 users visited in the last hour