# load methylation data object
data(methylationDataList)
# save the two datasets
saveBismark(methylationDataList[["WT"]],
"chr3test_a_thaliana_wt.CX_report")
saveBismark(methylationDataList[["met1-3"]],
"chr3test_a_thaliana_met13.CX_report")
# reload the two datasets and pool them
filenames <- c("chr3test_a_thaliana_wt.CX_report",
"chr3test_a_thaliana_met13.CX_report")
methylationDataPool <- readBismarkPool(filenames)
Run the code above in your browser using DataLab