[Bioperl-l] remoteblast xml problem

Hubert Prielinger hubert.prielinger at gmx.at
Sat Jun 3 01:36:42 UTC 2006


hi chris,
thanks but I never intended to run the remoteblast with so much, only a 
few of them, acutally I goal is to run the phiblast with regular 
expression, so that i just don't need that
file anymore.

another question for parsing the xml output....is there a xml parser 
available for blast xml output or how to start.....
I have looked up at the wikiperl and cpan Bio::SearchIO::blastxml, but 
I'm not sure how to start....sorry, I guess I'm too stupid....
is their maybe another introduction or an example.

thanks
Hubert


Chris Fields wrote:
> Yes, I see the same error you do.  But I have a similar script  
> (blastp, XML blast report, XML parsing, similar loop structure) that  
> works fine.  I'm trying to dissect the problem but I think it may be  
> something logically wrong here (something not so obvious) and not a  
> bug...
>
> What I'm trying to say is, when you send sequences using remoteblast  
> like, this you are essentially spamming the NCBI BLAST server with  
> ~1600 requests.  This script wasn't set up with that intent in mind;  
> you should really try to set up your own local blast database if  
> possible.  If you can't, try running this script in off-hours  
> (10pm-6am EST or something like that).
>
>
> Chris
>
> On Jun 2, 2006, at 7:49 PM, Hubert Prielinger wrote:
>
>   
>> hi,
>> input database: swissprot
>>         matrix: pam30
>>         count: 1
>>         gapcosts: 9 1
>>
>> I know that there are  a lot of sequences, but that doesn't matter,  
>> you can delete all of them except one, the amount of the sequences  
>> is not the problem, the script reads one line and submits  
>> it.....then the second line and so on.....I have tried it with only  
>> one sequence either and I got the same result.... the script run at  
>> that time for more than 20 minutes!!!!!! .....and that should be  
>> enough time to retrieve the results for ONE sequence, I guess
>>
>> regards
>> Hubert
>>
>>
>>
>> Chris Fields wrote:
>>     
>>> You need to add the input conditions as well (you have several  
>>> <STDIN> lines which may play a role; I would like to know what you  
>>> normally enter for those).
>>>
>>> How long did you let the script run?  I ran a quick check on your  
>>> sequences; you have almost 1600, so you have to expect that you'll  
>>> run into some problems here!  Most here (including me) would  
>>> suggest you try installing a local blast setup for something like  
>>> this.
>>>
>>> Chris
>>>
>>> On Jun 2, 2006, at 6:19 PM, Hubert Prielinger wrote:
>>>
>>>       
>>>> hi,
>>>> I have submitted the bug -> Bug 2017
>>>> with the script and input file, just start it from command line
>>>>
>>>> thank you very much
>>>> greetings
>>>>
>>>> Hubert
>>>>
>>>> Chris Fields wrote:
>>>>         
>>>>> Hubert,
>>>>>
>>>>> I have a script that's using blastxml and XML output which seems  
>>>>> to work.
>>>>> I'll try looking at it to get a better idea this weekend.
>>>>>
>>>>> Chris
>>>>>
>>>>>
>>>>>           
>>>>>> -----Original Message-----
>>>>>> From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-
>>>>>> bounces at lists.open-bio.org] On Behalf Of Hubert Prielinger
>>>>>> Sent: Friday, June 02, 2006 4:12 PM
>>>>>> To: Chris Fields; bioperl-l at bioperl.org; Chris Fields; 'Sendu  
>>>>>> Bala'
>>>>>> Subject: Re: [Bioperl-l] remoteblast xml problem
>>>>>>
>>>>>> hi,
>>>>>> sorry, but I have updated the remoteblast module and I have run  
>>>>>> several
>>>>>> attempts with the same results as before. It didn't work.
>>>>>> I didn't get any results.
>>>>>>
>>>>>> regards
>>>>>> Hubert
>>>>>>
>>>>>>
>>>>>> Chris Fields wrote:
>>>>>>
>>>>>>             
>>>>>>> Sendu, Hubert,
>>>>>>>
>>>>>>>
>>>>>>> Hubert, your code looks fine so Sendu's patch should fix the  
>>>>>>> problem
>>>>>>>
>>>>>>>               
>>>>>> (break
>>>>>>
>>>>>>             
>>>>>>> out of that infinite loop).  I applied Sendu's patch to  
>>>>>>> RemoteBlast in
>>>>>>>
>>>>>>>               
>>>>>> CVS;
>>>>>>
>>>>>>             
>>>>>>> it passed all tests in RemoteBlast.t.  Try updating from CVS  
>>>>>>> to see if
>>>>>>>
>>>>>>>               
>>>>>> it
>>>>>>
>>>>>>             
>>>>>>> works.
>>>>>>>
>>>>>>> Chris
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>               
>>>>>>>> -----Original Message-----
>>>>>>>> From: bioperl-l-bounces at lists.open-bio.org [mailto:bioperl-l-
>>>>>>>> bounces at lists.open-bio.org] On Behalf Of Sendu Bala
>>>>>>>> Sent: Friday, June 02, 2006 4:04 AM
>>>>>>>> To: bioperl-l at lists.open-bio.org
>>>>>>>> Subject: Re: [Bioperl-l] remoteblast xml problem
>>>>>>>>
>>>>>>>> Hubert Prielinger wrote:
>>>>>>>>
>>>>>>>>
>>>>>>>>                 
>>>>>>>>> hi,
>>>>>>>>> I have the following program and it worked quite well, for  
>>>>>>>>> retrieving
>>>>>>>>> remoteblast results in a textfile,
>>>>>>>>> now I have altered it to to xml, and it didn't work  
>>>>>>>>> anymore.....
>>>>>>>>> it takes all the parameter at the commandline, submits the  
>>>>>>>>> query, but
>>>>>>>>>
>>>>>>>>>                   
>>>>>> I
>>>>>>
>>>>>>             
>>>>>>>>> don't retrieve any results file anymore.....
>>>>>>>>>
>>>>>>>>> it seems that it hangs in a endless loop......
>>>>>>>>> the only output I get is:  $rc is not a ref! over and  
>>>>>>>>> over..... it
>>>>>>>>> doesn't enter the else term anymore....
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>                   
>>>>>>>> There is no problem with your code. The problem is with the  
>>>>>>>> NCBI server
>>>>>>>> and should be reported to them. You can visit the site and do  
>>>>>>>> a blast,
>>>>>>>> requesting xml format, and you will typically get one normal  
>>>>>>>> 'waiting'
>>>>>>>> message and the promise that it will be updated in x seconds,  
>>>>>>>> but
>>>>>>>> subsequent attempts to get progress information result in an  
>>>>>>>> xml error
>>>>>>>> page because the NCBI server doesn't actually send any data.
>>>>>>>>
>>>>>>>> Unfortunately the way that the bioperl code is written, it  
>>>>>>>> treats no
>>>>>>>> data as 'waiting' instead of an error. I've offered a patch  
>>>>>>>> to fix this
>>>>>>>> at this bug page:
>>>>>>>> http://bugzilla.bioperl.org/show_bug.cgi?id=2015
>>>>>>>> _______________________________________________
>>>>>>>> Bioperl-l mailing list
>>>>>>>> Bioperl-l at lists.open-bio.org
>>>>>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>>>>>>
>>>>>>>>
>>>>>>>>                 
>>>>>>> _______________________________________________
>>>>>>> Bioperl-l mailing list
>>>>>>> Bioperl-l at lists.open-bio.org
>>>>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>               
>>>>>> _______________________________________________
>>>>>> Bioperl-l mailing list
>>>>>> Bioperl-l at lists.open-bio.org
>>>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>>>>
>>>>>>             
>>>>>
>>>>>           
>>>> _______________________________________________
>>>> Bioperl-l mailing list
>>>> Bioperl-l at lists.open-bio.org
>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>>         
>>> Christopher Fields
>>> Postdoctoral Researcher
>>> Lab of Dr. Robert Switzer
>>> Dept of Biochemistry
>>> University of Illinois Urbana-Champaign
>>>
>>>
>>>
>>>
>>>       
>
> Christopher Fields
> Postdoctoral Researcher
> Lab of Dr. Robert Switzer
> Dept of Biochemistry
> University of Illinois Urbana-Champaign
>
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
>
>   




More information about the Bioperl-l mailing list