R 如何从不同的文件中提取相同的数据
为了从qiime2工件中提取频率表并将其写入tsv输出,您必须:R 如何从不同的文件中提取相同的数据,r,qiime,R,Qiime,为了从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
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 = ",")
rel\u freq\u table您可以编写一个函数,并使用lappy
将其应用于每个文件:
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)
read\u qiime2\u write\u tsv
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)