Dear Scott,
Having Icong in R would be nice but it is not the case
yet. Sorry...
However if someone has a code (in C for
example) to compute the maximum agreement subtree (MAST) between two
trees, I would be happy to do the rest in order to have the Icong index
within R (within the ape
Hi Damien Scott,
I don't know the formula of Icong or MAST, but I'd like to remind you
the function prop.part() that finds all the clades in a set of trees
and returns the compositions and frequencies, eg:
a - rtree(4)
b - rtree(4)
prop.part(c(a, b))
1: t2
2: t1
3: t3
4: t4
== 2