-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathGene_ontology.pm
executable file
·75 lines (59 loc) · 1.68 KB
/
Gene_ontology.pm
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
package Gene_ontology;
use strict;
sub new {
my $packagename = shift;
my ($go_id, $go_descript, $go_type) = @_;
my $self = { go_id=>$go_id,
go_type=>$go_type,
go_descript=>$go_descript,
go_evidences=>[], #list of go_evidence objects.
date => undef,
assignby=>undef,
qualifier => undef
};
bless ($self, $packagename);
return ($self);
}
sub set_go_info {
my ($self, $go_id, $descript, $type) = @_;
$self->{go_id} = $go_id;
$self->{go_descript} = $descript;
$self->{go_type} = $type;
return ($self);
}
sub get_evidence {
my $self = shift;
return (@{$self->{go_evidences}});
}
sub add_evidence {
my ($self, $ev_code, $evidence, $with_ev) = @_;
my $evidence_obj = Gene_ontology::Evidence->new($ev_code, $evidence, $with_ev);
push (@{$self->{go_evidences}}, $evidence_obj);
}
sub toString() {
my $self = shift;
my @evidence = $self->get_evidence();
my $text = "GO assignment: ID: $self->{go_id}, Type: $self->{go_type}, Descript: $self->{go_descript}\n, Qualifier: $self->{qualifier}, Assignedby: $self->{assignby}";
$text .= "\tEvidence:\n";
foreach my $evidence (@evidence) {
$text .= "\t\t" . $evidence->toString() . "\n";
}
return ($text);
}
###############################
package Gene_ontology::Evidence;
use strict;
sub new {
my ($packagename, $ev_code, $evidence, $with_ev) = @_;
my $self = { ev_code=>$ev_code,
evidence=>$evidence,
with_ev=>$with_ev };
bless ($self, $packagename);
return ($self);
}
sub toString() {
my ($self) = @_;
my $text = "ev_code: $self->{ev_code}, evidence: $self->{evidence}, with_ev: $self->{with_ev}";
return ($text);
}
1; #EOM