1

我有一些.vcf文件。我从我的目录中选择了这些文件，并希望将它们转换为其他两种格式。

我在这里使用if和有点困惑else if。我想这样做：如果该文件没有文件.bgz，我想将其转换为保留原始文件的文件。[i].vcf.bgz

如果已经有.bgz文件，但没有.bgz.tbi该文件的[i]文件.bgz，那么我想将 .bgz文件转换为文件，保留从文件中获得.bgz.tbi的原始文件。.bgz.vcf

有人可以帮我完成这个循环吗？它适用于 if 条件，但不知道如何从那里开始。

  path.file<-"/mypath/for/files/" 
  all.files <- list.files("/mypath/for/files")
  all.files <- all.files[grepl(".vcf$",all.files)]

for (i in 1:length(all.files)){
 if(!exists(paste0(all.files[i],".bgz"))){
  bgzip(paste0(path.file,all.files[i]), overwrite=FALSE)
  }else{(!exists(paste0(all.files[i],".bgz",".tbi"))){
  #if(!exists(paste0(all.files[i],".bgz",".tbi"))){
  indexTabix(paste0(paste0(path.file,all.files[i]),".bgz"), format="vcf")
  }
 }
4

1 回答 1

3

试试这个（未测试）：

#get VCF files with path
all.files <- list.files("/mypath/for/files", pattern = "*.vcf$",
                        full.names = TRUE)

for (i in all.files) {
  #make output names, so we don't mess about with paste
  file_bgz <- paste0(i, ".bgz")
  file_bgz_tbi <- paste0(i, ".bgz.tbi")

  #if bgz exists don't zip else zip
  if(!exists(file_bgz))
    bgzip(i, paste0(i, ".bgz"))

  #if tbi exists don't index else tabix
  if(!exists(file_bgz_tbi))
    indexTabix(file_bgz, format = "vcf")
  }
于 2016-04-21T11:35:00.857 回答