read files in parallel

readMethylRawList.para.R

library(foreach)
library(methylKit)
library(doMC)
registerDoMC(2)

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") )
ids=c("test1","test2","ctrl1","ctrl2")


# read methylRaw objects to a list in parallel
x <- foreach(i=1:length(file.list)) %dopar% read(file.list[[i]],sample.id=ids[i],
                                              assembly="hg18")

# make a new methylRawList object from the list x
myMethRawList=new("methylRawList",x,treatment=c(1,1,0,0))