Существует функция для получения списка файлов в данном каталоге.
Let sampleList <- c(12345, 4711, 1337)
. Предполагая, что ваша структура похожа на
O:/12345_vcf/secret1.vcf
O:/4711_vcf/foo.vcf
O:/1337_vcf/bar.vcf
, и вы не знаете названия ваших файлов, но в каждом каталоге есть только один .vcf.
for (i in sampleList){
directory <- paste0("O:/", i, "_vcf")
filename <- list.files(directory, ".vcf") # if there are more than one vcf's, there are better ways
nam <- readVcf(paste0(directory, filename, sep="/"))
}
Я не знаю Я не знаю значения hg19
, поэтому я проигнорировал его. Пожалуйста, используйте его, если необходимо.