This is an internal function to estimate the character value at the root node
ancState(phy, y)
a phylogeny in APE phylo format
A matrix of trait data corresponding to species in phy