How to evaluate best fit model in Easy Mode (BM, OU, etc)?

Can we implement Easy Mode / Expert Mode code to evaluate what evolutionary model best fits a dataset, e.g., BM, OU, etc.?
This example is in R-phylo-wiki
Can we calculate multiple models and then provide the user with a likelihood ratio test or comparison of AICc scores, enabling users to select the best fit model.
