Re: [R-sig-phylo] compar.ou

2013-11-12 Thread sandra goutte
Dear all, Following up this conversation; Using OUwie, i get reasonable values for theta under several models for all my traits, but one. In this case, this is a trait which can only take positive values and i get negative theta values for all the different models, including BM and OU1. the

[R-sig-phylo] Accommodating multiple data points per species in pGEE

2013-11-12 Thread Megan Welsford
Dear list users, I'm analyzing a data set consisting of several floral traits in dioecious *Leucadendron*. I would like to use pGEE to account for phylogentic relatedness among species. I'm particularly interested in interactions between plant sex and pollination system (wind or insect). However,

Re: [R-sig-phylo] match two tree edges

2013-11-12 Thread Jingchun Li
Thank you both! Liam's code worked for my case. Best, Jingchun On Mon, Nov 11, 2013 at 10:21 PM, Emmanuel Paradis emmanuel.para...@ird.frwrote: Hi, Have a look at the function makeNodeLabel with the option method = md5sum: it will create a label to each node depending on the tips

Re: [R-sig-phylo] compar.ou

2013-11-12 Thread Brian O'Meara
Often if you have a trait constrained to be positive, it's appropriate to log-transform it, which has the happy effect of making the theta for the trait when converted back to a non-log scale constrained to be positive (as well as probably being more appropriate for how your trait evolves). You