library(RnBeads.hg19)
data(small.example.object)
logger.start(fname=NA)
dm1 <- rnb.execute.computeDiffMeth(rnb.set.example,pheno.cols=c("Sample_Group"),region.types=c("genes","tiling"))
dm2 <- rnb.execute.computeDiffMeth(rnb.set.example,pheno.cols=c("Sample_Group","Treatment"),region.types=c("promoters"))
dm.join1 <- join.diffMeth(dm1,dm2)
#the following joint object is invalid, because some region type - comparison combinations are missing
is.valid(dm.join1)
dm3 <- rnb.execute.computeDiffMeth(rnb.set.example,pheno.cols=c("Treatment"),region.types=c("genes","tiling"))
dm.join2 <- join.diffMeth(dm.join1,dm3)
#after joining the missing information, the new object is valid
is.valid(dm.join2)
Run the code above in your browser using DataLab