User:Carl Boettiger/Notebook/Comparative Phylogenetics/2010/03/02

{| width="800"
 * style="background-color: #EEE"|[[Image:owwnotebook_icon.png|128px]] Comparative Phylogenetics
 * style="background-color: #F2F2F2" align="center"|  |Main project page
 * style="background-color: #F2F2F2" align="center"|  |Main project page


 * colspan="2"|
 * colspan="2"|

Power in trees

 * Back to Phylogenetics after a week with only population dynamics.

Exploration

 * data created by brown simulation on the star tree in ouch matches the expected variance, $$\sigma^2t$$. This variance is slightly reduced when simulated on the Felsenstein tree, as might be intuitively expected due to the increased correlations.  The degree of decrease should be analytically attainable, but isn't obvious to me at the moment.  Example using functions in felsenstein_tree.R, revision 210:

E_trait_var simulates reps replicates under the given model and computes the variance across the tips of each dataset generated, and averages this variance across the replicates.


 * Meanwhile, data created by the hansen simulation on the star tree falls significantly short of $$\sigma^2/2\alpha$$ even when rates are certainly high enough to have reached the stationary distribution.

not sure why this is. The simulations seems to think the steadystate variances is $$\sigma^2/2\alpha^2$$ instead. Using the deprecated function hansen.dev, I seem to recover the expected variance (averaging variance over 100 replicates):


 * ouch has a convenient function to bootstrap a model


 * ouch hansen trees give false impressions of accuracy around star (and hence nearly star) trees, where result is actually very sensitive to the original initialized search parameters for alpha and sigma.

Brainstorming

 * Current methods don't permit divergent selection models, for instance of parts of the tree. Again probably an information-limited arena, but interesting to think what one would do or do differently if you had geographic information or range overlaps along the tree (or at least for the tips).


 * }