一个循环函数 R 来执行任务

Question

I have file stored in directory named file1.csv, file2.csv, file3.csv....file100.csv. I have to run series of script with changing file name from 1,2,3...100. script is like

a1 <- data2haplohh(hap_file="file1.csv", polarize_vcf=FALSE)

followed by

a_1<-subset(a1, min_maf=0.05)

how can we run loop over these scipts in R. I am trying but getting errors.

Answer 1

A way to do in a lapply loop -

filenames <- list.files(pattern = 'file\\d+\\.csv', full.names = TRUE)

result <- lapply(filenames, function(x) {
  a1 <- data2haplohh(hap_file=x, polarize_vcf=FALSE)
  subset(a1, min_maf==0.05)  
})

result should be a list of length 100 if you have 100 files.

Answer 2

We can try like this

list2env(
    setNames(
        lapply(
            seq(100),
            function(i) {
                subset(
                    data2haplohh(hap_file = sprintf("file%s.csv", i), polarize_vcf = FALSE),
                    min_maf = 0.05
                )
            }
        ), 
        paste0("a_", seq(100))
    ),
    envir = .GlobalEnv
)