some news for ape 3.0-8

[ape.git] / R / ewLasso.R

diff --git a/R/ewLasso.R b/R/ewLasso.R
new file mode 100644 (file)
index 0000000..bd8f7fe
--- /dev/null
+++ b/R/ewLasso.R
@@ -0,0 +1,23 @@
+## ewLasso.R (2013-03-30)
+
+##   Lasso Tree
+
+## Copyright 2013 Andrei-Alin Popescu
+
+## This file is part of the R-package `ape'.
+## See the file ../COPYING for licensing issues.
+
+ewLasso <- function(X, phy)
+{
+    if (is.matrix(X)) X <- as.dist(X)
+    X[is.na(X)] <- -1
+    X[X < 0] <- -1
+    X[is.nan(X)] <- -1
+
+    N <- attr(X, "Size")
+    labels <- attr(X, "Labels")
+    if (is.null(labels)) labels <- as.character(1:N)
+    ans <- .C("ewLasso", as.double(X), as.integer(N),
+              phy$edge[, 1], phy$edge[, 2],
+              DUP = FALSE, NAOK = TRUE, PACKAGE = "ape")
+}