• Home
• About Us
  News Archive Partners Contact Us Acknowledgments
• Systems
  maya Other systems
• Research
  List of Projects Publications Posters Presentation Handout Supporting Materials
• Resources for HPCF Users
  System description Using your account Tutorial: Compiling C programs Scheduling rules Tutorial: Running programs Other Languages and Hardwares Other Packages General Productivity Usage rules Account request form Help request

Introduction

MrBayes is a program for Bayesian inference and model choice across a wide range of phylogenetic and evolutionary models. To run the MrBayes software interactively on the front-end node of maya, use the command We will now demonstrate running MrBayes on the compute nodes. This was adpated from MbWiki. First, make sure to update your switcher settings as follows This is necessary to match the compiler and MPI implementation originally used to configure MrBayes. Next create a .nex file as follows, with your contents instead of the placeholders. Now create a small script with commands to execute the .nex file Note that there are other ways to set up a call to MrBayes; for example, MCMC and likelihood options can be specified in a file separate from the data. The last step is to create a usual SLURM batch script. This script will run the MrBayes program on the compute nodes; the program will be run in parallel if multiple processes are requested. MyBayes will then run batch.txt, which will in turn execute our .nex file. Make sure you have read the how to run tutorial before attempting to use the batch system.

Setting checkpoints for long runs

For users needing very long runs of MrBayes, it is suggested to break up the work into several small jobs rather than one very long job. Long jobs have a higher probably of being interrupted by maintenance windows or unforeseen problems. Fortunately, MrBayes has a built in mechanism for creating checkpoints, where progress can be saved from one job and continued in a subsequent job.

To demonstrate this, consider the "primates.nex" example that comes with the MrBayes software. We will create two MB run scripts to analyze this data. The first script represents the initial run.

Notice that we set "checkfreq", which represents the number of generations before checkpointing. The second script continues where the first script left off. The only differences are the addition of the option "append=yes", which tells MrBayes to continue from the checkpoint, and that "ngen" has been increased to request additional generations. In this case, we want it to read the data from our 10-million generations run, and stop after 20 million generations. Note that all values, other than "ngen" and "append", must match between the two scripts; the run may fail otherwise. Running these two scripts directly from the command line yields the following This mechanism can be used to plan for very long runs, by setting "ngen" accordingly, and also to recover from unexpected failures. When carrying out real runs on maya, the user should place these calls into batch scripts as illustrated in the previous section.