search

popgenmethods / smcpp

SMC++ infers population history from whole-genome sequence data.
GNU General Public License v3.0
149 stars 32 forks source link

generation time differs #232

Open xuefenfei712 opened 1 year ago

xuefenfei712 commented 1 year ago

Hi, Thanks for your wonderful application, I use it to estimate Ne in my populations, however, the generation time for it was not fixed, for example before domestication g=10, after it will be 12+, if I want to use SMC++ plot to draw it, what should I do, thankyou!

Best Xue

Devashish13 commented 1 year ago

Hi @xuefenfei712 , As per my knowledge generation time(g) is constant, If you use SMC++ plot command.

Intuitively I think you plot using SMC++ taking time(in generations) and produce the .csv file containing Ne and time(in g). 1) Generate one column converting time in generation to actual time keeping g = 12+ up until domestication then change to g =10 2) Now plot Ne against actual time you have.

You need somebody expert to confirm this!!