Peter
unread,May 18, 2012, 10:47:33 AM5/18/12Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to Isolation with Migration
This IMa2 output header at the top of the 2Nm posterior estimates has
me confused because it seems inconsistent with the manual:
HISTOGRAM GROUP 4: POPULATION MIGRATION (2NM) POSTERIOR PROBABILITY
HISTOGRAMS
-----------------------------------------------------------------------------
curve height is an estimate of the posterior probability
each term is the product of a population parameter (e.g. q0) and a
migration rate (e.g.m0>1)
migration rates are in the coalescent (backwards in times), so that a
population migration rate of
q1m0>1 is the population rate (forward in time) at which population 1
receives migrants from population 0
IMa2 Manual p. 12: For example, m0>1 is the name of the migration rate
parameter for migration of genes from population 0 to 1, backwards in
time, in the coalescent direction. Interpreted forward in time, in the
usual way that most people think about migration, m0>1 is the rate at
which population 0 receives genes from population 1.
I'm stuck on the last line of the header from the output because if
q0m0>1 is the population migration rate for migrants going from
population 0 to population 1 backwards in time, I feel like that
header should read:
q1m1>0 is the population rate (forward in time) at which population 1
receives migrants from population 0.