Modify getMarginalBranchRateMatrix.R to return consistent output structure, and return appropriately named object when using bamm trait data