man/bionj.Rd

   1 \name{BIONJ}
   2 \alias{bionj}
   3 \title{
   4   Tree Estimation Based on an Improved Version of the NJ Algorithm
   5 }
   6 \description{
   7   This function performs the BIONJ algorithm of Gascuel (1997).
   8 }
   9 \usage{
  10 bionj(X)
  11 }
  12 \arguments{
  13   \item{X}{a distance matrix; may be an object of class \code{"dist"}.}
  14 }
  15 \value{
  16   an object of class \code{"phylo"}.
  17 }
  18 \references{
  19   Gascuel, O. (1997) BIONJ: an improved version of the NJ algorithm
  20   based on a simple model of sequence data.
  21   \emph{Molecular Biology and Evolution}, \bold{14:}, 685--695.
  22 }
  23 \author{
  24   original C code by Hoa Sien Cuong and Olivier Gascuel; adapted and
  25   ported to \R by Vincent Lefort \email{vincent.lefort@lirmm.fr}
  26 }
  27 \seealso{
  28   \code{\link{nj}}, \code{\link{fastme}}, \code{\link{mvr}},
  29   \code{\link{SDM}}, \code{\link{dist.dna}}
  30 }
  31 \examples{
  32 ### From Saitou and Nei (1987, Table 1):
  33 x <- c(7, 8, 11, 13, 16, 13, 17, 5, 8, 10, 13,
  34        10, 14, 5, 7, 10, 7, 11, 8, 11, 8, 12,
  35        5, 6, 10, 9, 13, 8)
  36 M <- matrix(0, 8, 8)
  37 M[lower.tri(M)] <- x
  38 M <- t(M)
  39 M[lower.tri(M)] <- x
  40 dimnames(M) <- list(1:8, 1:8)
  41 tr <- bionj(M)
  42 plot(tr, "u")
  43 ### a less theoretical example
  44 data(woodmouse)
  45 trw <- bionj(dist.dna(woodmouse))
  46 plot(trw)
  47 }
  48 \keyword{models}
  49