我已经在r中编写了一个非常快速的高级脚本,以便与NCBI blast API进行交互。然而,有时候,结果url需要一段时间才能加载,我的脚本会抛出一个错误,直到url准备就绪。有没有一种优雅的方式(即tryCatch选项)来处理错误,直到结果返回或在指定时间后超时? library(rvest)
## Definitive set of blast API instructions ca
我正在使用以下格式生成数百个不同文件的管道(我在字段中写入X,我不在乎): id1 X X X X X X X X X evalue1 X
id2 X X X X X X X X X evalue2 X
...
我要过滤这些文件,对于每一个ID,采取基础上,安勤(越小越好)最好的结果,但如果最好的安勤重复使用相同的ID不计该ID。 作为例子,如果输入文件是: id1 X X X X X
我需要创建一个程序,可以将一些文本文件称为fasta文件并将其转换为sequence_name,Domain_names,Domain of Domain,Domain of Domain。 所以一个FASTA文件只是一个文本文件,它看起来像这样 >MICE_8
ATTCGATCGATCGATTTCGATCGATCGATCGATCGGGATCGATCGATCGATCGATC
>MICE_59
我想要使用登录号码获取Get taxonomic层次结构。我发现这样做最简单的方法是 library('taxize')
for (year in c("AY744148.1","AY656167.1","AY656168.1")){print(paste(year))}
classification(genbank2uid(id = year), db = "ncbi")
R的给我下面的