When I run the test data, I get this error (warning?), but still seem to get the correct answer. Should I be worried? I am using nlopt-2.4.2.
./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