# NOT RUN {
data(methylKit)
methylRawDB.obj=methRead( system.file("extdata","test1.txt.bgz",package="methylKit"),
sample.id="test1", assembly="hg18",
dbtype = "tabix",dbdir = "methylDB")
methylBaseDB.obj=unite(methylRawList.obj,save.db=TRUE,dbdir="methylDB")
# selects first hundred rows, returns a methylRaw object
subset1=select(methylRawList.obj[[1]],1:100)
subset1=select(methylRawDB.obj,1:100)
# selects first hundred rows, returns a methylBase object
subset2=select(methylBase.obj,1:100)
subset2=select(methylBaseDB.obj,1:100)
# selects first hundred rows, returns a methylDiff object
subset3=select(methylDiff.obj,1:100)
# remove Database again
rm(methylBaseDB.obj)
rm(methylRawDB.obj)
unlink("methylDB",recursive=TRUE)
# }
Run the code above in your browser using DataLab