我在一个文本文件中有很多序列。我使用“read.fasta”函数导入这些序列。我使用“for loop”为每个序列创建核苷酸频率表,并使用“write.table”进行输出。但是它为每个序列创建一个文件(许多输出文件和每个文件都有一个序列表)。我搜索一个命令来创建一个包含所有表的文件。
注意:“mydata.txt”是一个包含许多fasta格式序列的文件
file=read.fasta(file="mydata")
for (i in seq_along(file))
{
NucleotidesFrequency=table(file[[i]])
print(NucleotidesFrequency)
write.table(NucleotidesFrequency, paste(i, (".txt"), sep=""),sep="\t")
}