Closed drabosky closed 10 years ago
Can you git pull and try again? Without seeing the event data I'm guessing that it was a constant rate root process and that R was storing the lambda shift parameter as an integer. I modified getEventData to cast all lambda and mu parameters to double.
haven't looked into this but had found that plot.bammdata did not work when I input a BAMM run that had 0 rate shifts (e.g. only the root process for each sample). the R message was:
Error in dtRates(ephy, 0.01) : REAL() can only be applied to a 'numeric', not a 'integer'
we have probably not checked these but all BAMMtools fxns need to work with "single regime" output (or return an appropriate error message).