Ola pessoal, eu estou tentando procurar e extrair algumas sequências de DNA pelo R<br>Em um livro que estou lendo fui guiado até o pacote seqinr, ele tem uma função chamada query() pra fazer isso.<br><br>A principio tudo parece bem simples e intuitivo.<br>
Por exemplo, eu tenho interesse na sequência 16s rRNA de varias especies de sapinhos.<br><br>Ai se procuro pra uma especie tenho:<br><br>> library(seqinr)<br>> s <- choosebank("genbank")<br>> query("SP","SP=Rhinella rubescens AND K=16S RIBOSOMAL RNA",s$socket)<br>
> SP$nelem<br>[1] 1<br>> getName(SP)<br>[1] "GU907196.RR1"<br><br>Blz eu achei 1 sequência especifica pro que eu quero, e o nomezinho pra acessar ela.<br><br>Procuro outra<br>> query("SP","SP=Leptodactylus syphax AND K=16S RIBOSOMAL RNA",s$socket)<br>
> print(SP$nelem)<br>[1] 2<br>> print(getName(SP))<br>[1] "JF789923" "JF789924"<br><br>Blz, 2 sequencias,ai eu fui ao site do genbank, que é a base da dados da uma olhada e estao la blz<br><br>So que muitas especies que procuro  o R para.<br>
Por exemplo:<br>> query("SP","SP=Leptodactylus marmoratus AND K=16S RIBOSOMAL RNA",s$socket)<br><br>O R fica parado e tenho que fechar ele.<br>Eu não entendo o que esta dando errado, principalmente pq algumas especies funcionam e outras não.<br>
Essa especie em especifico num era pra achar nada, não achei registro dela procurando pelo site do genbank.<br>Mas a hora que vou usar o R num vai.<br>Outro exemplo é "Leptodactylus syphax" que devia acha 2 sequencias mas ele para também.<br>
<br>E a maior parte das especies num vai.<br>Eu uso windows 7 como SO. Estou com o R 2.15.0, internet residencial.<br><br>Alguem tem experiencia com isso? Sabe dizer o que esto fazendo de errado?<br>E alguem conhece alguma outra forma de fazer o mesmo tipo de busca usando o R? O task view indica o Biocondutor que tem o pacote SRAdb que deve fazer a mesma coisa, mas não consegui usar. <br>
<br>Eu vi que o livro que é de 2011 a função mudou a ordem dos argumentos, e teve um update no pacote, mas não sei se estou fazendo algo errado, se tenho que fazer mais alguma coisa. Os manuais que vi não tem um faq com problemas comum e soluções.<br>
<br>Minha lista de especies caso alguem queria ver é:<br>lista.spp<-c("Rhinella rubescens", "Rhinella fernandezae", "Elosia rustica", <br>"Crossodactylus gaudichaudii", "Leptodactylus pustulatus", "Leptodactylus syphax", <br>
"Dendropsophus cachimbo", "Leptodactylus marmoratus", "Physalaemus soaresi", <br>"Hylodes nasus", "Scinax fuscomarginatus", "Leptodactylus petersii", <br>"Chiasmocleis capixaba", "Ceratophrys aurita", "Pleurodema diplolister", <br>
"Cystignatus gigas", "Scinax trilineatus", "Elosis nasus", "Dryadophis bifossatus"<br>)<br><br>Se alguem puder dar uma luz, não sei nem por onde começar a procurar pra entender o que ta acontecendo. <br>
Vou quebrando a cabeça por enquanto, Boa noite a todos.<br><br clear="all"><br>-- <br><div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://augustoribas.heliohost.org" target="_blank">http://augustoribas.heliohost.org</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div><br>