Re: [R-sig-phylo] p-value in pgls.Ives

2018-07-30 Thread Iara Reinaldo Coriolano
Thank you Liam! I'm gonna check the post! All the best, Iara Iara Reinaldo Coriolano Bacharela em Ciências Biológicas pela Universidade Federal do Ceará Mestranda pelo Programa de Pós-graduação em Ecologia - UFRGS Laboratório de Ecologia de Populações - UFRGS - http://ferrazlab.org/ 2018-07-30

Re: [R-sig-phylo] p-value in pgls.Ives

2018-07-30 Thread Liam J. Revell
Dear Iara. In fact, it is now possible to conduct a null hypothesis test with this method in phytools. I posted about it in my blog in 2017: http://blog.phytools.org/2017/11/bivariate-phylogenetic-regression-with.html. All the best, Liam Liam J. Revell, Associate Professor of Biology

Re: [R-sig-phylo] Simulating tree from a given tree

2018-07-30 Thread Liam J. Revell
Dear Miguel & Eduardo. This is correct, except that you should keep in mind that pbtree (like most, but not all, such simulators) generates trees from the first speciation event. Thus, the correct procedure must first involve drawing a wait time from an exponential distribution, deciding if

Re: [R-sig-phylo] Simulating tree from a given tree

2018-07-30 Thread Eduardo Ascarrunz
Maybe you could use a birth-death model. You can estimate the speciation rate and extinction rate from your original tree, and use those parameters to run simulations of trees for 5 My with the pbtree function from phytools. Run as many simulations as you have living tips in your original tree,

[R-sig-phylo] p-value in pgls.Ives

2018-07-30 Thread Iara Reinaldo Coriolano
Hello all, I'm using the function pgls.Ives() from the package phytools to incorporate within-species variance in both variables. I would like to know whether a statistic to evaluate the significance of the values were incorporated after the 2015 post asking about it. Thanks in advance, Iara