tag:blogger.com,1999:blog-8499895524521663926.post8512938765987339656..comments2019-09-19T00:31:54.140-04:00Comments on Phylogenetic Tools for Comparative Biology: Testing a hypothesis that phylogenetic signal is different from 1.0Liam Revellhttp://www.blogger.com/profile/04314686830842384151noreply@blogger.comBlogger7125tag:blogger.com,1999:blog-8499895524521663926.post-25252648545919319242018-10-13T11:24:17.949-04:002018-10-13T11:24:17.949-04:00Hi Liam and phytools users,
I have a question abo...Hi Liam and phytools users,<br /><br />I have a question about the interpretation of phylogenetic signal (specifically, pagel's lambda). I calculated the phylogenetic signal over multiple trees and for each tree I obtained the value of lambda and to evaluate if lambda is significantly different from 0 (with test=TRUE in phylosig). Additionally I performed a likelihood ratio test to evaluate in each case if the phylogenetic signal is different from 1 (as in this post). <br /><br />Overall I got that the phylogenetic signal over all trees has a mean of 0.89 (0.8-0.96) and is always significantly different from 0. However, the p-value when comparing 1 is sometimes much lower than 0.05 and sometimes higher. Would you still consider this as a 'high' phylogenetic signal even if deviates from 1 in some cases? I appreciate all your opinions and insight!<br /><br />Thank you all very much for your time and help!<br /><br />Best,<br /><br />LauraLaura Cespedeshttps://www.blogger.com/profile/16625533422053049422noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-59916223337690340902017-04-06T17:55:01.226-04:002017-04-06T17:55:01.226-04:00Hello Mr. Revell,
I work with Ecology and Evoluti...Hello Mr. Revell,<br /><br />I work with Ecology and Evolution in Tropical Forests and I'm still few familiar with Pagel's Lambda and the function provided by "phytools" to calculate it and I have a doubt in some results I got here.<br /><br />For instance, in one of results I got: Lamba =1.00033; logL<br />= -769.7819; logL0 = -1313.418 and P-value of pchisq to the likelihood test = 1.928349e-238.<br /><br />In principle, I thought these values a bit wierd. So, I checked everything In my phylo tree and trait matrix and apparentely everything was correct, even because I perfectly ran the phylo signal test using Blomberg's K.<br /><br />Thus, my doubt is centered exactly in the fact Lambda and p value above. Lamba is subtly > 1 and p - value is signiticant. In my opinion the lambda is quite close to entire 1 and p value very high. I mean, lambda seems more consistent with Brownian Motion (Lambda = 1) and p-value is significant. Perhaps, it's a stupid misleading I'm doing, but one lamba so close to one like should present significant p-value in the likelihood test? <br /><br />I read some publications to better interpret this result, but I remained with this doubt: could I consider this result like significant phylo signal or more congruent with Brownian Motion?<br /><br />Thank you very much,<br /><br />ÉcioÉcio Dinizhttps://www.blogger.com/profile/05689419984378141651noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-11266621595901967032014-12-13T17:18:24.129-05:002014-12-13T17:18:24.129-05:00Nope, it did not. Thanks Liam!Nope, it did not. Thanks Liam!Marcelo Arayahttps://www.blogger.com/profile/11306025292691573119noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-24659718381405030332014-12-08T11:28:58.828-05:002014-12-08T11:28:58.828-05:00Hola Marcelo.
Since this looks like it is trying ...Hola Marcelo.<br /><br />Since this looks like it is trying to invert the Hessian matrix, it may just be involved in computing standard errors for the parameter estimates (and not the optimization). Did the function still return an optimized λ in this case.<br /><br />Saludos, LiamLiam Revellhttps://www.blogger.com/profile/04314686830842384151noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-62551962358311568462014-12-07T09:30:52.281-05:002014-12-07T09:30:52.281-05:00thanks Liam. I was able to run this for most varia...thanks Liam. I was able to run this for most variables but in some cases I got this error when applying the brownie.lite function:<br /><br />Error in solve.default(model1$hessian) :<br />Lapack routine dgesv: system is exactly singular: U[2,2] = 0<br /><br />Do you know what could be the problem?<br /><br />MarceloMarcelo Arayahttps://www.blogger.com/profile/11306025292691573119noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-2610955996176955552014-12-07T09:29:37.903-05:002014-12-07T09:29:37.903-05:00This comment has been removed by the author.Marcelo Arayahttps://www.blogger.com/profile/11306025292691573119noreply@blogger.comtag:blogger.com,1999:blog-8499895524521663926.post-44370504778575580472014-12-02T17:40:12.485-05:002014-12-02T17:40:12.485-05:00Note that a couple of caveats apply - particularly...Note that a couple of caveats apply - particularly to testing hypothesis about λ using the likelihood ratio:<br /><br />(1) The likelihood-ratio test, in which we use a Χ^2distribution for hypothesis testing, depends on asymptotic properties of likelihoods which do not hold for small sample sizes (and small trees).<br /><br />(2) There is also a problem when we are at the bounds of a particular parameter - and this will be true near 0.0 & 1.0 for λ.<br /><br />For more details, see Carl Boettiger's paper which I believe can be found <a href="http://dx.doi.org/10.1111/j.1558-5646.2011.01574.x" rel="nofollow">here</a>.Liam Revellhttps://www.blogger.com/profile/04314686830842384151noreply@blogger.com