new alex()

[ape.git] / man / nj.Rd

diff --git a/man/nj.Rd b/man/nj.Rd
index 51983589a8ad4238a578c374181ae71176f8b6e0..88d2258c5893c5da384b74136ac357635d4b66d9 100644 (file)
--- a/man/nj.Rd
+++ b/man/nj.Rd
@@ -18,6 +18,10 @@ nj(X)
   Saitou, N. and Nei, M. (1987) The neighbor-joining method: a new
   method for reconstructing phylogenetic trees. \emph{Molecular Biology
     and Evolution}, \bold{4}, 406--425.
+
+  Studier, J. A. and Keppler, K. J. (1988) A note on the
+  neighbor-joining algorithm of Saitou and Nei. \emph{Molecular Biology
+    and Evolution}, \bold{5}, 729--731.
 }
 \author{Emmanuel Paradis}
 \seealso{
@@ -31,9 +35,10 @@ x <- c(7, 8, 11, 13, 16, 13, 17, 5, 8, 10, 13,
        10, 14, 5, 7, 10, 7, 11, 8, 11, 8, 12,
        5, 6, 10, 9, 13, 8)
 M <- matrix(0, 8, 8)
-M[row(M) > col(M)] <- x
-M[row(M) < col(M)] <- x
-rownames(M) <- colnames(M) <- 1:8
+M[lower.tri(M)] <- x
+M <- t(M)
+M[lower.tri(M)] <- x
+dimnames(M) <- list(1:8, 1:8)
 tr <- nj(M)
 plot(tr, "u")
 ### a less theoretical example