[R-sig-eco] Phylogenetic Independent Contrasts in R

JOSÉ Luiz luizecologia at gmail.com
Fri Nov 4 23:14:20 CET 2016


Dear list members,

I repeat to fail in running a phylogenetic independent Contrasts analysis
in the ape package. I'm new in this kind of analysis. Apparently, the
length of phenotypic and of phylogenetic data do not match.

I'm using the following script:

phy=read.newick(file.choose())#open my tree with adjusted branches
phytree=collapse.singles(phy)
phytree=multi2di(phytree)# remove polytomies in the tree randomly - only
dichotomous trees work in this package
is.binary.tree(phytree)

names(phytree)
[1] "edge"        "Nnode"       "tip.label"   "edge.length" "node.label"

data = read.csv(file.choose(), sep = ';')#trait data

names(data)
 [1] "LDMC"    "Lhum"    "Larea"   "SLA"     "Shum"    "Sdens"   "Bark"
 "Vdiam"   "Vdens"   "Vindex"
[11] "Slength" "Sdiam"   "Litter"  "Vlitter"

pic(data, phytree, scaled = TRUE, var.contrasts = FALSE)

Error in pic(data, phytree, scaled = TRUE, var.contrasts = FALSE) :
  length of phenotypic and of phylogenetic data do not match

the length of PHYTREE is 5
the length of DATA is 14

I don't understand how pic doesn't use species as length. I thank anyone
who can help.

José Luiz Silva

Brazil

	[[alternative HTML version deleted]]



More information about the R-sig-ecology mailing list