Hi Vanderlei,
Try this which uses ape function prop.part():
temp<-prop.part(tree)
X<-matrix(0,nrow=length(tree$tip),ncol=length(temp),dimnames=list(tree$tip.label,tree$node.label))
for(i in 1:ncol(X)) X[temp[[i]],i]<-1
- Liam
--
Liam J. Revell
University of Massachusetts Boston
web: http://fac
Good morning!
I need to create a binary matrix with all node of a phylogenetic tree and the
presence of each taxo in their respective node.
Example:
require(ape)
y<-read.tree(text="(E,((H,I)D,(F,G)C)B)A;")
y
plot(y, show.node=TRUE)
I need to create a binary matrix as follows:
A