 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 two (fake) Swiss towns                                           
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 +                                                                +
 +   POPULATION SIZE, MIGRATION, DIVERGENCE, ASSIGNMENT, HISTORY  +
 +   Bayesian inference using the structured coalescent           +
 +                                                                +
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  Compiled for a PARALLEL COMPUTER ARCHITECTURE
  One master and 7 compute nodes are available.
  Compiled for a SYMMETRIC multiprocessors (GrandCentral)
  PDF output enabled [Letter-size]
  Version 4.2.8   [June-24-2016]
  Program started at   Wed Jul 27 18:17:53 2016
                                                                        


Options in use:
---------------

Analysis strategy is BAYESIAN INFERENCE

Proposal distribution:
Parameter group          Proposal type
-----------------------  -------------------
Population size (Theta)  Metropolis sampling
Migration rate      (M)  Metropolis sampling


Prior distribution (Proposal-delta will be tuned to acceptance frequence 0.440000):
Parameter group            Prior type   Minimum    Mean(*)    Maximum    Delta
-------------------------  ------------ ---------- ---------- ---------- ----------
Population size (Theta_1)      Uniform  0.000000   0.050000   0.100000   0.010000 
Population size (Theta_2)      Uniform  0.000000   0.050000   0.100000   0.010000 
Migration 1 to 2 (M)      Uniform  0.000000  500.000000 1000.00000 100.000000




Inheritance scalers in use for Thetas (specified scalars=1)
1.00 1.00 1.00 
[Each Theta uses the (true) inheritance scalar of the first locus as a reference]


Pseudo-random number generator: Mersenne-Twister                                
Random number seed (with internal timer)           2639607152

Start parameters:
   First genealogy was started using a random tree
   Start parameter values were generated
Connection matrix:
m = average (average over a group of Thetas or M,
s = symmetric migration M, S = symmetric 4Nm,
0 = zero, and not estimated,
* = migration free to vary, Thetas are on diagonal
d = row population split off column population
D = split and then migration
   1 Aadorf         * c 
   2 Bern           * * 



Mutation rate is constant for all loci

Markov chain settings:
   Long chains (long-chains):                              1
      Steps sampled (inc*samples*rep):            2500000000
      Steps recorded (sample*rep):                    500000
   Combining over replicates:                             10
   Static heating scheme
      4 chains with  temperatures
       1.00, 1.50, 3.00,1000000.00
      Swapping interval is 1
   Burn-in per replicate (samples*inc):            250000000

Print options:
   Data file:                                  twoswisstowns
   Haplotyping is turned on:                              NO
   Output file (ASCII text):         outfile-twoswisstowns_c
   Output file (PDF):            outfile-twoswisstowns_c.pdf
   Posterior distribution:                         bayesfile
   Print data:                                            No
   Print genealogies:                                     No

Summary of data:
Title:                                two (fake) Swiss towns
Data file:                                     twoswisstowns
Datatype:                                     Haplotype data
Number of loci:                                            3
Mutationmodel:
 Locus  Sublocus  Mutationmodel   Mutationmodel parameter
-----------------------------------------------------------------
     1         1 Felsenstein 84  [Bf:0.26 0.27 0.22 0.25, t/t ratio=2.000]
     1         2 Felsenstein 84  [Bf:0.24 0.27 0.22 0.27, t/t ratio=2.000]
     2         1 Felsenstein 84  [Bf:0.26 0.23 0.25 0.25, t/t ratio=2.000]
     3         1 Felsenstein 84  [Bf:0.26 0.24 0.24 0.26, t/t ratio=2.000]


Sites per locus
---------------
Locus    Sites
     1     500 500
     2     300
     3     700

Population                   Locus   Gene copies    
----------------------------------------------------
  1 Aadorf                       1        10
  1 Aadorf                       1        10
  1                              2        10
  1                              3        10
  2 Bern                         1        10
  2 Bern                         1        10
  2                              2        10
  2                              3        10
    Total of all populations     1        40
                                 2        20
                                 3        20

