Question: User-defined tree in MrBayes
0
gravatar for victoria_aleks
2.5 years ago by
victoria_aleks30 wrote:

Dear all, I have a primitive question, I am afraid. I want to run MrBayes on amino acid sequences with the predefined tree, as random trees run does not reach conversion. Nevertheless, I can not find the command on how to do it. All the manuals and tutorials I have found say "User-defined trees can be read into MrBayes by executing a Nexus file with a ”trees” block" . That is followed by 'startvals' and 'npert'. I understand how to set the parameters for using the tree, but how to introduce the tree to Mr Bayes?

shortly: HOW can I READ IN the existing tree in MrBayes?

*'usertree' command resulted in the error message: "Usertree command deprecated. Define the tree in a treeblock and use 'Startvals' instead. Error in command "Usertree".

Thank you!

bayesian phylogeny mrbayes • 1.4k views
ADD COMMENTlink modified 2.5 years ago by h.mon31k • written 2.5 years ago by victoria_aleks30
2
gravatar for h.mon
2.5 years ago by
h.mon31k
Brazil
h.mon31k wrote:

See What do I do when it is difficult to get convergence from MrBayes FAQ:

begin trees;
  tree funkytree = ((runer:19.19959,newborn:6.80041):0.84600,((deacon:11.99700, neomorph:12.00300):7.52973,((protomorph:100.85930,facehugger:47.14069):20.59201, chestbuster:18.87953):2.09460):3.87382,praetorian:25.46154);
end trees;
begin mrbayes;
   startvals tau = funkytree V = funkytree;
end;

edit1: corrected extra semi-colon on tree.

edir2: as per victoria_aleks comment, I edited the tree block such as 'startvals tau = mytree' is in one line.

ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by h.mon31k

thank you! but then it says: "Could not find taxon 'Prot_ID' in list of taxa". what would that mean?

also, when I type: ... startvals tau = mytree -- it doesnt find the command 'tau', neither 'topology'

when I type: ... startvals tau = mytree --it doesnt know what is 'mytree' :(

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by victoria_aleks30

Probably formatting issues - for example, I just noticed my tree above has an extra semi-colon: praetorian:25.46154); ;. Also, it seems Prot_ID is found either at the tree and / or at the data blocks, but not at the taxa block.

Start with a small file - use the wikipedia toy example and build from that. Once you get the formatting right, then edit your real-world data set and try it.

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by h.mon31k

thanks, good idea. I spent already quite some time on all these small things...:) if i may ask: do you use MB often? My dataset is quite big (180 protein seq., ~100000 amino acids), and I have read that Mr Bayes is not good for big datasets...

ADD REPLYlink written 2.5 years ago by victoria_aleks30
1

It has been a long time since I've used MrBayes, but indeed seems like a big data set, and should take a long time to run. Probably ExaBayes is a better option for this data set.

ADD REPLYlink written 2.5 years ago by h.mon31k

it finally worked when I wrote 'startvals tau = mytree' in one line :) thank you again!

ADD REPLYlink written 2.5 years ago by victoria_aleks30

Glad it worked. Please don't forget to follow up on your threads: if an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they work. Upvote|Bookmark|Accept

Please do the same for your previous posts as well.

ADD REPLYlink written 2.5 years ago by h.mon31k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1293 users visited in the last hour