annotate variant_effect_predictor/Bio/Seq/SeqFactory.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: SeqFactory.pm,v 1.8 2002/10/25 22:49:04 lapp Exp $
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
2 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
3 # BioPerl module for Bio::Seq::SeqFactory
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::Seq::SeqFactory - Instantiates a new Bio::PrimarySeqI (or derived class) through a factory
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::Seq::SeqFactory;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
20 my $factory = new Bio::Seq::SeqFactory;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
21 my $seq = $factory->create(-seq => 'WYRAVLC',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
22 -id => 'name');
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
23
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
24 # If you want the factory to create Bio::Seq objects instead
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
25 # of the default Bio::PrimarySeq objects, use the -type parameter:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
26
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
27 my $factory = new Bio::Seq::SeqFactory(-type => 'Bio::Seq');
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
28
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
29
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
30 =head1 DESCRIPTION
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
31
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
32 This object will build L<Bio::PrimarySeqI> and L<Bio::SeqI> objects
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
33 generically.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
34
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
35 =head1 FEEDBACK
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
36
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
37 =head2 Mailing Lists
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
38
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
39 User feedback is an integral part of the evolution of this and other
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
40 Bioperl modules. Send your comments and suggestions preferably to
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
41 the Bioperl mailing list. Your participation is much appreciated.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
42
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
43 bioperl-l@bioperl.org - General discussion
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
44 http://bioperl.org/MailList.shtml - About the mailing lists
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
45
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
46 =head2 Reporting Bugs
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
47
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
48 Report bugs to the Bioperl bug tracking system to help us keep track
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
49 of the bugs and their resolution. Bug reports can be submitted via
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
50 email or the web:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
51
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
52 bioperl-bugs@bioperl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
53 http://bugzilla.bioperl.org/
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
54
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
55 =head1 AUTHOR - Jason Stajich
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
56
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
57 Email jason@bioperl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
58
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
59 Describe contact details here
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
60
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
61 =head1 CONTRIBUTORS
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
62
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
63 Additional contributors names and emails here
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
64
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
65 =head1 APPENDIX
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
66
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
67 The rest of the documentation details each of the object methods.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
68 Internal methods are usually preceded with a _
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
69
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
70 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
71
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
72
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
73 # Let the code begin...
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
74
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
75
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
76 package Bio::Seq::SeqFactory;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
77 use vars qw(@ISA);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
78 use strict;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
79
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
80 use Bio::Root::Root;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
81 use Bio::Factory::SequenceFactoryI;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
82
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
83 @ISA = qw(Bio::Root::Root Bio::Factory::SequenceFactoryI);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
84
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
85 =head2 new
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
86
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
87 Title : new
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
88 Usage : my $obj = new Bio::Seq::SeqFactory();
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
89 Function: Builds a new Bio::Seq::SeqFactory object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
90 Returns : Bio::Seq::SeqFactory
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
91 Args : -type => string, name of a PrimarySeqI derived class
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
92 This is optional. Default=Bio::PrimarySeq.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
93
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
94 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
95
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
96 sub new {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
97 my($class,@args) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
98 my $self = $class->SUPER::new(@args);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
99 my ($type) = $self->_rearrange([qw(TYPE)], @args);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
100 if( ! defined $type ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
101 $type = 'Bio::PrimarySeq';
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
102 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
103 $self->type($type);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
104 return $self;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
105 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
106
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
107
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
108 =head2 create
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
109
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
110 Title : create
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
111 Usage : my $seq = $seqbuilder->create(-seq => 'CAGT', -id => 'name');
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
112 Function: Instantiates new Bio::SeqI (or one of its child classes)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
113 This object allows us to genericize the instantiation of sequence
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
114 objects.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
115 Returns : Bio::PrimarySeq object (default)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
116 The return type is configurable using new(-type =>"...").
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
117 Args : initialization parameters specific to the type of sequence
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
118 object we want. Typically
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
119 -seq => $str,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
120 -display_id => $name
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
121
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
122 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
123
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
124 sub create {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
125 my ($self,@args) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
126 return $self->type->new(-verbose => $self->verbose, @args);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
127 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
128
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
129 =head2 type
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
130
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
131 Title : type
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
132 Usage : $obj->type($newval)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
133 Function:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
134 Returns : value of type
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
135 Args : newvalue (optional)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
136
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
137
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
138 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
139
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
140 sub type{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
141 my ($self,$value) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
142 if( defined $value) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
143 eval "require $value";
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
144 if( $@ ) { $self->throw("$@: Unrecognized Sequence type for SeqFactory '$value'");}
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
145
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
146 my $a = bless {},$value;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
147 unless( $a->isa('Bio::PrimarySeqI') ||
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
148 $a->isa('Bio::Seq::QualI') ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
149 $self->throw("Must provide a valid Bio::PrimarySeqI or Bio::Seq::QualI or child class to SeqFactory Not $value");
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
150 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
151 $self->{'type'} = $value;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
152 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
153 return $self->{'type'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
154 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
155
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
156 1;