Add ability of getMarginalBranchRateMatrix.R to handle multiphylo input in addition to BAMM-data objects.