modification of maxlik.betasplit function from apTreeshape package that uses lbeta instead of beta; does not currently work when confidence.interval = 'profile'