X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=man%2Fnj.Rd;h=88d2258c5893c5da384b74136ac357635d4b66d9;hb=d1546ec66ff1a8ea123adefebe14f6316c23705f;hp=51983589a8ad4238a578c374181ae71176f8b6e0;hpb=21eb56120c84786502f24ff9c27b39d5badfe1f7;p=ape.git

diff --git a/man/nj.Rd b/man/nj.Rd
index 5198358..88d2258 100644
--- 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