如何在R中编写csv文件以获取以下数据?

时间:2017-06-22 08:36:55

标签: r parsing export-to-csv taxonomy blast

我想使用入藏号获取Get taxonomic层次结构。我发现最简单的方法是

library('taxize')
for (year in c("AY744148.1","AY656167.1","AY656168.1")){print(paste( year))}
classification(genbank2uid(id = year), db = "ncbi")

R正在给我以下输出,  但是我无法从这个

写一个csv文件

enter image description here

我想在csv文件中得到一个结果,如下图所示   需要输出  enter image description here

请你在这方面提供亲切的帮助 问候 Ali Zohaib

1 个答案:

答案 0 :(得分:0)

也许你是在追求这样的事情?

library('taxize')

IDs <- c("AY744148.1","AY656167.1","AY656168.1") # your sequence of ids to iterate over
output <- vector("list", length(years)) # prepare a list to write to

# for each ID, run the function and store the result into a list
for (i in 1:length(years)) {
  # output is a list of length one, so we need to subset it using "[["
  output[[i]] <- classification(genbank2uid(id = IDs[i]), db = "ncbi")[[1]]
  }

# take list elements of output and merge them row-wise (rbind == row bind)
output <- do.call(rbind, output)

                   name    rank    id
1        other sequences no rank 28384
2   artificial sequences no rank 81077
3                vectors no rank 29278
4  Cloning vector pBR322 species 47470
5        other sequences no rank 28384
6   artificial sequences no rank 81077
7                vectors no rank 29278
8  Cloning vector pGEM-3 species 90108
9        other sequences no rank 28384
10  artificial sequences no rank 81077
11               vectors no rank 29278
12 Cloning vector pGEM-3 species 90108

您可以使用

编写结果
write.table(output, file = "endfile.txt", row.names = FALSE, quote = FALSE)