Hi Nick, I hope you are ok. I have run six different models (DEC, DIVALIKE, BAYAREALIKEA and their J+ versions) and I got that BAYAREALIKE+J was the best fit model. However, this one shows quite strange results.
A large number of nodes have as their most probable ancestral range ranges composed of 3 areas. This seems strange to me because many of the extant species do not have that range of 3 areas (only a few). In addition, these are areas that, with the geographic information that we currently have, do not make much sense. I thought that this would be a result of the BAYAREALIKE model allowing only widespread sympatry and not other processes such as subset sympatry and vicariance, therefore, it is assumed that the ancestral ranks are identical to those of the descendants (although they are few, there are some that have a range of 3 areas).
I work with species that are distributed on different islands (each species on no more than one island) and that have probably dispersed by sea between them. Before testing the BAYARELIKE (y +J) model, I performed the model comparison between DEC (y +J) and DIVALIKE (y +J). In this case, the DEC model had the best fit, but both DEC and the other models yielded results that did make a lot of sense, taking into account the geographical and biological history of the group and its distribution. My question is, is it possible to exclude a model because, based on how it's built, it doesn't "fit" my biological group? Is there any information regarding this about the BAYAREALIKE?
Thank you so much!