| line |
true |
false |
branch |
|
241
|
1 |
19 |
if defined $dir |
|
272
|
20 |
0 |
unless $self->{'_summary'} and not $self->{'_summary'}{'multidata'} |
|
276
|
16 |
4 |
if ($seqtype eq 'CODONML' or $seqtype eq 'AAML') { } |
|
|
2 |
2 |
elsif ($seqtype eq 'BASEML') { } |
|
|
2 |
0 |
elsif ($seqtype eq 'YN00') { } |
|
279
|
8 |
65 |
if ($seqtype eq 'CODONML' and /^pairwise comparison, codon frequencies:/) { } |
|
|
1 |
64 |
elsif ($seqtype eq 'AAML' and /^ML distances of aa seqs\.$/) { } |
|
|
9 |
55 |
elsif (/^Model\s+(\d+)/ or not $has_model_line and /^TREE/ and $seqtype eq 'CODONML' and not $self->{'_summary'}{'version'} =~ /4/) { } |
|
|
2 |
53 |
elsif (/for each branch/) { } |
|
|
3 |
50 |
elsif (/^TREE/) { } |
|
|
0 |
50 |
elsif (/Heuristic tree search by stepwise addition$/) { } |
|
|
0 |
50 |
elsif (/Heuristic tree search by NNI perturbation$/) { } |
|
|
0 |
50 |
elsif (/^stage 0:/) { } |
|
316
|
0 |
2 |
unless (defined $data{'-trees'}) |
|
322
|
0 |
2 |
unless ($tree and ref $tree and $tree->isa('Bio::Tree::Tree')) |
|
339
|
14 |
30 |
@nodes_L < 2 ? : |
|
343
|
0 |
44 |
unless ($n) |
|
346
|
30 |
14 |
unless ($n->is_Leaf and $n->id) |
|
407
|
2 |
1 |
if (/^Distances:/) { } |
|
|
1 |
0 |
elsif (/^TREE/) { } |
|
423
|
2 |
1 |
if (/^Estimation by the method|\(B\) Yang & Nielsen \(2000\) method/) |
|
433
|
20 |
0 |
if (%data) { } |
|
474
|
0 |
20 |
$self->{'_already_parsed_seqs'} ? : |
|
478
|
20 |
165 |
if (/^($SEQTYPES) \s+ # seqtype: CODONML, AAML, BASEML, CODON2AAML, YN00, etc
(?: \(in \s+ ([^\)]+?) \s* \) \s* )? # version: "paml 3.12 February 2002"; not present < 3.1 or YN00
(\S+) \s* # tree filename
(?: (.+?) )? # model description (not there in YN00)
\s* $ # trim any trailing space
/ox) { } |
|
|
6 |
159 |
elsif (/\s+?\d+?\s+?\d+?/ and $self->{'_already_parsed_seqs'} != 1) { } |
|
|
1 |
158 |
elsif (/^Data set \d$/) { } |
|
|
0 |
158 |
elsif (/^Before\s+deleting\s+alignment\s+gaps/) { } |
|
|
5 |
153 |
elsif (@lines >= 3 and $self->{'_already_parsed_seqs'} != 1) { } |
|
|
0 |
300 |
elsif (/Printing out site pattern counts/ and $self->{'_already_parsed_seqs'} != 1) { } |
|
490
|
8 |
12 |
unless (defined $self->{'_summary'}{'model'}) |
|
493
|
6 |
2 |
if ($model_line =~ /^Model:/) |
|
499
|
18 |
2 |
if defined $self->{'_summary'}{'model'} |
|
501
|
0 |
20 |
if $self->{'_summary'}{'seqtype'} eq 'AAMODEL' |
|
519
|
2 |
3 |
unless (/^\n$/) |
|
529
|
0 |
20 |
unless (defined $self->{'_summary'}{'seqtype'}) |
|
536
|
14 |
6 |
if ($seqtype eq 'CODONML') { } |
|
|
2 |
4 |
elsif ($seqtype eq 'AAML') { } |
|
|
0 |
4 |
elsif ($seqtype eq 'CODON2AAML') { } |
|
|
2 |
2 |
elsif ($seqtype eq 'BASEML') { } |
|
|
2 |
0 |
elsif ($seqtype eq 'YN00') { } |
|
585
|
10 |
29 |
if (/^((?:Codon frequencies)|(?:Site-class models))\s*:\s+(.+)/) { } |
|
|
4 |
25 |
elsif (/^\s+$/) { } |
|
|
16 |
9 |
elsif (/^ns\s+=\s+/ or /^Frequencies/) { } |
|
604
|
2 |
1493 |
if (/^Nei|\(A\) Nei/) |
|
605
|
14 |
1479 |
if $done |
|
606
|
557 |
922 |
if /^\s+/ |
|
608
|
866 |
56 |
unless $okay or /^Codon position x base \(3x4\) table\, overall/ |
|
610
|
42 |
14 |
if (s/^position\s+(\d+):\s+//) |
|
619
|
14 |
28 |
if $pos == 3 |
|
624
|
8 |
29 |
if (/^Nei\s\&\sGojobori|\(A\)\sNei-Gojobori/) |
|
628
|
8 |
21 |
if $done |
|
629
|
8 |
13 |
if (/^Codon frequencies under model, for use in evolver/) |
|
631
|
8 |
128 |
if /^\s+$/ |
|
645
|
0 |
2 |
unless $self->{'_summary'}{'seqs'} |
|
647
|
1 |
33 |
if (/^TREE/ or /^AA distances/) |
|
651
|
1 |
32 |
if $done |
|
652
|
12 |
20 |
if /^\s+$/ or /^\(Ambiguity/ |
|
653
|
2 |
18 |
if (/^Frequencies\./) { } |
|
|
0 |
18 |
elsif (not $okay) { } |
|
|
2 |
16 |
elsif (not $header) { } |
|
|
2 |
14 |
elsif (/^\#\s+constant\s+sites\:\s+
(\d+)\s+ # constant sites
\(\s*([\d\.]+)\s*\%\s*\)/x) { } |
|
|
1 |
13 |
elsif (/^ln\s+Lmax\s+\(unconstrained\)\s+\=\s+(\S+)/x) { } |
|
705
|
0 |
3 |
unless $self->{'_summary'}{'seqs'} |
|
708
|
2 |
21 |
if $done |
|
709
|
0 |
21 |
if (/^TREE/) |
|
710
|
2 |
19 |
if (/^\s+$/) |
|
711
|
0 |
2 |
if $seen |
|
714
|
3 |
16 |
if (/^(AA|ML) distances/) |
|
720
|
16 |
0 |
if ($okay) |
|
726
|
0 |
5 |
if ($type eq 'ML' and not @names and @vl) |
|
736
|
0 |
35 |
unless @vl |
|
743
|
3 |
13 |
if scalar @names == $numseqs |
|
753
|
0 |
86 |
unless defined $v |
|
771
|
0 |
18 |
if exists $self->{'_summary'}{'patterns'} |
|
774
|
11 |
88 |
if (/^Codon\s+(usage|position)/ or /Model/) { } |
|
|
52 |
36 |
elsif ($patternct) { } |
|
|
18 |
18 |
elsif (/^ns\s+\=\s*(\d+)\s+ls\s+\=\s*(\d+)/) { } |
|
|
7 |
11 |
elsif (/^\# site patterns \=\s*(\d+)/) { } |
|
781
|
7 |
45 |
if /^\s+$/ |
|
812
|
8 |
18 |
if $self->{'_already_parsed_seqs'} |
|
815
|
5 |
86 |
if (/^(Printing|After|TREE|Codon)/) |
|
819
|
13 |
73 |
if /^\s+$/ and @seqs > 0 |
|
820
|
6 |
67 |
if /^\s+$/ |
|
821
|
1 |
66 |
if /^\d+\s+$/ |
|
826
|
13 |
53 |
if (not @firstseq) { } |
|
851
|
13 |
5 |
if (@seqs > 0) |
|
867
|
8 |
16 |
if /^\s+$/ |
|
870
|
1 |
15 |
if (/^\(A\)\sNei-Gojobori\s\(1986\)\smethod/) |
|
873
|
1 |
8 |
if ($_ =~ /.*\d+?\.\d+?\s*\(.*/) |
|
886
|
14 |
2 |
if ($self->{'_summary'}{'seqtype'} eq 'CODONML') |
|
893
|
2 |
14 |
if ($self->{'_summary'}{'seqtype'} eq 'YN00') |
|
894
|
1 |
1 |
if ($firstseq) |
|
904
|
16 |
0 |
if /^\s+$/ and exists $self->{'_summary'}{'ngmatrix'} |
|
905
|
0 |
82 |
if /^\s+$/ or /^NOTE:/i |
|
909
|
12 |
70 |
if ($self->{'_summary'}{'seqtype'} eq 'YN00') { } |
|
917
|
14 |
68 |
unless defined $rest |
|
919
|
12 |
70 |
if ($self->{'_summary'}{'seqtype'} eq 'YN00') |
|
934
|
2 |
14 |
if ($self->{'_summary'}{'seqtype'} eq 'YN00' and @seqs) |
|
946
|
1 |
7 |
if ($self->{'_summary'}{'model'} =~ /kappa = (\d+?\.\d+?) fixed/) |
|
950
|
8 |
672 |
if (/^pairwise comparison, codon frequencies\:\s*(\S+)\./) { } |
|
|
112 |
560 |
elsif (/^(\d+)\s+\((\S+)\)\s+\.\.\.\s+(\d+)\s+\((\S+)\)/) { } |
|
|
112 |
448 |
elsif (/^lnL\s+\=\s*(\-?\d+(\.\d+)?)/) { } |
|
|
112 |
336 |
elsif (/^t\=\s*(\d+(\.\d+)?)\s+/) { } |
|
|
336 |
0 |
elsif (/^\s+$/) { } |
|
|
0 |
0 |
elsif (/^\s+(\d+\.\d+)\s+(\d+\.\d+)\s+(\d+\.\d+)/) { } |
|
962
|
112 |
0 |
if (defined($_ = $self->_readline)) |
|
968
|
1 |
111 |
if ($omega eq '') |
|
995
|
112 |
0 |
defined $t && length $t ? : |
|
|
112 |
0 |
defined $omega && length $omega ? : |
|
1023
|
2 |
9 |
if /^seq\.\s+seq\./ |
|
1026
|
43 |
8 |
if (/^\s+(\d+)\s+ # seq #
(\d+)\s+ # seq #
(\d+(\.\d+))\s+ # S
(\d+(\.\d+))\s+ # N
(\d+(\.\d+))\s+ # t
(\d+(\.\d+))\s+ # kappa
(\d+(\.\d+))\s+ # omega
\-??(\d+(\.\d+))\s+ # dN
\+\-\s+
\-??(\d+(\.\d+))\s+ # dN SE
\-??(\d+(\.\d+))\s+ # dS
\+\-\s+
\-??(\d+(\.\d+))\s+ # dS SE
/ox) { } |
|
|
4 |
4 |
elsif (/^\s+$/) { } |
|
|
1 |
3 |
elsif (/^\(C\) LWL85, LPB93 & LWLm methods/) { } |
|
1074
|
0 |
175 |
if $done |
|
1075
|
13 |
162 |
if (s/^TREE\s+\#\s*\d+:\s+//) { } |
|
|
11 |
151 |
elsif (/^Node\s+\&/ or /^\s+N37/ or /^(CODONML|AAML|YN00|BASEML)/ or /^\*\*/ or /^Detailed output identifying parameters/) { } |
|
|
13 |
138 |
elsif (/^tree\s+length\s+\=\s+(\S+)/) { } |
|
|
13 |
125 |
elsif (/^\s*lnL\(.+np\:\s*(\d+)\)\:\s+(\S+)/) { } |
|
|
26 |
99 |
elsif (/^\(/) { } |
|
|
3 |
96 |
elsif (/^SEs for parameters/) { } |
|
|
13 |
83 |
elsif (/^\s*\d+\.\.\d+/) { } |
|
1107
|
26 |
0 |
if ($tree) |
|
1110
|
13 |
13 |
if ($okay > 0) |
|
1113
|
13 |
0 |
if (not %match and @ids) |
|
1121
|
93 |
27 |
if ($match{$child}) { } |
|
1128
|
0 |
13 |
if ($self->verbose > 1) |
|
1137
|
3 |
10 |
if (defined $self->{'_SEs'}) |
|
1147
|
21 |
12 |
@nodes_L < 2 ? : |
|
1151
|
0 |
33 |
unless ($n) |
|
1185
|
5 |
115 |
if $done |
|
1186
|
47 |
68 |
if /^\s+$/ |
|
1187
|
0 |
9 |
unless $okay or /^Model\s+\d+/ or /^TREE/ |
|
1189
|
7 |
61 |
if (/^Model\s+(\d+)/) { } |
|
|
7 |
54 |
elsif (/^Time used\:\s+(\S+)/) { } |
|
|
9 |
45 |
elsif (m[^kappa\s+\(ts/tv\)\s+\=\s+(\S+)]) { } |
|
|
9 |
36 |
elsif (/^TREE/) { } |
|
|
1 |
35 |
elsif (m[^omega\s+\(dn/ds\)\s+\=\s+(\S+)]i) { } |
|
|
6 |
29 |
elsif (/^(Naive Empirical Bayes)|(Bayes Empirical Bayes)|(Positively\sselected\ssites)/i) { } |
|
|
15 |
14 |
elsif (/^dN/i) { } |
|
|
0 |
14 |
elsif (/^Parameters in beta:/) { } |
|
|
0 |
14 |
elsif (/^Parameters in beta\&w\>1:/) { } |
|
|
3 |
11 |
elsif (/^alpha\s+\(gamma\)\s+\=\s+(\S+)/) { } |
|
1190
|
0 |
7 |
if ($okay) { } |
|
1214
|
9 |
0 |
if (defined $data{'-trees'} and scalar @{$data{'-trees'};}) |
|
1247
|
7 |
8 |
if (/K\=(\d+)/) { } |
|
|
8 |
0 |
elsif (/for each branch/) { } |
|
1250
|
7 |
5 |
unless ($_ =~ /^\s+$/) |
|
1255
|
1 |
6 |
if (/^site class/) { } |
|
1290
|
0 |
8 |
unless (defined $data{'-trees'}) |
|
1296
|
0 |
8 |
unless ($tree and ref $tree and $tree->isa('Bio::Tree::Tree')) |
|
1313
|
40 |
72 |
@nodes_L < 2 ? : |
|
1317
|
0 |
112 |
unless ($n) |
|
1321
|
72 |
40 |
unless ($n->is_Leaf and $n->id) |
|
1335
|
0 |
0 |
if (/p\=\s+(\S+)\s+q\=\s+(\S+)/) { } |
|
1353
|
0 |
0 |
if (/p0\=\s+(\S+)\s+p\=\s+(\S+)\s+q\=\s+(\S+)/) { } |
|
1362
|
0 |
0 |
if (/\(p1\=\s+(\S+)\)\s+w\=\s*(\S+)/) { } |
|
1382
|
3 |
0 |
if (s/^r\s+\(\s*\d+\)\:\s+//) |
|
1386
|
3 |
0 |
if (s/^f\s*\:\s+//) |
|
1410
|
36 |
80 |
if /^\s+$/ or /^\s+Pr\(w\>1\)/ |
|
1411
|
4 |
76 |
if (/^Time used/ or /^TREE/) |
|
1415
|
5 |
71 |
if (/^Naive Empirical Bayes/i) { } |
|
|
2 |
69 |
elsif (/^Bayes Empirical Bayes/i) { } |
|
|
4 |
65 |
elsif (/^Positively selected sites/) { } |
|
|
9 |
56 |
elsif ($okay and /^\s+(\d+)\s+(\S+)\s+(\-?\d+(?:\.\d+)?)(\**)\s+(\-?\d+(?:\.\d+)?)\s+\+\-\s+(\-?\d+(?:\.\d+)?)/) { } |
|
|
8 |
48 |
elsif ($okay and /^\s+(\d+)\s+(\S+)\s+(\-?\d*(?:.\d+))(\**)\s+(\-?\d+(?:\.\d+)?)/) { } |
|
|
8 |
42 |
elsif ($okay and /^\s+(\d+)\s+(\S)\s+([\d\.\-\+]+)(\**)/) { } |
|
1429
|
0 |
9 |
unless defined $signif |
|
1437
|
0 |
8 |
unless defined $signif |
|
1442
|
0 |
8 |
unless defined $signif |
|
1455
|
38 |
98 |
if /^\s+$/ |
|
1456
|
0 |
98 |
unless $okay or /^\s+branch\s+t/ |
|
1457
|
10 |
88 |
if (/^\s+branch\s+(.+)/) { } |
|
|
78 |
10 |
elsif (/^\s*(\d+\.\.\d+)/) { } |
|
1486
|
0 |
2 |
unless $self->{'_summary'}{'seqs'} |
|
1488
|
0 |
28 |
if (/^TREE/ or /^Distances/) |
|
1489
|
2 |
26 |
if $done |
|
1490
|
10 |
16 |
if /^\s+$/ or /^\(Ambiguity/ |
|
1491
|
2 |
14 |
if (/^Frequencies\./) { } |
|
|
0 |
14 |
elsif (not $okay) { } |
|
|
2 |
12 |
elsif (not $header) { } |
|
|
2 |
10 |
elsif (/^\#\s+constant\s+sites\:\s+
(\d+)\s+ # constant sites
\(\s*([\d\.]+)\s*\%\s*\)/ox) { } |
|
|
2 |
8 |
elsif (/^ln\s+Lmax\s+\(unconstrained\)\s+\=\s+(\S+)/ox) { } |
|
1534
|
0 |
2 |
unless $self->{'_summary'}{'seqs'} |
|
1537
|
0 |
13 |
if (/^TREE/) |
|
1538
|
1 |
12 |
if $done |
|
1539
|
2 |
10 |
if /^This matrix is not used in later/ |
|
1540
|
2 |
8 |
if (/^\s+$/) |
|
1541
|
0 |
2 |
if $seen |
|
1544
|
2 |
6 |
if (/^Distances:(\S+)\s+\(([^\)]+)\)\s+\(alpha set at (\-?\d+\.\d+)\)/) |
|
1551
|
6 |
0 |
if ($okay) |
|
1555
|
4 |
2 |
if (defined $vl) |
|
1563
|
0 |
4 |
unless ($i) |
|
1570
|
2 |
4 |
if scalar @names == $numseqs |
|
1610
|
0 |
0 |
if (/^Rate\s+parameters:\s+/) { } |
|
|
0 |
0 |
elsif (/^Base\s+frequencies:\s+/) { } |
|
|
0 |
0 |
elsif (m[^Rate\s+matrix\s+Q,\s+Average\s+Ts/Tv\s+(\([^\)+]+\))?\s*\=\s+
(\-?\d+\.\d+)]x) { } |
|
|
0 |
0 |
elsif (/^alpha\s+\(gamma,\s+K=\s*(\d+)\s*\)\s*\=\s*(\-?\d+\.\d+)/) { } |
|
|
0 |
0 |
elsif (s/^(r|f):\s+//) { } |
|
1627
|
0 |
0 |
if /^\s+$/ |
|
1628
|
0 |
0 |
if (/^alpha/) |
|
1652
|
0 |
1 |
unless $self->{'_dir'} and -d $self->{'_dir'} and -r $self->{'_dir'} |
|
1655
|
0 |
1 |
unless -e $rstfile and not -z $rstfile |
|
1666
|
1 |
615 |
if (/^TREE\s+\#\s+(\d+)/) { } |
|
|
1 |
614 |
elsif (/^Prob\sof\sbest\scharacter\sat\seach\snode,\slisted\sby\ssite/) { } |
|
|
1 |
613 |
elsif (/^Check\sroot\sfor\sdirections\sof\schange\./ or /^Summary\sof\schanges\salong\sbranches\./) { } |
|
|
0 |
613 |
elsif (/^Overall\s+accuracy\s+of\s+the\s+(\d+)\s+ancestral\s+sequences:/) { } |
|
|
3 |
610 |
elsif (/^List of extant and reconstructed sequences/o) { } |
|
1668
|
1 |
9 |
if (/tree\s+with\s+node\s+labels\s+for/) |
|
1681
|
11 |
1 |
if (defined(my $blen = $n->branch_length)) |
|
1696
|
4 |
131 |
if /^Site/i or /^\s+$/ |
|
1697
|
130 |
1 |
if (s/^\s+(\d+)\s+(\d+)\s+([^:]+)\s*:\s*(.+)//) { } |
|
|
1 |
0 |
elsif (/^Summary\sof\schanges\salong\sbranches\./) { } |
|
1735
|
0 |
1 |
unless ($tree) |
|
1747
|
36 |
57 |
if /^\s+$/ |
|
1748
|
1 |
56 |
if (/^List\sof\sextant\sand\sreconstructed\ssequences/) { } |
|
|
11 |
45 |
elsif (/^Branch\s+(\d+):\s+(\d+)\.\.(\d+)\s+/) { } |
|
|
45 |
0 |
elsif (/^\s+(\d+)\s+([A-Z*])\s+(\S+)\s+\-\>\s+([A-Z*])\s+(\S+)?/) { } |
|
1756
|
0 |
11 |
unless ($node) |
|
1762
|
11 |
0 |
if (/\(n=\s*(\S+)\s+s=\s*(\S+)\)/) |
|
1773
|
0 |
45 |
unless ($node) |
|
1783
|
28 |
17 |
if defined $dprob |
|
1803
|
0 |
0 |
if (@overall_seq != @overall_site or @overall_seq != @seqs) |
|
1823
|
2 |
309 |
if /^Overall accuracy of the/ |
|
1824
|
61 |
248 |
if (/^\s+$/) |
|
1825
|
0 |
61 |
if $seqcount and $seqcount == @seqs |
|
1828
|
0 |
248 |
if (/^\s*(\d+)\s+(\d+)\s+$/) |
|
1833
|
20 |
228 |
if (/^node/) |
|
1837
|
1 |
19 |
if (not @firstseq) { } |