We used the R package BioGeoBEARS  for ancestral range estimation (ARE) on the MCC tree from the BEAST run under the optimal clock model and tree speciation prior. Recently, Ree & Sanmartín  demonstrated that the likelihood-based models with the +J parameter are invalid because of errors in the estimation of likelihoods. Here we compared the following three models of biogeographical estimation in the maximum likelihood (ML) framework: dispersal-extinction cladogenesis (DEC) model , dispersal–vicariance analysis (DIVA)  and BayArea model . The fit for the different models was assessed using the Akaike information criterion scores.