【发布时间】:2017-01-09 20:07:34
【问题描述】:
以下是我的代码。我正在尝试获取以.idat 结尾的所有文件(~20000)的列表,并使用函数illuminaio::readIDAT 读取每个文件。
library(illuminaio)
library(parallel)
library(data.table)
# number of cores to use
ncores = 8
# this gets all the files with .idat extension ~20000 files
files <- list.files(path = './',
pattern = "*.idat",
full.names = TRUE)
# function to read the idat file and create a data.table of filename, and two more columns
# write out as csv using fwrite
get.chiptype <- function(x)
{
idat <- readIDAT(x)
res <- data.table(filename = x, nSNPs = nrow(idat$Quants), Chip = idat$ChipType)
fwrite(res, file.path = 'output.csv', append = TRUE)
}
# using mclapply call the function get.chiptype on all 20000 files.
# use 8 cores at a time
mclapply(files, FUN = function(x) get.chiptype(x), mc.cores = ncores)
读取和写入有关 1200 个文件的信息后,我收到以下消息:
Warning message:
In mclapply(files, FUN = function(x) get.chiptype(x), mc.cores = ncores) :
all scheduled cores encountered errors in user code
我该如何解决?
【问题讨论】:
-
什么是destdir和destfile
-
它们只是将写入 data.table 的目录和文件名。我会删除它。
-
你仍然得到错误?
-
这可能不是问题,但我会警惕从并行进程追加到单个文件。我不是专家,但这似乎是一个麻烦的秘诀。你知道他们是否以某种方式锁定了文件,所以他们一次只能写一个?
-
这不是问题。这些只是我的源文件和目标文件的路径。这就是并行处理的问题。
标签: r parallel-processing mclapply