Olá Aline, muito legal teu trabalho e os resultados. Achei bem interessante o boxplot porque mostra que há maior variância nos ambientes lóticos, isso concorda com a sua previção em relação à heterogeneidade dos ambientes. Acho que seria legal se vcs comentassem a respeito disso tb e não só sobre a média, pode ser um aspecto interessante do trabalho.  Uma dúvida que fiquei é, normalmente a regressão logística é feita quando temos um preditor contínuo. No seu caso era um fator com dois níveis, vcs pensaram em fazer outros testes, como o threshold do Felsenstein ou mesmo correlações entre traits binários no corHMM ou correlações no BayesTrait? 

Olá, Aline. Na verdade, eu só gostaria de elogiar o seu trabalho! Achei super interessante mesmo e você apresentou super bem! Sucesso no mestrado :)