# NOT RUN {
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") )
methylRawListDB.obj <- methRead(file.list,
sample.id = list("test1","test2","ctrl1","ctrl2"),
assembly = "hg18",treatment = c(1,1,0,0),
dbtype = "tabix",dbdir = "methylDB")
#applying functions designed for methylRawDB on methylRawListDB object
lapply(methylRawListDB.obj,"getAssembly")
# remove Database again
rm(methylRawListDB.obj)
unlink("methylDB",recursive=TRUE)
# }
Run the code above in your browser using DataLab