[Bioperl-l] Bio ::seqIO ::tigr
Hilmar Lapp
hlapp at gmx.net
Wed Jan 28 11:55:15 EST 2004
I suspect you have an old version of bioperl-db, or a version mix-up.
You need to download and install the latest revision from CVS for
bioperl-db.
Note that if the root of the problem is with the pir parser then
load_seqdatabase.pl will not cure it, as it just uses any Bio::SeqIO
compliant parser to provide the input sequences. If the parser is
broken then there won't be input ... It just saves you the round-trip
(and possible errors associated with it) of going through swissprot
format.
-hilmar
On Wednesday, January 28, 2004, at 02:07 AM, matthieu CONTE wrote:
> Ok , I try directly with "load_seqdatabase.pl" but there is another
> problem.....
>
> [conte at bearn scripts]$ perl load_seqdatabase.pl -dbuser biosql
> -dbpass biosql -format tigr tigr
> /home/conte/pipeline_orthologues/data/orysa_tigr.txt
>
> Can't locate object method "get_BioDatabaseAdaptor" via package
> "Bio::DB::BioSQL::DBAdaptor" at load_seqdatabase.pl line 84.
>
> Indeed this method does not exist in Bio::DB::BioSQL::DBAdaptor....
>
>
>
>
> Matthieu CONTE
> M. Sc. in Bioinformatics from SIB
>
> 00 33 06.68.90.28.70
> m_conte at hotmail.com
>
>
>
>
>
>> From: Hilmar Lapp <hlapp at gmx.net>
>> To: "matthieu CONTE" <m_conte at hotmail.com>
>> CC: bioperl-l at bioperl.org
>> Subject: Re: [Bioperl-l] Bio ::seqIO ::tigr Date: Tue, 27 Jan 2004
>> 09:31:39 -0800
>>
>> A question aside: why do you want to convert to swissprot in order to
>> load into biosql? (load_seqdatabase.pl can use any SeqIO reader.)
>>
>> -hilmar
>>
>> On Tuesday, January 27, 2004, at 02:50 AM, matthieu CONTE wrote:
>>
>>> I currently trying to use the Bio ::seqIO ::tigr module.
>>> My objective is to download the whole rice genome form Tigr ( adress
>>> below)and to integrate it in my BioSQL DB.
>>> For this I am trying to convert the tigr format in swiss format with
>>> the script below
>>>
>>>
>>> use Bio::SeqIO;
>>>
>>> my $in = Bio::SeqIO->new(-file
>>> =>'</home/conte/pipeline_orthologues/data/orysa_tigr.txt', -format
>>> =>'tigr');
>>>
>>> my $out = Bio::SeqIO->new(-file =>
>>> '>/home/conte/pipeline_orthologues/data/orysa_swiss.txt' ,
>>> -format=>'swiss');
>>>
>>> print $out $_ while <$in>;
>>>
>>> I obtain:
>>>
>>> ------------ EXCEPTION -------------
>>> MSG: [19]Required <AUTHOR_LIST> missing
>>> STACK Bio::SeqIO::tigr::throw
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:1338
>>> STACK Bio::SeqIO::tigr::_process_header
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:700
>>> STACK Bio::SeqIO::tigr::_process_assembly
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:535
>>> STACK Bio::SeqIO::tigr::_process_tigr
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:453
>>> STACK Bio::SeqIO::tigr::_process
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:420
>>> STACK Bio::SeqIO::tigr::_initialize
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO/tigr.pm:90
>>> STACK Bio::SeqIO::new
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO.pm:358
>>> STACK Bio::SeqIO::new
>>> /usr/local/ActivePerl-5.8/lib/site_perl/5.8.0/Bio/SeqIO.pm:378
>>> STACK toplevel get_bioseq_tigr.pl:8
>>>
>>> Could you please tell me if there is a problem with the parser or
>>> with the input data format of Tigr?
>>>
>>> Thanks in advance
>>>
>>>
>>>
>>>
>>> Matthieu CONTE
>>> m_conte at hotmail.com
>>>
>>> _________________________________________________________________
>>> MSN Messenger : discutez en direct avec vos amis !
>>> http://www.msn.fr/msger/default.asp
>>>
>>> _______________________________________________
>>> Bioperl-l mailing list
>>> Bioperl-l at portal.open-bio.org
>>> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>>>
>>>
>> --
>> -------------------------------------------------------------
>> Hilmar Lapp email: lapp at gnf.org
>> GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
>> -------------------------------------------------------------
>>
>>
>
> _________________________________________________________________
> MSN Messenger : discutez en direct avec vos amis !
> http://www.msn.fr/msger/default.asp
>
>
--
-------------------------------------------------------------
Hilmar Lapp email: lapp at gnf.org
GNF, San Diego, Ca. 92121 phone: +1-858-812-1757
-------------------------------------------------------------
More information about the Bioperl-l
mailing list