MCbaseRandomTree.dat

   1  0        * 0: paml format (mc.paml); 1:paup format (mc.nex)\r
   2 1234567   * random number seed (odd number)\r
   3 \r
   4 50 10000 100  * <# seqs>  <# nucleotide sites>  <# replicates>\r
   5 \r
   6 0.1 0.2 0.3 1.5   * birth rate, death rate, sampling fraction, and mutation rate (tree height)\r
   7 \r
   8 7          * model: 0:JC69, 1:K80, 2:F81, 3:F84, 4:HKY85, 5:T92, 6:TN93, 7:REV\r
   9 10 5 1 2 3 * kappa or rate parameters in model\r
  10 0.5  4     * <alpha>  <#categories for discrete gamma>\r
  11 \r
  12 0.1 0.2 0.3 0.4    * base frequencies\r
  13   T   C   A   G\r
  14 \r
  15 \r
  16 //end of file\r
  17 ==================================================\r
  18 The rest of this data file are notes, ignored by evolver.\r
  19 Change the values of the parameters, but do not delete them.\r
  20 evolver simulates nucleotide sequences under the REV+Gamma model\r
  21 and its simpler forms.\r
  22 \r
  23 More notes:\r
  24   Parameter kappa or rate parameters in the substituton model: \r
  25     For TN93, two kappa values are required, while for REV, 5 values\r
  26     (a,b,c,d,e) are required (see Yang 1994 for the definition of these\r
  27     parameters).  \r
  28     The kappa parameter is defined differently under HKY85 (when k=1 means\r
  29     no transition bias) and under F84 (when k=0 means no bias).\r
  30     JC69 and F81 are considered species cases of HKY85, so use 1 for kappa\r
  31     for those two models.  Notation is from my two papers in JME in 1994.\r
  32   \r
  33   Use equal base frequencies (0.25) for JC69 and K80.\r
  34   Use 0 for alpha to have the same rate for all sites.\r
  35   Use 0 for <#categories for discrete gamma> to use the continuous gamma\r
  36 \r
  37 =========!! Check screen output carefully !! =====\r
  38 \r
  39 \r
  40 A few interesting trees: \r
  41 For TipDate models:\r
  42 ((A@1970: 0.3, B@1990: 0.5) @1940 :0.4, (C@1940 : 0.1, (D@1960: 0.1, (E@1980:0.1, F@2000: 0.3) @1970 :0.2) @1950 : 0.2) @1930 :0.3) @1900;\r
  43 (((A :0.1, B :0.2) :0.12, C :0.3) :0.123, D :0.4, E :0.5) ;\r
  44 (((human: 0.1, chimpanzee: 0.1): 0.1, gorilla: 0.2): 0.1, orangutan: 0.3);\r
  45 ((a:0.1, b:0.2):0.12, (c:0.3, d:0.4):0.34, ((e:0.5, f:0.6):0.56, (g:0.7, h:0.8):0.78):0.69);\r