Mercurial > repos > mahtabm > ensembl
diff variant_effect_predictor/Bio/Map/Position.pm @ 0:1f6dce3d34e0
Uploaded
author | mahtabm |
---|---|
date | Thu, 11 Apr 2013 02:01:53 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/variant_effect_predictor/Bio/Map/Position.pm Thu Apr 11 02:01:53 2013 -0400 @@ -0,0 +1,202 @@ +# $Id: Position.pm,v 1.9 2002/10/22 07:45:16 lapp Exp $ +# +# BioPerl module for Bio::Map::Position +# +# Cared for by Jason Stajich <jason@bioperl.org> +# +# Copyright Jason Stajich +# +# You may distribute this module under the same terms as perl itself + +# POD documentation - main docs before the code + +=head1 NAME + +Bio::Map::Position - A single position of a Marker in a Map + +=head1 SYNOPSIS + + use Bio::Map::Position; + my $position = new Bio::Map::Position(-map => $map, + -marker => $marker + -value => 100 + ); + +=head1 DESCRIPTION + +This object is an implementation of the PositionI interface that +handles the specific values of a position. This allows an element +(e.g. Marker) to have multiple positions within a map and still be +treated as a single entity. + +This does not directly handle the concept of a relative map in which +no known exact positions exist but markers are just ordered relative +to one another - in that case arbitrary values must be assigned for +position values. + +No units are assumed here - units are handled by context of which Map +a position is placed in or the subclass of this Position. + +=head1 FEEDBACK + +=head2 Mailing Lists + +User feedback is an integral part of the evolution of this and other +Bioperl modules. Send your comments and suggestions preferably to +the Bioperl mailing list. Your participation is much appreciated. + + bioperl-l@bioperl.org - General discussion + http://bioperl.org/MailList.shtml - About the mailing lists + +=head2 Reporting Bugs + +Report bugs to the Bioperl bug tracking system to help us keep track +of the bugs and their resolution. Bug reports can be submitted via +email or the web: + + bioperl-bugs@bioperl.org + http://bugzilla.bioperl.org/ + +=head1 AUTHOR - Jason Stajich + +Email jason@bioperl.org + +Describe contact details here + +=head1 CONTRIBUTORS + +Lincoln Stein, lstein@cshl.org +Heikki Lehvaslaiho, heikki@ebi.ac.uk +Chad Matsalla, bioinformatics1@dieselwurks.com + +=head1 APPENDIX + +The rest of the documentation details each of the object methods. +Internal methods are usually preceded with a _ + +=cut + + +# Let the code begin... + + +package Bio::Map::Position; +use vars qw(@ISA); +use strict; + +# Object preamble - inherits from Bio::Root::Root + +use Bio::Root::Root; +use Bio::Map::PositionI; + +@ISA = qw(Bio::Root::Root Bio::Map::PositionI ); + +=head2 new + + Title : new + Usage : my $obj = new Bio::Map::Position(); + Function: Builds a new Bio::Map::Position object + Returns : Bio::Map::Position + Args : -map a <Bio::Map::MapI> object + -marker a <Bio::Map::MarkerI> object + -value string or number + +=cut + +sub new { + my($class,@args) = @_; + my $self = $class->SUPER::new(@args); + + my ($map, $marker, $value) = + $self->_rearrange([qw( MAP + MARKER + VALUE + )], @args); + + $map && $self->map($map); + $marker && $self->marker($marker); + $value && $self->value($value); + + return $self; +} + +=head2 map + + Title : map + Usage : my $id = map->$map; + Function: Get/Set the map the position is in. + Returns : L<Bio::Map::MapI> + Args : [optional] new L<Bio::Map::MapI> + +=cut + +sub map { + my ($self,$map) = @_; + if( defined $map ) { + $self->throw("This is [$map], not a Bio::Map::MapI object") + unless $map->isa('Bio::Map::MapI'); + $self->{'_map'} = $map; + } + return $self->{'_map'}; +} + +=head2 marker + + Title : marker + Usage : my $id = marker->$marker; + Function: Get/Set the marker the position is in. + Returns : L<Bio::Map::MarkerI> + Args : [optional] new L<Bio::Map::MarkerI> + +=cut + +sub marker { + my ($self,$marker) = @_; + if( defined $marker ) { + $self->thow("This is [$marker], not a Bio::Map::MarkerI object") + unless $marker->isa('Bio::Map::MarkerI'); + $self->{'_marker'} = $marker; + } + return $self->{'_marker'}; +} + +=head2 value + + Title : value + Usage : my $pos = $position->value; + Function: Get/Set the value for this postion + Returns : scalar, value + Args : [optional] new value to set + +=cut + +sub value { + my ($self,$value) = @_; + if( defined $value ) { + $self->{'_value'} = $value; + } + return $self->{'_value'}; +} + +=head2 numeric + + Title : numeric + Usage : my $num = $position->numeric; + Function: Read-only method that is guarantied to return a numeric + representation for this position. + Returns : numeric (int or real) + Args : none + +=cut + +sub numeric { + my ($self) = @_; + my $num = $self->{'_value'} || 0; + + # expand this to cover scientific notation, too! + $self->throw("This value [$num] is not numeric!") + unless $num && $num =~ /^[+-]?[\d.]+$/; + return $num; +} + +1;