One way to do it with pGLS and root to tip path length may be to include the number of nodes along a root to tip path as a control for node density effects (maybe with some transformation). I don't know how well that will work, but I am also experimenting with similar questions.
Another way is to take the tip rates from a relaxed clock estimate of divergence times. CoEvol is much more powerful but as you note, not super practical. If you have a good number of species pairs, that may be the least problematic way to do it-- it sounds like you have a large tree so maybe that is optimal. there is an R package called sisters that may help with it. https://rdrr.io/github/bomeara/sisters/ Best, Jake Berv On Thu, Nov 17, 2022, 5:38 AM Karla Shikev <karlashi...@gmail.com> wrote: > Dear all, > > I want to test for a relationship between interspecific rates of evolution > and a continuous predictor variable. I didn't want to use species pairs > because they end up looking a lot of information. The method implemented in > Colevol is very elegant, but is computationally unfeasible with large > trees. > > One alternative could be to optimize BLs on a fixed topology and then use > root-to-tip distances as a measure of rate of evolution, but I'm concerned > about node density artifacts. > > Any suggestions would be greatly appreciated. > > Karla > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-phylo mailing list - R-sig-phylo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo > Searchable archive at > http://www.mail-archive.com/r-sig-phylo@r-project.org/ > [[alternative HTML version deleted]] _______________________________________________ R-sig-phylo mailing list - R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/