[Bioperl-l] RE: Bio::Tools::EPCR.pm patch request
Jason Stajich
jason at cgt.duhs.duke.edu
Tue Mar 16 11:39:39 EST 2004
oops - my fault - I missed that accidently - had trouble applying the
patch so was doing it by hand. will do.
-jason
On Mon, 15 Mar 2004, Cook, Malcolm wrote:
> Jason,
>
> Wow! You the man! That was fast. I see it in the cvs change log
> already...
>
> But...
>
> Would you mind reinstating the tiny additional change I made to only
> produce the note tag value if $rest is present. I'm using it in an
> applicatoin where there are no addition columns in my e-PCR input file,
> and am getting an empty note attribute on each GFF line....
>
> There may be a better way of coding it than:
>
> $markerfeature->add_tag_value('Note',$rest) if $rest;
>
> Thanks,
>
> Malcolm
>
> >-----Original Message-----
> >From: Jason Stajich [mailto:jason at cgt.duhs.duke.edu]
> >Sent: Monday, March 15, 2004 3:29 PM
> >To: Cook, Malcolm
> >Cc: Bioperl
> >Subject: Re: Bio::Tools::EPCR.pm patch request
> >
> >
> >done!
> >
> >I added the ability to the update the primary,source,groupclass fields
> >using get/set methods as well. I also capitalized the 'Note'
> >field to be
> >in line with the way GFF3 looks.
> >
> >-j
> >On Mon, 15 Mar 2004, Cook, Malcolm wrote:
> >
> >> Jason,
> >>
> >> I'd like it if you would accept and apply the following diff
> >to the very
> >> useful Bio::Tools::EPCR module. I hope you agree that it
> >provides more
> >> flexibility in defining the structure of the features that
> >are produced
> >> (and, ultimately for me, the format of the GFF which gets
> >written before
> >> being loaded into Bio::DB::GFF database).
> >>
> >> I have made the behaviour completely backwards compatible.
> >>
> >> Sincerely,
> >>
> >> Malcolm Cook - mec at stowers-institute.org - 816-926-4449
> >> Database Applications Manager - Bioinformatics
> >> Stowers Institute for Medical Research - Kansas City, MO USA
> >>
> >> Index: EPCR.pm
> >> ===================================================================
> >> RCS file: /home/repository/bioperl/bioperl-live/Bio/Tools/EPCR.pm,v
> >> retrieving revision 1.8
> >> diff -c -r1.8 EPCR.pm
> >> *** EPCR.pm 2002/12/01 00:05:21 1.8
> >> --- EPCR.pm 2004/03/15 20:27:04
> >> ***************
> >> *** 92,110 ****
> >> =head2 new
> >>
> >> Title : new
> >> ! Usage : my $epcr = new Bio::Tools::EPCR(-file => $file);
> >> ! Function: Initializes a new EPCR parser
> >> Returns : Bio::Tools::EPCR
> >> Args : -fh => filehandle
> >> OR
> >> -file => filename
> >>
> >> =cut
> >>
> >> sub new {
> >> my($class, at args) = @_;
> >>
> >> my $self = $class->SUPER::new(@args);
> >> $self->_initialize_io(@args);
> >>
> >> return $self;
> >> --- 92,120 ----
> >> =head2 new
> >>
> >> Title : new
> >> ! Usage : my $epcr = new Bio::Tools::EPCR(-file => $file,
> >-primary =>
> >> $fprimary, -source => $fsource, -groupclass => $fgroupclass);
> >> ! Function: Initializes a new EPCR parser.
> >> Returns : Bio::Tools::EPCR
> >> Args : -fh => filehandle
> >> OR
> >> -file => filename
> >>
> >> + -primary => a string to be used as the common value for
> >> each features '-primary' tag. Defaults to 'sts'. (This in
> >turn maps to
> >> the GFF 'type' tag (aka 'method')).
> >> +
> >> + -source => a string to be used as the common
> >value for each
> >> features '-source' tag. Defaults to 'e-PCR'. (This in turn
> >maps to the
> >> GFF 'source' tag)
> >> +
> >> + -groupclass => a string to be used as the name
> >of the tag
> >> which will hold the sts marker namefirst attribute.
> >Defaults to 'name'.
> >> +
> >> =cut
> >>
> >> sub new {
> >> my($class, at args) = @_;
> >>
> >> my $self = $class->SUPER::new(@args);
> >> + my ($primary, $source, $groupclass) =
> >$self->_rearrange([qw(PRIMARY
> >> SOURCE GROUPCLASS)], at args);
> >> + $self->{'_primary'} = $primary || 'sts';
> >> + $self->{'_source'} = $source || 'e-PCR';
> >> + $self->{'_groupclass'} = $groupclass || 'name';
> >> $self->_initialize_io(@args);
> >>
> >> return $self;
> >> ***************
> >> *** 139,152 ****
> >> my $markerfeature = new Bio::SeqFeature::Generic (
> >'-start' =>
> >> $start,
> >> '-end' =>
> >> $end,
> >> '-strand' =>
> >> $strand,
> >> ! '-source' =>
> >> 'e-PCR',
> >> ! '-primary' =>
> >> 'sts',
> >> '-seq_id' =>
> >> $seqname,
> >> '-tag' => {
> >> ! 'name'=>
> >> $mkrname,
> >> ! 'note'=>
> >> $rest,
> >> });
> >> return $markerfeature;
> >> }
> >>
> >> 1;
> >> --- 149,163 ----
> >> my $markerfeature = new Bio::SeqFeature::Generic (
> >'-start' =>
> >> $start,
> >> '-end' =>
> >> $end,
> >> '-strand' =>
> >> $strand,
> >> ! '-source' =>
> >> $self->{'_source'},
> >> ! '-primary' =>
> >> $self->{'_primary'},
> >> '-seq_id' =>
> >> $seqname,
> >> '-tag' => {
> >> !
> >> $self->{'_groupclass'}=> $mkrname,
> >> });
> >> + $markerfeature->add_tag_value('note',$rest) if $rest;
> >> return $markerfeature;
> >> +
> >> }
> >>
> >> 1;
> >>
> >
> >--
> >Jason Stajich
> >Duke University
> >jason at cgt.mc.duke.edu
> >
>
--
Jason Stajich
Duke University
jason at cgt.mc.duke.edu
More information about the Bioperl-l
mailing list