[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