Modify and merge two R functions: getSpeciesRateThroughTime.R and getSpeciesRateThroughTimeReturnAll.R