我想使用入藏号获取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文件 请你在这方面提供亲切的帮助 问候 Ali Zohaib答案 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)