various updates...[21~

[ape.git] / R / DNA.R

diff --git a/R/DNA.R b/R/DNA.R
index b5f472be1ed6f53115220335dee71e00747f2052..f5c0a5947376b36d3f289ffa601b62a09b42020d 100644 (file)
--- a/R/DNA.R
+++ b/R/DNA.R
@@ -1,4 +1,4 @@
-## DNA.R (2008-01-19)
+## DNA.R (2008-02-01)
 
 ##   Manipulations and Comparisons of DNA Sequences
 
@@ -240,7 +240,7 @@ nuc.div <- function(x, variance = FALSE, pairwise.deletion = FALSE)
     ans <- sum(dist.dna(x, "raw", pairwise.deletion = pairwise.deletion))/
         (n*(n - 1)/2)
     if (variance) {
-        var <- (n + 1)*ans/(3*(n + 1)*s) + 2*(n^2 + n + 3)*ans/(9*n*(n - 1))
+        var <- (n + 1)*ans/(3*(n + 1)*dim(x)[2]) + 2*(n^2 + n + 3)*ans/(9*n*(n - 1))
         ans <- c(ans, var)
     }
     ans