PhyloEM function

Model Estimation with Detection of Shifts