Add as.phylo.bammdata() to getRateThroughTimeMatrix.R to account for lack of phylo class in ephy object.