Mercurial > repos > mahtabm > ensemb_rep_gvl
diff variant_effect_predictor/Bio/Tree/AlleleNode.pm @ 0:2bc9b66ada89 draft default tip
Uploaded
author | mahtabm |
---|---|
date | Thu, 11 Apr 2013 06:29:17 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/variant_effect_predictor/Bio/Tree/AlleleNode.pm Thu Apr 11 06:29:17 2013 -0400 @@ -0,0 +1,186 @@ +# $Id: AlleleNode.pm,v 1.4 2002/10/22 07:45:24 lapp Exp $ +# +# BioPerl module for Bio::Tree::AlleleNode +# +# 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::Tree::AlleleNode - DESCRIPTION of Object + +=head1 SYNOPSIS + +Give standard usage here + +=head1 DESCRIPTION + +Describe the object here + +=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 + +Additional contributors names and emails here + +=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::Tree::AlleleNode; +use vars qw(@ISA); +use strict; + +use Bio::Tree::Node; + +@ISA = qw(Bio::Tree::Node ); + +=head2 new + + Title : new + Usage : my $obj = new Bio::Tree::AlleleNode(); + Function: Builds a new Bio::Tree::AlleleNode object + Returns : Bio::Tree::AlleleNode + Args : + + +=cut + +sub new { + my($class,@args) = @_; + + my $self = $class->SUPER::new(@args); + + my ($alleles) = $self->_rearrange([qw(ALLELES)], @args); + $self->{'_data'} = {}; + if( defined $alleles ) { + if( ref($alleles) !~ /HASH/i ) { + $self->warn("Must specify a valid HASH reference for the -alleles value...Ignoring initializing input"); + + } else { + foreach my $mkr ( keys %{$alleles} ) { + $self->add_alleles($mkr,@{$alleles->{$mkr}}); + } + } + } + return $self; +} + +=head2 add_alleles + + Title : add_alleles + Usage : $node->add_alleles($mkr,@alleles); + Function: Adds allele(s) for $mkr, @alleles can be a single or + multiple alleles. If the same marker is added more than one, the + previous value will be overwritten with a warning. + Returns : none + Args : $marker => marker name + @alleles => alleles for the marker + + +=cut + +sub add_alleles{ + my ($self,$marker,@alleles) = @_; + if( ! defined $marker || $marker eq '' ) { + $self->warn("must specify a valid marker name for add_alleles"); + return; + } + if( $self->{'_data'}->{$marker} ) { + $self->warn("Overwriting value of $marker"); + } + $self->{'_data'}->{$marker} = []; # reset the array ref + foreach my $a ( sort @alleles ) { + next if ! defined $a; # skip undef alleles + push @{$self->{'_data'}->{$marker}},$a; + } +} + +=head2 get_alleles + + Title : get_alleles + Usage : my @alleles = $node->get_alleles($marker); + Function: Return the alleles for a marker $marker + Returns : Array of Alleles for a marker or empty array + Args : $marker name + +=cut + +sub get_alleles{ + my ($self,$marker) = @_; + if( defined $self->{'_data'}->{$marker} ) { + return @{$self->{'_data'}->{$marker}}; + } + return (); +} + +=head2 get_marker_names + + Title : get_marker_names + Usage : my @names =$node->get_marker_names(); + Function: Return the names of the markers that have been added to this node + Returns : List of Marker Names + Args : none + +=cut + +sub get_marker_names{ + my ($self) = @_; + return keys %{$self->{'_data'}}; +} + +=head2 purge_markers + + Title : purge_markers + Usage : $node->purge_markers; + Function: Reset the markers and alleles + Returns : none + Args : none + + +=cut + +sub purge_markers{ + my ($self) = @_; + $self->{'_data'} = {}; + return; +} + + +1;