I'm trying to run a scenario with two populations. I think it might have got something to do with the conditions set on the time and Ne parameters, but I'm stuck. Thanks in advance!
Scenario:
N1 N2
0 sample 1
0 sample 2
t1-db varNe 1 N1b
t2 varNe 2 N2
t3-db varNe 2 N2b
t4 merge 1 2
Conditions:
N2>N1b, N2>N1, N1bN2, t4>t2, t4>t3, t4>t1, t3>t2, t2>t1
Error:
Something happened during the reftable generation :
Program of thread 'guanacos reference table generation' exited (with return code 1) unsuccessfully.
I will try to save current RNGs' states.
I have saved current RNGs' states into /Users/ana_agapito/Documents/Documentos/ubb/paper MSS/lamaguanicoe2/guanacos_2018_8_3-1/RNG_state_0000.bin
Content of progress file :
Just wondering what t1-db varNe 1 N1b mean?
And did you specify 0 samples from pop 1 and pop 2? This might cause the error as well since there are no samples that will coalesce.
Hi Ana,
In the scenario you specify "t2 varNe 2 N2" but "N2" is already the effective population size specified on line 1. So this line is not really changing the Ne of population 2 like you might want it to. Try changing it's name to see if this solves your problem. Also make sure the naming of all variables is consistent. You have t1-db in the scenario but just t1 in the conditions.
Also, for your conditions you can also remove some redundancies. For example, t3>t2 and t4>t3 imply that t4>t2 and so you could remove this condition. For the same reasons N2>N1b and t4>t1 are redundant conditions and could be removed as well.
Hope this helps!
- Badges
- Science topic
- Similar topics
- Biological Science
- Kinesiology
- Running