Re: [R-sig-phylo] Estimating ancestral states of a trait under selection

2013-04-20 Thread Eliot Miller
Hi Brian, Liam and Graham, Thanks a lot for these three suggestions! I'm actually going to go ahead and try all three approaches. Appreciate the help! Cheers, Eliot On Fri, Apr 19, 2013 at 4:17 PM, Slater, Graham slat...@si.edu wrote: Hi Eliot and others, You can also use the code

[R-sig-phylo] Estimating ancestral states of a trait under selection

2013-04-19 Thread Eliot Miller
Hello all, I've read some past posts and tried to figure this one out myself. It appears that the short answer is if there are not some extinct lineages in the tree with known trait values, it's difficult to do. I have a data set of the current climatic niches of a continental radiation of taxa.

Re: [R-sig-phylo] Estimating ancestral states of a trait under selection

2013-04-19 Thread Brian O'Meara
We have a method in R, TreEvo, that can allow you to do things like have a model where the optimum shifts through time according to some factor (like external data on climate) or where there are constraints that change through time: it's up on R-forge, and you can install the package, but it

Re: [R-sig-phylo] Estimating ancestral states of a trait under selection

2013-04-19 Thread Liam J. Revell
If you want to do a Bayesian reconstruction for continuous traits with a strong prior on the root - as Brian suggests - you can do this in the phytools function anc.Bayes. The prior means variances need to be provided in anc.Bayes(...,control=list(pr.var,pr.mean)). Unfortunately, you need to

Re: [R-sig-phylo] Estimating ancestral states of a trait under selection

2013-04-19 Thread Slater, Graham
Hi Eliot and others, You can also use the code contained in the fitContinuousMCMC package available here - https://www.eeb.ucla.edu/gradstudents/slater/Graham/code.html And described in this paper - http://onlinelibrary.wiley.com/doi/10./j.1558-5646.2012.01723.x/abstrac t To place a