tag:blogger.com,1999:blog-8499895524521663926.post7015603177917320294..comments2018-09-18T10:24:43.482-04:00Comments on Phylogenetic Tools for Comparative Biology: Fitting a variable-process model of discrete character evolution on the tree using phytoolsLiam Revellhttp://www.blogger.com/profile/04314686830842384151noreply@blogger.comBlogger6125tag:blogger.com,1999:blog-8499895524521663926.post-64013835572946724962018-05-23T10:42:55.853-04:002018-05-23T10:42:55.853-04:00It was thinking about whether I could utilize this... It was thinking about whether I could utilize this review on my other site, I will connect it back to your site though.Great Thanks. <a href="https://www.kplokusa.com" rel="nofollow">click here</a><br />Matiashttps://www.blogger.com/profile/02810391078111197553noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-7849918202195116612017-12-06T22:40:58.189-05:002017-12-06T22:40:58.189-05:00To add to your discussion. I used corHMM's ray...To add to your discussion. I used corHMM's raydisk function to do an extension of the covarion model with very large Q-matrices. You can do LRT there and not only type I error but full power analyses. The pruning algorithm in corHMM is really efficient allowing you to work with large phylogenies. What I added to raydisk was an efficient yet decently accurate way to calculate exp{Qt}. I know phytools uses expm from package matrix but in the case of large and sparse Q matrices (>20) is nor efficient nor accurate.Rosanahttps://www.blogger.com/profile/05017413675490511626noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-88617348450948220742017-12-06T15:53:26.230-05:002017-12-06T15:53:26.230-05:00Incidentally (and this is a cool feature added by ...Incidentally (and this is a cool feature added by Jeremy Beaulieu, not me), the rayDISC function in corHMM is one of the coolest unknown tools in phylogenetics to me. Once you convert multiple characters (say, three binary traits) into a single multistate character (000 = state 0, ..., 111 = state 7) you can create very flexible rate matrices. So you can say "I can go from 000 to 001 at the same rate as 000 to 010, but 000 -> 100 is forbidden, and...."). I bring it up here because if you had a case where you thought one trait affected evolution of another, rather than stochastically mapping one and then using that to fit regimes, you probably get a better fit by inferring the mapping and the rates jointly (so you have a matrix that has the 0->1 rate for trait 1 depending on the state of trait 2, etc.). [see https://academic.oup.com/sysbio/article/62/2/339/1668230 for why the mapping one trait then estimating rate of the other might be less good than doing it jointly. ;-)]Brian O'Mearahttps://www.blogger.com/profile/02450575101806778825noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-17821992172721293302017-12-06T15:46:04.678-05:002017-12-06T15:46:04.678-05:00Beaulieu et al.'s corhmm is an extension of th...Beaulieu et al.'s corhmm is an extension of the covarion model that does what John suggests. It does rate reconstruction and ancestral state estimates at nodes (and of course it's then straightforward do make a stochastic character mapping version). <br /><br />It's an interesting comparison of Liam's approach vs corhmm. Liam's requires mapping a regime on the tree, and having the rate matrix change based on the regime mapping. Corhmm doesn't require the mapping ahead of time, in effect allowing it to be inferred (the placement of the hidden state). It's like Brownie (or brownie.lite) and OUwie requiring pre-mapped regimes rather than auteur and SURFACE inferring the regime shifts [for standardized testing fans, Liam new method is to corHMM as OUwie is to SURFACE]. I can see biological use cases for both the pre-mapped regime case and the hidden regime case.Brian O'Mearahttps://www.blogger.com/profile/02450575101806778825noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-23100863737566637052017-12-05T22:17:23.496-05:002017-12-05T22:17:23.496-05:00Hi John.
I believe this is what Beaulieu et al. ...Hi John. <br /><br />I believe this is what Beaulieu et al. (<a href="https://donoghuelab.yale.edu/sites/default/files/beaulieu_systbiol_2013_0.pdf" rel="nofollow">2013</a>) published, and they specifically point out that it is a generalization of the covarion model. <br /><br />Please correct me if I'm wrong, of course, but this is not the same, so far as I can tell - in the sense that the regimes we are computing the likelihood over are fixed. That means, for instance, that we might have one time period (from t=0 to t=0.25, say) on the tree that we permit to evolve under one process, than a second under a different process, and so on. That is - the regime 'paintings' are observed or specifically hypothesized a priori. By contrast in the covarion model & its guild we have an unobserved set of regimes & we compute the likelihood by integrating over the probability that each datum comes from each regime.<br /><br />Note that this is not better or more sophisticated than the covarion or 'hidden-rates' model - in fact, it is less! However, it does suit a specific class of hypothesis that is common in phylogenetic comparative biology - for instance that trait evolution differs between clades or among geological eras, and allows us to contrast that hypothesis against one in which it does not differ or varies in a different way.<br /><br />Thanks for pointing out the relationship to the covarion model. This is a very important literature, of course!<br /><br />- LiamLiam Revellhttps://www.blogger.com/profile/04314686830842384151noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-87092644678930735852017-12-05T21:56:54.078-05:002017-12-05T21:56:54.078-05:00This seems like an overly complicated solution. Wh...This seems like an overly complicated solution. Why not just make a covarion-like model? That is to say, you embed different 2 X 2 rate matrices into a (2 X N) X (2 X N) rate matrix, where N is the number of "regimes." The 2 X 2 rate matrices go along the diagonal. Along the off diagonal 2 X 2 areas, you have a rate of switching from one regime to another. As an aside, this is not new at all. It's an extension of the covariant model. My colleagues and I did something identical, though computationally more intensive, for selection regimes: we embedded three codon models into a 183 X 183 rate matrix, allowing switching among selection regimes. <br /><br />Guindon, S., A. G. Rodrigo, K. A. Dyer, and J. P. Huelsenbeck. 2004. Modeling the site-specific variation of selection patterns along lineages. {\it Proceedings of the National Academy of Sciences, U.S.A.} 101(35):12957--12962.John Huelsenbeckhttps://www.blogger.com/profile/02319608745646179087noreply@blogger.com