dportik / dadi_pipeline

An accessible and flexible tool for fitting demographic models with dadi using custom or published models (available here), conducting goodness of fit tests, and plotting.
GNU Lesser General Public License v3.0
63 stars 30 forks source link

nu1,nu2,nu3 #20

Closed pengyan19 closed 2 years ago

pengyan19 commented 2 years ago

how did i transform the nu1 to real Ne?,i found that nu1<1, and nu2>1,but i found that nu1+nu2=nuA, so in three population, nu1,nu2,nuA,nu3,how did i transform ?

dportik commented 2 years ago

Hello @pengyan19 , The dadi user group contains discussions about transforming parameters, and I would encourage you to double check this information there too.

Real population size = nu * Nref -nu is the inferred unscaled population size

Nref requires an additional calculation involving theta to get the reference population estimate:

Theta = 4 Nref mu * L

So you must first solve for Nref (based on theta), then use Nref value to scale nu.

Example: Theta = 4 Nref mu L 87.08 = Nref 4 0.00000001 289,873 87.08 = Nref * 0.01159492 7,510 = Nref