Dropping random tips and adjusting branch lengths and branching order for ML phylogeny (many times)
1
0
Entering edit mode
6 months ago

Hi all,

I would like to obtain many "sub-samples" of a Maximum Likelihood tree by randomly removing tips (about 20-30% of them from the complete tree). However, unlike the drop.tips function that comes with ape, I was hoping that the branch lengths and branching order of the resulting tree would be re-adjusted so that it too has maximized likelihood. This would be the same result as a maximum likelihood tree constructed from a smaller subset of sequences - although ideally faster, as it must be done many (~100). Are there any tree-building tools that accomplish this? Perhaps through some efficient use of information obtained during the construction of the complete tree?

*I'm not super intrested in Bootstrap support values (although if there's a reasonably fast way to get those too, I'd be impressed and happy to use it)

*I'm aiming to work with sets of sequences numbering in the low to mid 1000s (~1.3 Kb long)

thanks for any help that can be provided!

sequence R • 164 views
ADD COMMENT
0
Entering edit mode
6 months ago
Mensur Dlakic ★ 10k

Why not estimate the trees from scratch each time you drop some sequences? You may be interested in FastTree which I think is named aptly. There is another program called VeryFastTree but I haven't tried it.

ADD COMMENT

Login before adding your answer.

Traffic: 1478 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6