new standard RAxML version with marginal ancestral state computation

388 views
Skip to first unread message

Alexis

unread,
Sep 16, 2011, 10:49:50 AM9/16/11
to raxml
Dear All,

I just pushed a new version of standard RAxML to the GIT that computes
marginal ancestral states (see Ziheng Yang's book for details).

For marginal ancestral state reconstructions you need to pass a rooted
tree to RAxML.

To invoke, type for instance:

./raxmlHPC-SSE3 -f A -t testTree -s testData -m BINGAMMA -n T1

with the small test alignment:

4 4
t1 0101
t2 0111
t3 0100
t4 0110

and a test tree:

((t1,t2),(t3,t4));

RAxML will then output a rooted binary tree again with inner node
labels and two files
containing the marginal probabilities for each inner root label as
well as a guess for the ancestral sequence.

The code should work for all data types, partitioned datasets etc.

Once again, latest code updates are available at:

https://github.com/stamatak/standard-RAxML

Have fun,

Alexis



Reply all
Reply to author
Forward
0 new messages