From f0ba00490fd11146cab4a4039b0055217d1b376f Mon Sep 17 00:00:00 2001
From: paradis <paradis@6e262413-ae40-0410-9e79-b911bd7a66b7>
Date: Mon, 18 Mar 2013 10:29:37 +0000
Subject: [PATCH] some changes done before the HD failure of my laptop... I
 presume

git-svn-id: https://svn.mpl.ird.fr/ape/dev/ape@210 6e262413-ae40-0410-9e79-b911bd7a66b7
---
 NEWS    |  6 ++++++
 R/pic.R | 18 ++++++++----------
 2 files changed, 14 insertions(+), 10 deletions(-)

diff --git a/NEWS b/NEWS
index 9bc3f6c..adf3082 100644
--- a/NEWS
+++ b/NEWS
@@ -21,6 +21,12 @@ BUG FIXES
       distance matrix and itself.
 
 
+OTHER CHANGES
+
+    o The files CDAM.global.R and CDAM.post.R have been renamed
+      CADM.global.R and CADM.post.R
+
+
 
 		CHANGES IN APE VERSION 3.0-7
 
diff --git a/R/pic.R b/R/pic.R
index bc54067..2735f8f 100644
--- a/R/pic.R
+++ b/R/pic.R
@@ -1,26 +1,24 @@
-## pic.R (2012-11-20)
+## pic.R (2013-02-18)
 
 ##   Phylogenetically Independent Contrasts
 
-## Copyright 2002-2012 Emmanuel Paradis
+## Copyright 2002-2013 Emmanuel Paradis
 
 ## This file is part of the R-package `ape'.
 ## See the file ../COPYING for licensing issues.
 
 pic <- function(x, phy, scaled = TRUE, var.contrasts = FALSE, rescaled.tree = FALSE)
 {
-    if (!inherits(phy, "phylo"))
-      stop("object 'phy' is not of class \"phylo\"")
-    if (is.null(phy$edge.length))
-      stop("your tree has no branch lengths")
+    if (!inherits(phy, "phylo")) stop("object 'phy' is not of class \"phylo\"")
+    if (is.null(phy$edge.length)) stop("your tree has no branch lengths")
     nb.tip <- length(phy$tip.label)
     nb.node <- phy$Nnode
     if (nb.node != nb.tip - 1)
-      stop("'phy' is not rooted and fully dichotomous")
+        stop("'phy' is not rooted and fully dichotomous")
     if (length(x) != nb.tip)
-      stop("length of phenotypic and of phylogenetic data do not match")
+        stop("length of phenotypic and of phylogenetic data do not match")
     if (any(is.na(x)))
-      stop("missing data in 'x': you may consider removing the species with missing data from your tree with the function 'drop.tip'.")
+        stop("missing data in 'x': you may consider removing the species with missing data from your tree with the function 'drop.tip'.")
 
     phy <- reorder(phy, "postorder")
     phenotype <- numeric(nb.tip + nb.node)
@@ -65,7 +63,7 @@ pic.ortho <- function(x, phy, var.contrasts = FALSE, intra = FALSE)
 {
     n <- length(x)
     m <- n - 1L # number of nodes
-    phy <- reorder(phy, "pruningwise")
+    phy <- reorder(phy, "postorder")
     xx <- unlist(lapply(x, mean)) # 'x' in Felsenstein's paper
     xx <- c(xx, numeric(m))
     delta.v <- numeric(n + m)
-- 
2.39.2