Description Usage Arguments Value Author(s) Examples
Plotting the results of DMC identifation stored in
a BSDMC-class
object
1 2 3 4 |
object |
A |
region |
An integer vector of length two specifying which subset of the object to be plotted |
nSplit |
A integer value specifying the number of subsets must be done for plotting the results of DMC identification |
parList |
A list specifying plots parameters, see |
Plot
Farhad Shokoohi <shokoohi@icloud.com>
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | set.seed(1980)
nr <- 1000
nc <- 4
metht <- matrix(as.integer(runif(nr * nc, 0, 100)), nr)
methc <- matrix(rbinom(n = nr * nc, c(metht), prob = runif(nr * nc)), nr, nc)
methl <- methc / metht
r1 <- GRanges(rep("chr1", nr), IRanges(1:nr, width = 1), strand = "*")
names(r1) <- 1:nr
cd1 <- DataFrame(
Group = rep(c("G1", "G2"), each = nc / 2),
row.names = LETTERS[1:nc]
)
OBJ1 <- cBSDMC(
rowRanges = r1, methReads = methc, totalReads = metht,
methLevels = methl, colData = cd1
)
OBJ2 <- findDMCFB(OBJ1,
bwa = 10, bwb = 10, nBurn = 50, nMC = 50, nThin = 1,
alpha = 0.05, nCores = 2, pSize = 500, sfiles = FALSE
)
plotDMCFB(OBJ2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.