Skip to content

Commit

Permalink
sync with main trunk
Browse files Browse the repository at this point in the history
svn path=/bioperl-live/branches/branch-1-6/; revision=16147
  • Loading branch information
cjfields committed Sep 22, 2009
1 parent 89fe0dd commit 086fcff
Show file tree
Hide file tree
Showing 19 changed files with 1,620 additions and 1,615 deletions.
123 changes: 76 additions & 47 deletions Bio/Align/DNAStatistics.pm
Expand Up @@ -73,21 +73,37 @@ in brackets are the pattern which will match
=over 3
=item JukesCantor [jc|jukes|jukescantor|jukes-cantor]
=item *
=item Uncorrected [jcuncor|uncorrected]
JukesCantor [jc|jukes|jukescantor|jukes-cantor]
=item F81 [f81|felsenstein]
=item *
=item Kimura [k2|k2p|k80|kimura]
Uncorrected [jcuncor|uncorrected]
=item Tamura [t92|tamura|tamura92]
=item *
=item F84 [f84|felsenstein84]
F81 [f81|felsenstein]
=item TajimaNei [tajimanei|tajima\-nei]
=item *
=item JinNei [jinnei|jin\-nei] (not implemented)
Kimura [k2|k2p|k80|kimura]
=item *
Tamura [t92|tamura|tamura92]
=item *
F84 [f84|felsenstein84]
=item *
TajimaNei [tajimanei|tajima\-nei]
=item *
JinNei [jinnei|jin\-nei] (not implemented)
=back
Expand All @@ -104,7 +120,7 @@ several pre-requisites for the alignment.
=item 1
DNA alignment must be based on protein alignment. Use the subroutine
L<aa_to_dna_aln> in Bio::Align::Utilities to achieve this.
L<Bio::Align::Utilities/aa_to_dna_aln> to achieve this.
=item 2
Expand Down Expand Up @@ -140,49 +156,53 @@ comparisons in an MSA. The statistics returned are:
=over 3
=item S_d
=item *
Number of synonymous mutations between the 2 sequences.
S_d - Number of synonymous mutations between the 2 sequences.
=item N_d
=item *
Number of non-synonymous mutations between the 2 sequences.
N_d - Number of non-synonymous mutations between the 2 sequences.
=item S
=item *
Mean number of synonymous sites in both sequences.
S - Mean number of synonymous sites in both sequences.
=item N
=item *
mean number of synonymous sites in both sequences.
N - mean number of synonymous sites in both sequences.
=item P_s
=item *
proportion of synonymous differences in both sequences given by P_s = S_d/S.
P_s - proportion of synonymous differences in both sequences given by
P_s = S_d/S.
=item P_n
=item *
proportion of non-synonymous differences in both sequences given by P_n = S_n/S.
P_n - proportion of non-synonymous differences in both sequences given
by P_n = S_n/S.
=item D_s
=item *
estimation of synonymous mutations per synonymous site (by Jukes-Cantor).
D_s - estimation of synonymous mutations per synonymous site (by
Jukes-Cantor).
=item D_n
=item *
estimation of non-synonymous mutations per non-synonymous site (by Jukes-Cantor).
D_n - estimation of non-synonymous mutations per non-synonymous site (by
Jukes-Cantor).
=item D_n_var
=item *
estimation of variance of D_n .
D_n_var - estimation of variance of D_n .
=item D_s_var
=item *
estimation of variance of S_n.
D_s_var - estimation of variance of S_n.
=item z_value
=item *
calculation of z value.Positive value indicates D_n E<gt> D_s,
z_value - calculation of z value.Positive value indicates D_n E<gt> D_s,
negative value indicates D_s E<gt> D_n.
=back
Expand All @@ -191,25 +211,25 @@ The statistics returned by calc_average_KaKs are:
=over 3
=item D_s
=item *
Average number of synonymous mutations/synonymous site.
D_s - Average number of synonymous mutations/synonymous site.
=item D_n
=item *
Average number of non-synonymous mutations/non-synonymous site.
D_n - Average number of non-synonymous mutations/non-synonymous site.
=item D_s_var
=item *
Estimated variance of Ds from bootstrapped alignments.
D_s_var - Estimated variance of Ds from bootstrapped alignments.
=item D_n_var
=item *
Estimated variance of Dn from bootstrapped alignments.
D_n_var - Estimated variance of Dn from bootstrapped alignments.
=item z_score
=item *
calculation of z value. Positive value indicates D_n E<gt>D_s,
z_score - calculation of z value. Positive value indicates D_n E<gt>D_s,
negative values vice versa.
=back
Expand All @@ -222,7 +242,6 @@ the book, and reproduce those results. If people like having this sort
of analysis in BioPerl other methods for estimating Ds and Dn can be
provided later.
Much of the DNA distance code is based on implementations in EMBOSS
(Rice et al, www.emboss.org) [distmat.c] and PHYLIP (J. Felsenstein et
al) [dnadist.c]. Insight also gained from Eddy, Durbin, Krogh, &
Expand All @@ -232,26 +251,36 @@ Mitchison.
=over 3
=item D_JukesCantor
=item *
D_JukesCantor
"Phylogenetic Inference", Swoffrod, Olsen, Waddell and Hillis, in
Mol. Systematics, 2nd ed, 1996, Ch 11. Derived from "Evolution of
Protein Molecules", Jukes & Cantor, in Mammalian Prot. Metab., III,
1969, pp. 21-132.
=item D_Tamura
=item *
D_Tamura
K Tamura, Mol. Biol. Evol. 1992, 9, 678.
=item D_Kimura
=item *
D_Kimura
M Kimura, J. Mol. Evol., 1980, 16, 111.
=item JinNei
=item *
JinNei
Jin and Nei, Mol. Biol. Evol. 82, 7, 1990.
=item D_TajimaNei
=item *
D_TajimaNei
Tajima and Nei, Mol. Biol. Evol. 1984, 1, 269.
Expand Down
16 changes: 8 additions & 8 deletions Bio/AnalysisI.pm
Expand Up @@ -198,8 +198,8 @@ sub describe { shift->throw_not_implemented(); }
The analysis input data are named, and can be also associated with a
default value, with allowed values and with few other attributes. The
names are important for feeding the service with the input data (the
inputs are given to methods C<create_job>, C<run>, and/or C<wait_for>
as name/value pairs).
inputs are given to methods C<create_job>, C<Bio::AnalysisI|run>, and/or
C<Bio::AnalysisI|wait_for> as name/value pairs).
Here is a (slightly shortened) example of an input specification:
Expand Down Expand Up @@ -324,8 +324,8 @@ tool.
Call this method if you wish to "stage the scene" - to create a job
with all input data but without actually running it. This method is
called automatically from other methods (C<run> and C<wait_for>) so
usually you do not need to call it directly.
called automatically from other methods (C<Bio::AnalysisI|run> and
C<Bio::AnalysisI|wait_for>) so usually you do not need to call it directly.
The input data and prameters for this execution can be specified in
various ways:
Expand Down Expand Up @@ -459,24 +459,24 @@ sub id { shift->throw_not_implemented(); }

# -----------------------------------------------------------------------------

=head2 run
=head2 Bio::AnalysisI::JobI::run
Usage : $job->run
Returns : itself
Args : none
It starts previously created job. The job already must have all input
data filled-in. This differs from the method of the same name of the
C<Bio::Tools::Run::Analysis> object where the C<run> method creates
also a new job allowing to set input data.
C<Bio::Tools::Run::Analysis> object where the C<Bio::AnalysisI::JobI::run> method
creates also a new job allowing to set input data.
=cut

sub run { shift->throw_not_implemented(); }

# -----------------------------------------------------------------------------

=head2 wait_for
=head2 Bio::AnalysisI::JobI::wait_for
Usage : $job->wait_for
Returns : itself
Expand Down
2 changes: 1 addition & 1 deletion Bio/Assembly/Tools/ContigSpectrum.pm
Expand Up @@ -785,7 +785,7 @@ sub average {
}


=head2 average
=head2 score
Title : score
Usage : my $score = $csp->score();
Expand Down

0 comments on commit 086fcff

Please sign in to comment.