kr-colab / im_clam Goto Github PK

./im_clam -s stateSpaceFiles/testSS_33 -m stateSpaceFiles/testSS_33_mats -d exampleInput_33 -u 0.00001 -seqLen 100000
.___   _____             .__
|   | /     \       ____ |  | _____    _____
|   |/  \ /  \    _/ ___\|  | \__  \  /     \
|   /    Y    \   \  \___|  |__/ __ \|  Y Y  \
|___\____|__  /____\___  >____(____  /__|_|  /
            \/_____/   \/          \/      \/


im_clam -- Isolation with Migration Composite Likelihood Analysis using Markov chains


setup time:0.000000 secs


Parameter estimation run mode

now optimizing....

initial parameter guess:
5.249593 9.150234 11.652572 11.325893 5.873512 

nlopt failed with code -1!
for scaling:
u: 1.000000e-05 gen: 20.000000 N0:2484.745858 meanTreeLength:6.353064 seqLen:100000
Composite Likelihood estimates of IM params (scaled by 1/theta_pop1):
theta_pop2	theta_anc	mig_1->2	mig_2->1	t_div
1.013740	0.904389	0.100406	0.098707	2.103962	

Composite Likelihood estimates of IM params (unscaled):
theta_pop2	theta_anc	mig_1->2	mig_2->1	t_div
2518.887310	2247.175844	0.000010	0.000010	1045.562027	

Uncertainty estimates of IM params (scaled by 1/theta_pop1):
theta_pop2	theta_anc	mig_1->2	mig_2->1	t_div
0.017684	0.614021	0.003796	0.004887	0.355798	

Uncertainty estimates of IM params (unscaled):
theta_pop2	theta_anc	mig_1->2	mig_2->1	t_div
43.939041	1525.685700	0.000000	0.000000	176.813465	


likelihood: -125869.001381

Expected AFS:
0.000000 0.165561 0.089683 0.211452 
0.163769 0.004994 0.003885 0.017372 
0.089038 0.003865 0.002239 0.009124 
0.212445 0.017415 0.009159 0.000000 
///////////
total run time:68.260000 secs
 Liklihood Func. Evals: 1974