annotate variant_effect_predictor/Bio/Map/Position.pm @ 0:21066c0abaf5 draft

Uploaded
author willmclaren
date Fri, 03 Aug 2012 10:04:48 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
1 # $Id: Position.pm,v 1.9 2002/10/22 07:45:16 lapp Exp $
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
2 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
3 # BioPerl module for Bio::Map::Position
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
4 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
5 # Cared for by Jason Stajich <jason@bioperl.org>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
6 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
7 # Copyright Jason Stajich
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
8 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
9 # You may distribute this module under the same terms as perl itself
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
10
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
11 # POD documentation - main docs before the code
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
12
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
13 =head1 NAME
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
14
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
15 Bio::Map::Position - A single position of a Marker in a Map
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
16
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
17 =head1 SYNOPSIS
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
18
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
19 use Bio::Map::Position;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
20 my $position = new Bio::Map::Position(-map => $map,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
21 -marker => $marker
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
22 -value => 100
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
23 );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
24
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
25 =head1 DESCRIPTION
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
26
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
27 This object is an implementation of the PositionI interface that
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
28 handles the specific values of a position. This allows an element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
29 (e.g. Marker) to have multiple positions within a map and still be
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
30 treated as a single entity.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
31
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
32 This does not directly handle the concept of a relative map in which
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
33 no known exact positions exist but markers are just ordered relative
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
34 to one another - in that case arbitrary values must be assigned for
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
35 position values.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
36
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
37 No units are assumed here - units are handled by context of which Map
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
38 a position is placed in or the subclass of this Position.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
39
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
40 =head1 FEEDBACK
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
41
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
42 =head2 Mailing Lists
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
43
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
44 User feedback is an integral part of the evolution of this and other
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
45 Bioperl modules. Send your comments and suggestions preferably to
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
46 the Bioperl mailing list. Your participation is much appreciated.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
47
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
48 bioperl-l@bioperl.org - General discussion
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
49 http://bioperl.org/MailList.shtml - About the mailing lists
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
50
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
51 =head2 Reporting Bugs
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
52
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
53 Report bugs to the Bioperl bug tracking system to help us keep track
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
54 of the bugs and their resolution. Bug reports can be submitted via
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
55 email or the web:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
56
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
57 bioperl-bugs@bioperl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
58 http://bugzilla.bioperl.org/
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
59
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
60 =head1 AUTHOR - Jason Stajich
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
61
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
62 Email jason@bioperl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
63
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
64 Describe contact details here
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
65
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
66 =head1 CONTRIBUTORS
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
67
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
68 Lincoln Stein, lstein@cshl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
69 Heikki Lehvaslaiho, heikki@ebi.ac.uk
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
70 Chad Matsalla, bioinformatics1@dieselwurks.com
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
71
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
72 =head1 APPENDIX
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
73
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
74 The rest of the documentation details each of the object methods.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
75 Internal methods are usually preceded with a _
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
76
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
77 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
78
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
79
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
80 # Let the code begin...
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
81
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
82
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
83 package Bio::Map::Position;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
84 use vars qw(@ISA);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
85 use strict;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
86
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
87 # Object preamble - inherits from Bio::Root::Root
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
88
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
89 use Bio::Root::Root;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
90 use Bio::Map::PositionI;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
91
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
92 @ISA = qw(Bio::Root::Root Bio::Map::PositionI );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
93
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
94 =head2 new
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
95
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
96 Title : new
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
97 Usage : my $obj = new Bio::Map::Position();
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
98 Function: Builds a new Bio::Map::Position object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
99 Returns : Bio::Map::Position
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
100 Args : -map a <Bio::Map::MapI> object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
101 -marker a <Bio::Map::MarkerI> object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
102 -value string or number
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
103
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
104 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
105
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
106 sub new {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
107 my($class,@args) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
108 my $self = $class->SUPER::new(@args);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
109
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
110 my ($map, $marker, $value) =
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
111 $self->_rearrange([qw( MAP
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
112 MARKER
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
113 VALUE
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
114 )], @args);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
115
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
116 $map && $self->map($map);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
117 $marker && $self->marker($marker);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
118 $value && $self->value($value);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
119
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
120 return $self;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
121 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
122
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
123 =head2 map
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
124
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
125 Title : map
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
126 Usage : my $id = map->$map;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
127 Function: Get/Set the map the position is in.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
128 Returns : L<Bio::Map::MapI>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
129 Args : [optional] new L<Bio::Map::MapI>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
130
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
131 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
132
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
133 sub map {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
134 my ($self,$map) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
135 if( defined $map ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
136 $self->throw("This is [$map], not a Bio::Map::MapI object")
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
137 unless $map->isa('Bio::Map::MapI');
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
138 $self->{'_map'} = $map;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
139 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
140 return $self->{'_map'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
141 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
142
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
143 =head2 marker
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
144
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
145 Title : marker
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
146 Usage : my $id = marker->$marker;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
147 Function: Get/Set the marker the position is in.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
148 Returns : L<Bio::Map::MarkerI>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
149 Args : [optional] new L<Bio::Map::MarkerI>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
150
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
151 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
152
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
153 sub marker {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
154 my ($self,$marker) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
155 if( defined $marker ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
156 $self->thow("This is [$marker], not a Bio::Map::MarkerI object")
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
157 unless $marker->isa('Bio::Map::MarkerI');
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
158 $self->{'_marker'} = $marker;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
159 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
160 return $self->{'_marker'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
161 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
162
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
163 =head2 value
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
164
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
165 Title : value
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
166 Usage : my $pos = $position->value;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
167 Function: Get/Set the value for this postion
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
168 Returns : scalar, value
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
169 Args : [optional] new value to set
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
170
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
171 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
172
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
173 sub value {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
174 my ($self,$value) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
175 if( defined $value ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
176 $self->{'_value'} = $value;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
177 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
178 return $self->{'_value'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
179 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
180
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
181 =head2 numeric
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
182
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
183 Title : numeric
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
184 Usage : my $num = $position->numeric;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
185 Function: Read-only method that is guarantied to return a numeric
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
186 representation for this position.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
187 Returns : numeric (int or real)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
188 Args : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
189
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
190 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
191
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
192 sub numeric {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
193 my ($self) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
194 my $num = $self->{'_value'} || 0;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
195
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
196 # expand this to cover scientific notation, too!
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
197 $self->throw("This value [$num] is not numeric!")
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
198 unless $num && $num =~ /^[+-]?[\d.]+$/;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
199 return $num;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
200 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
201
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
202 1;