arives / rr2

R2 for correlated data; see https://academic.oup.com/sysbio/advance-article/doi/10.1093/sysbio/syy060/5098616
GNU General Public License v3.0
18 stars 7 forks source link

phyloglm() alpha > 50 #14

Closed matthieu-haudiquet closed 4 years ago

matthieu-haudiquet commented 4 years ago

Hi,

Thanks for this great package.

I noticed this warning when running the rr2:R2() function on the results of phyloglm() :

Warning message:
In R2.lik.phyloglm(mod, mod.r) :
  In mod, estimate of alpha >50, so model refit with glm()

But the problem is that an alpha > 50 does not always mean that there is no phylogenetic signal (as stated in the doc of rr2), see : https://github.com/lamho86/phylolm/issues/28

Best