2017-06-22 74 views
-1

我想要使用登录号码获取Get taxonomic层次结构。我发现这样做最简单的方法是如何在R中编写csv文件以获取以下数据?

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

R的给我下面的输出, 但是我无法从这个

enter image description here

写一个CSV文件我想在csv文件的结果如下图所示 需要输出 enter image description here

您需要这方面的帮助 问候 阿里Zohaib

+1

您将需要从申请家庭为'()'循环或者是'sapply'或一些其他功能的结果。 –

+0

谢谢先生,for()函数解决了输入数据的问题。现在我正面临如何从输出中编写csv文件。我已更新我的问题 –

回答

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) 
+0

谢谢,先生,它的工作 –

相关问题