如何解决如何从不同的文件中提取相同的数据
为了从 qiime2 工件中提取频率表并将其写入 tsv 输出，您必须：
rel_freq_table <- read_qza("05.TaxonomyClassification/rel_freq/1_rel_freq_table.qza")
rel_freq_table <- data.frame(rel_freq_table$data)
rel_freq_table <- setDT(rel_freq_table,keep.rownames = TRUE )
colnames(rel_freq_table)[1]<- c( "Feature.ID" )
write.table(rel_freq_table,file = "1_rel_freq_table.tsv",row.names = FALSE,dec = ",")
现在，我有七个 qiime2 工件
rel_freqs <- c("05.TaxonomyClassification/rel_freq/1_rel_freq_table.qza","05.TaxonomyClassification/rel_freq/2_rel_freq_table.qza","05.TaxonomyClassification/rel_freq/3_rel_freq_table.qza","05.TaxonomyClassification/rel_freq/4_rel_freq_table.qza","05.TaxonomyClassification/rel_freq/5_rel_freq_table.qza","05.TaxonomyClassification/rel_freq/6_rel_freq_table.qza","05.TaxonomyClassification/rel_freq/7_rel_freq_table.qza")
我想自动化这个过程，所以我只需要运行一次命令。我尝试过循环和 lapply 但没有运气，有人知道吗？
解决方法
您可以编写一个函数并使用 lapply 将其应用到每个文件：
read_qiime2_write_tsv <- function(filename) {
rel_freq_table <- read_qza(filename)
rel_freq_table <- data.frame(rel_freq_table$data)
rel_freq_table <- setDT(rel_freq_table,keep.rownames = TRUE )
colnames(rel_freq_table)[1]<- c( "Feature.ID" )
write.table(rel_freq_table,file = paste0(tools::file_path_sans_ext(basename(filename)),'.tsv'),row.names = FALSE,dec = ",")
}
lapply(rel_freqs,read_qiime2_write_tsv)
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。