ancr function

Compute marginal or joint ancestral state estimates