fix bug in prop.clades

[ape.git] / R / CDAM.global.R

diff --git a/R/CDAM.global.R b/R/CDAM.global.R
index 184bffaf0cfc18372d41abc68e32efd9052fd645..5644ca5ad3af09ebfb47dd79f1f950a9bbec1ab6 100644 (file)
--- a/R/CDAM.global.R
+++ b/R/CDAM.global.R
@@ -143,7 +143,7 @@
        Rmat <- apply(Y,2,rank)
 
        ## Correction factors for tied ranks (eq. 3.3)
-       t.ranks <- apply(Rmat, 2, function(x) summary(as.factor(x)))
+       t.ranks <- apply(Rmat, 2, function(x) summary(as.factor(x), maxsum=nd))
        TT <- sum(unlist(lapply(t.ranks, function(x) sum((x^3)-x))))
        # if(!silent) cat("TT = ",TT,'\n')