inst/scripts/create.data.R

In methylKit: DNA methylation analysis from high-throughput bisulfite sequencing results

library(methylKit)
file.list=list( system.file("extdata", "test1.myCpG.txt", package = "methylKit"),
                system.file("extdata", "test2.myCpG.txt", package = "methylKit"),
                system.file("extdata", "control1.myCpG.txt", package = "methylKit"),
                system.file("extdata", "control2.myCpG.txt", package = "methylKit") )

methylRawList.obj=methRead( file.list,
                sample.id=list("test1","test2","ctrl1","ctrl2"),assembly="hg18",treatment=c(1,1,0,0))


methylBase.obj=unite(methylRawList.obj)

methylDiff.obj=calculateDiffMeth(methylBase.obj)

save(methylRawList.obj,methylBase.obj,methylDiff.obj,file="methylKit.RData")