[adegenet-commits] r621 - pkg/man
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Mon May 10 17:47:18 CEST 2010
Author: jombart
Date: 2010-05-10 17:47:18 +0200 (Mon, 10 May 2010)
New Revision: 621
Modified:
pkg/man/seqTrack.Rd
Log:
Have to test the doc. Example not done yet (have to format the data).
Modified: pkg/man/seqTrack.Rd
===================================================================
--- pkg/man/seqTrack.Rd 2010-05-10 15:34:05 UTC (rev 620)
+++ pkg/man/seqTrack.Rd 2010-05-10 15:47:18 UTC (rev 621)
@@ -108,12 +108,42 @@
}
\details{
=== Maximum parsimony genealogies ===
+ Maximum parsimony genealogies can be obtained easily using this
+ implementation of seqTrack. One has to provide in \code{x} a matrix of
+ genetic distances. The most straightforward distance is the number of
+ differing nucleotides. See \code{\link[ape]{dist.dna}} in the ape
+ package for a wide range of genetic distances between aligned
+ sequences. The argument \code{best} should be set to "min" (its
+ default value), so that the identified genealogy minimizes the total
+ number of mutations. If \code{x} contains number of mutations, then
+ \code{mu} and \code{haplo.length} should also be provided for
+ resolving ties in equally parsimonious ancestors using maximum
+ likelihood.
+
+ === Likelihood of observed genetic differentiation ===
+ The probability of oberving a given number of mutations between a
+ sequence and its ancestor can be computed using
+ \code{\link{get.likelihood.seqTrack}}. Note that this is only possible
+ if \code{x} contained number of mutations.
+
- === Likelihood of observed genetic differentiation ===
=== Converting seqTrack objects to graphs ===
+ seqTrack objects can be converted to \linkS4class{graphNEL} objects,
+ which can in turn be plotted and manipulated using classical graph
+ tools. Simply use 'as(x, "graphNEL")' where 'x' is a seqTrack
+ object. This functionality requires the \code{graph} package. Note
+ that this is to be installed from Bioconductor, likely using the following
+ command lines:
+ source("http://bioconductor.org/biocLite.R")
+ biocLite("graph")
-
+ Also note that the R package Rgraphviz (also on Bioconductor) provides
+ nice ways of plotting graphs (replace 'graph' with 'Rgraphviz' in the
+ previous command lines to install this package).
}
\seealso{
- \code{\link{get.likelihood.seqTrack}}
+ \code{\link[ape]{dist.dna}} to compute pairwise genetic distances in aligned sequences.
+}
+\example{
+
}
\ No newline at end of file
More information about the adegenet-commits
mailing list