Note
Go to the end to download the full example code.
ASW inference - Three pulses model#
This example implements inference for the ASW population under a three pulses model of admixture, using the tracts package. Inference is performed using autosomal and X chromosome data, allowing for the specification of sex-biased admixture.
To implement this example, we use the following driver file:
samples:
directory: ../ASW_TrioPhased
individual_names: [
"NA19625","NA19700","NA19701","NA19703","NA19704","NA19707","NA19711","NA19712","NA19713","NA19818","NA19819",
"NA19834","NA19835","NA19900","NA19901","NA19904","NA19908","NA19909","NA19913","NA19914","NA19916","NA19917",
"NA19920","NA19921","NA19922","NA19923","NA19982","NA19984","NA20126","NA20127","NA20274","NA20276","NA20278",
"NA20281","NA20282","NA20287","NA20289","NA20291","NA20294","NA20296","NA20298","NA20299","NA20314","NA20317",
"NA20318","NA20320","NA20321","NA20332","NA20334","NA20339","NA20340","NA20342","NA20346","NA20348","NA20351",
"NA20355","NA20356","NA20357","NA20359","NA20362","NA20412"]
male_names : [
"NA19700","NA19703","NA19711","NA19818","NA19834","NA19900","NA19904","NA19908","NA19916","NA19920",
"NA19922","NA19982","NA19984","NA20126","NA20278","NA20281","NA20291","NA20298","NA20318","NA20340",
"NA20342","NA20346","NA20348","NA20351","NA20356","NA20362"] #see Readme_dataprocessing.md for how this was generated
filename_format: "{name}_{label}_final.bed"
labels: [A, B] #If this field is omitted, 'A' and 'B' will be used by default
chromosomes: 1-22 #The chromosomes to use for analysis. Can be specified as a list or a range
allosomes: [X]
models:
model_filename: ../models/ppx_xxp_pxx.yaml
ad_model_autosomes: M
ad_model_allosomes: DC
start_params:
t1: 9:12
REUR: 0.8
RAFR: 0.9
REUR2: 0.2
t2: 5:8
t3: 1:4
REUR_sex_bias: -0.2:0.2
REUR2_sex_bias: -0.2:0.2
RAFR_sex_bias: -0.2:0.2
optim:
repetitions: 5
seed: 100
maximum_iterations: 1000
npts: 50
exclude_tracts_below_cm: 2
unknown_labels_for_smoothing: ["UNK", "centromere","miscall"] # segments with these labels will be smoother over, that is, will be filled with neighbouring ancestries up to their midpoints.
output:
output_directory: ./output_three_pulses/
output_filename_format: "ASW_test_output_{label}"
log_filename: 'ASW_three_pulses.log'
verbose_log: 1
verbose_screen: 30
log_scale: True
Complete results from this analysis are saved in the output directory specified in the driver file. Below, we display the optimal parameters estimated from this analysis, as well as the plots illustrating the inferred tract length distributions, compared to the observed histograms, for every source population and chromosome type (autosomes and X chromosome).
Optimal parameters#
parameter |
value |
|---|---|
REUR |
0.7594493513864244 |
REUR_sex_bias |
-0.999 |
t1 |
8.20018058951109 |
RAFR |
0.8832892626465892 |
RAFR_sex_bias |
-0.999 |
t2 |
6.963541367195792 |
REUR2 |
0.14602038131657813 |
REUR2_sex_bias |
-0.999 |
t3 |
5.993208977682569 |
targetpop_NAT_rate |
0.24043037328926875 |
targetpop_NAT_sex_bias |
1.0 |
likelihood -961.058 |
Optimal migration matrices#
Tract length histograms#
Autosomal admixture#
X chromosome admixture in females#
X chromosome admixture in males#
------------------------------------------------------------------------------------------------
Running tracts 2.0 with driver file: /home/runner/work/tracts/tracts/example/documentation_examples/ASW/ASW_three_pulses.yaml
------------------------------------------------------------------------------------------------
Results will be written to: output_three_pulses.
Using log file: output_three_pulses/ASW_three_pulses.log.
excluding_tracts_below set to 2.0 cM.
Re-optimization will be performed until convergence or maximum 5 times.
Ancestries: EUR, NAT, AFR
Data autosome proportions: [0.19578862 0.03825495 0.76595643]
Data allosome proportions: [0.16839124 0.03818939 0.79341937]
Model parameters and bounds:
--------------------------------------------
Parameter | Lower bound | Upper bound
--------------------------------------------
REUR | 1e-09 | 1
REUR_sex_bias | -1 | 1
t1 | 1 | inf
RAFR | 1e-09 | 1
RAFR_sex_bias | -1 | 1
t2 | 1 | inf
REUR2 | 1e-09 | 1
REUR2_sex_bias | -1 | 1
t3 | 1 | inf
--------------------------------------------
Multiple starting parameters will be generated and used for multiple optimization runs.
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Starting parameters for step 1 optimization
---------------------------------------------------------------------------------------------
Run | REUR | t1 | RAFR | t2 | REUR2 | t3
---------------------------------------------------------------------------------------------
1 | 0.8 | 9.867 | 0.9 | 6.789 | 0.2 | 3.064
2 | 0.8 | 9.854 | 0.9 | 6.8 | 0.2 | 2.506
3 | 0.8 | 11.99 | 0.9 | 5.966 | 0.2 | 3.074
4 | 0.8 | 9.395 | 0.9 | 5.824 | 0.2 | 3.504
5 | 0.8 | 10.64 | 0.9 | 6.023 | 0.2 | 2.971
---------------------------------------------------------------------------------------------
Optimization run #1
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -1278.6 , array([ 0.798945 , 0 , 9.8085 , 0.892502 , 0 , 6.92645 , 0.201858 , 0 , 3.55563 ]), Autosomes
60 , -937.233 , array([ 0.796271 , 0 , 9.83821 , 0.884552 , 0 , 7.40269 , 0.196754 , 0 , 4.10466 ]), Autosomes
90 , -789.812 , array([ 0.788992 , 0 , 9.17757 , 0.87948 , 0 , 7.67343 , 0.173544 , 0 , 4.65663 ]), Autosomes
120 , -730.679 , array([ 0.782085 , 0 , 8.86059 , 0.88044 , 0 , 7.70366 , 0.160901 , 0 , 4.90585 ]), Autosomes
150 , -720.888 , array([ 0.781072 , 0 , 8.86113 , 0.880239 , 0 , 7.70523 , 0.157726 , 0 , 4.96008 ]), Autosomes
180 , -714.705 , array([ 0.77997 , 0 , 8.85467 , 0.879915 , 0 , 7.70596 , 0.155684 , 0 , 4.98644 ]), Autosomes
206 , -711.155 , array([ 0.779737 , 0 , 8.85345 , 0.879825 , 0 , 7.7066 , 0.155182 , 0 , 4.99359 ]), Autosomes
Optimization completed.
-----------------------
Optimization run #2
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -1448.67 , array([ 0.79862 , 0 , 10.0701 , 0.892076 , 0 , 7.28883 , 0.20148 , 0 , 2.76194 ]), Autosomes
60 , -1028.04 , array([ 0.79667 , 0 , 9.91729 , 0.872188 , 0 , 7.9864 , 0.192084 , 0 , 3.15195 ]), Autosomes
90 , 1.49127e+30 , array([ 0.795679 , 0 , 9.72853 , 0.868025 , 0 , 8.44057 , 0.1892 , 0 , 3.28867 ]), OOB (oob=-0.01491266033413119)
120 , -879.829 , array([ 0.79575 , 0 , 9.74611 , 0.867777 , 0 , 8.425 , 0.188943 , 0 , 3.29491 ]), Autosomes
145 , -879.472 , array([ 0.795754 , 0 , 9.74615 , 0.867794 , 0 , 8.42464 , 0.188848 , 0 , 3.29576 ]), Autosomes
Optimization completed.
-----------------------
Optimization run #3
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -1344.44 , array([ 0.797258 , 0 , 12.0966 , 0.893864 , 0 , 6.61638 , 0.198304 , 0 , 3.49862 ]), Autosomes
60 , -1067.31 , array([ 0.792484 , 0 , 11.9045 , 0.883506 , 0 , 6.82478 , 0.192461 , 0 , 4.09496 ]), Autosomes
90 , -854.959 , array([ 0.785365 , 0 , 10.5917 , 0.87842 , 0 , 6.98906 , 0.176887 , 0 , 4.91177 ]), Autosomes
120 , -757.413 , array([ 0.782 , 0 , 9.7073 , 0.880135 , 0 , 7.38696 , 0.160937 , 0 , 5.00078 ]), Autosomes
150 , -696.14 , array([ 0.778346 , 0 , 8.95998 , 0.876457 , 0 , 7.65304 , 0.14431 , 0 , 5.00881 ]), Autosomes
180 , -679.413 , array([ 0.777402 , 0 , 8.8325 , 0.876516 , 0 , 7.65229 , 0.142612 , 0 , 5.08528 ]), Autosomes
208 , -679.258 , array([ 0.777219 , 0 , 8.83431 , 0.876593 , 0 , 7.65991 , 0.142664 , 0 , 5.08901 ]), Autosomes
Optimization completed.
-----------------------
Optimization run #4
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -1208.41 , array([ 0.796431 , 0 , 9.74832 , 0.894172 , 0 , 6.06708 , 0.202341 , 0 , 4.10765 ]), Autosomes
60 , -905.931 , array([ 0.790579 , 0 , 9.67757 , 0.887345 , 0 , 6.55233 , 0.195692 , 0 , 4.81247 ]), Autosomes
90 , -775.659 , array([ 0.780459 , 0 , 9.20241 , 0.885858 , 0 , 6.73737 , 0.176541 , 0 , 5.5511 ]), Autosomes
120 , -677.526 , array([ 0.764803 , 0 , 8.40064 , 0.883027 , 0 , 6.87889 , 0.149179 , 0 , 6.05145 ]), Autosomes
150 , -668.039 , array([ 0.761553 , 0 , 8.2793 , 0.883065 , 0 , 6.8886 , 0.146112 , 0 , 5.99589 ]), Autosomes
166 , -667.613 , array([ 0.761489 , 0 , 8.28143 , 0.883072 , 0 , 6.88913 , 0.146071 , 0 , 5.99373 ]), Autosomes
Optimization completed.
-----------------------
Optimization run #5
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -1345.02 , array([ 0.79753 , 0 , 10.806 , 0.892393 , 0 , 6.63758 , 0.197754 , 0 , 3.41564 ]), Autosomes
60 , -981.378 , array([ 0.792196 , 0 , 10.4417 , 0.877369 , 0 , 6.99194 , 0.191882 , 0 , 4.17112 ]), Autosomes
90 , -813.59 , array([ 0.786918 , 0 , 9.90972 , 0.878202 , 0 , 7.45719 , 0.175261 , 0 , 4.61679 ]), Autosomes
120 , -711.575 , array([ 0.7782 , 0 , 8.9548 , 0.87847 , 0 , 7.49016 , 0.150685 , 0 , 4.96187 ]), Autosomes
150 , -696.175 , array([ 0.778181 , 0 , 8.80015 , 0.877417 , 0 , 7.50103 , 0.148487 , 0 , 4.99452 ]), Autosomes
160 , -696.055 , array([ 0.778205 , 0 , 8.80218 , 0.877417 , 0 , 7.50035 , 0.148483 , 0 , 4.99397 ]), Autosomes
Optimization completed.
-----------------------
In Step 1: Results from multiple optimization runs with different starting parameters:
-------------------------------------------------------------------------------------------------------------------------------------------------------------
Run | LogLik | REUR | REUR_sex_bias | t1 | RAFR | RAFR_sex_bias | t2 | REUR2 | REUR2_sex_bias | t3
-------------------------------------------------------------------------------------------------------------------------------------------------------------
1 | -711.155 | 0.7797 | 0 | 8.853 | 0.8798 | 0 | 7.707 | 0.1552 | 0 | 4.994
2 | -879.472 | 0.7958 | 0 | 9.746 | 0.8678 | 0 | 8.425 | 0.1888 | 0 | 3.296
3 | -679.258 | 0.7772 | 0 | 8.834 | 0.8766 | 0 | 7.66 | 0.1427 | 0 | 5.089
4 | -667.613 | 0.7615 | 0 | 8.281 | 0.8831 | 0 | 6.889 | 0.1461 | 0 | 5.994
5 | -696.055 | 0.7782 | 0 | 8.802 | 0.8774 | 0 | 7.5 | 0.1485 | 0 | 4.994
-------------------------------------------------------------------------------------------------------------------------------------------------------------
Selecting best parameters from step 1 and proceeding to step 2 optimization.
----------------------------------------------------------------------------------------------------------------
Step 2 : Optimizing allosomal likelihood over parameters : ['REUR_sex_bias', 'RAFR_sex_bias', 'REUR2_sex_bias'].
----------------------------------------------------------------------------------------------------------------
Starting parameters for step 2 optimization (non-sex-bias parameters are fixed to the best step 1 estimates).
----------------------------------------------------------------------------------------------------------------------------------------------
Run | REUR | REUR_sex_bias | t1 | RAFR | RAFR_sex_bias | t2 | REUR2 | REUR2_sex_bias | t3
----------------------------------------------------------------------------------------------------------------------------------------------
1 | 0.7615 | 0.03862 | 8.281 | 0.8831 | 0.1895 | 6.889 | 0.1461 | 0.1641 | 5.994
2 | 0.7615 | 0.1926 | 8.281 | 0.8831 | 0.03241 | 6.889 | 0.1461 | 0.1983 | 5.994
3 | 0.7615 | -0.002335 | 8.281 | 0.8831 | -0.04257 | 6.889 | 0.1461 | 0.1197 | 5.994
4 | 0.7615 | -0.04409 | 8.281 | 0.8831 | -0.167 | 6.889 | 0.1461 | -0.1941 | 5.994
5 | 0.7615 | 0.009648 | 8.281 | 0.8831 | 0.1824 | 6.889 | 0.1461 | 0.147 | 5.994
----------------------------------------------------------------------------------------------------------------------------------------------
Optimization run #1
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -216.069 , array([ 0.761489 , -0.00979548 , 8.28143 , 0.883072 , 0.0726866 , 6.88913 , 0.146071 , 0.182862 , 5.99373 ]), Female allosomes
30 , -96.005 , array([ 0.761489 , -0.00979548 , 8.28143 , 0.883072 , 0.0726866 , 6.88913 , 0.146071 , 0.182862 , 5.99373 ]), Male allosomes
60 , -213.58 , array([ 0.761489 , -0.0624032 , 8.28143 , 0.883072 , -0.0654063 , 6.88913 , 0.146071 , 0.192606 , 5.99373 ]), Female allosomes
60 , -96.354 , array([ 0.761489 , -0.0624032 , 8.28143 , 0.883072 , -0.0654063 , 6.88913 , 0.146071 , 0.192606 , 5.99373 ]), Male allosomes
90 , -211.358 , array([ 0.761489 , -0.11614 , 8.28143 , 0.883072 , -0.199499 , 6.88913 , 0.146071 , 0.198174 , 5.99373 ]), Female allosomes
90 , -96.7154 , array([ 0.761489 , -0.11614 , 8.28143 , 0.883072 , -0.199499 , 6.88913 , 0.146071 , 0.198174 , 5.99373 ]), Male allosomes
120 , -209.405 , array([ 0.761489 , -0.163753 , 8.28143 , 0.883072 , -0.3292 , 6.88913 , 0.146071 , 0.196415 , 5.99373 ]), Female allosomes
120 , -97.0704 , array([ 0.761489 , -0.163753 , 8.28143 , 0.883072 , -0.3292 , 6.88913 , 0.146071 , 0.196415 , 5.99373 ]), Male allosomes
150 , -207.728 , array([ 0.761489 , -0.216736 , 8.28143 , 0.883072 , -0.445356 , 6.88913 , 0.146071 , 0.195417 , 5.99373 ]), Female allosomes
150 , -97.4108 , array([ 0.761489 , -0.216736 , 8.28143 , 0.883072 , -0.445356 , 6.88913 , 0.146071 , 0.195417 , 5.99373 ]), Male allosomes
180 , -206.341 , array([ 0.761489 , -0.271058 , 8.28143 , 0.883072 , -0.545654 , 6.88913 , 0.146071 , 0.189117 , 5.99373 ]), Female allosomes
180 , -97.7047 , array([ 0.761489 , -0.271058 , 8.28143 , 0.883072 , -0.545654 , 6.88913 , 0.146071 , 0.189117 , 5.99373 ]), Male allosomes
210 , -205.223 , array([ 0.761489 , -0.324216 , 8.28143 , 0.883072 , -0.631624 , 6.88913 , 0.146071 , 0.157265 , 5.99373 ]), Female allosomes
210 , -97.89 , array([ 0.761489 , -0.324216 , 8.28143 , 0.883072 , -0.631624 , 6.88913 , 0.146071 , 0.157265 , 5.99373 ]), Male allosomes
240 , -204.34 , array([ 0.761489 , -0.381459 , 8.28143 , 0.883072 , -0.699847 , 6.88913 , 0.146071 , 0.107806 , 5.99373 ]), Female allosomes
240 , -97.9741 , array([ 0.761489 , -0.381459 , 8.28143 , 0.883072 , -0.699847 , 6.88913 , 0.146071 , 0.107806 , 5.99373 ]), Male allosomes
270 , -203.656 , array([ 0.761489 , -0.433527 , 8.28143 , 0.883072 , -0.752237 , 6.88913 , 0.146071 , 0.0525431 , 5.99373 ]), Female allosomes
270 , -97.9633 , array([ 0.761489 , -0.433527 , 8.28143 , 0.883072 , -0.752237 , 6.88913 , 0.146071 , 0.0525431 , 5.99373 ]), Male allosomes
300 , -203.062 , array([ 0.761489 , -0.487273 , 8.28143 , 0.883072 , -0.797776 , 6.88913 , 0.146071 , -0.0132023 , 5.99373 ]), Female allosomes
300 , -97.9329 , array([ 0.761489 , -0.487273 , 8.28143 , 0.883072 , -0.797776 , 6.88913 , 0.146071 , -0.0132023 , 5.99373 ]), Male allosomes
330 , -202.593 , array([ 0.761489 , -0.539398 , 8.28143 , 0.883072 , -0.832782 , 6.88913 , 0.146071 , -0.0907237 , 5.99373 ]), Female allosomes
330 , -97.8389 , array([ 0.761489 , -0.539398 , 8.28143 , 0.883072 , -0.832782 , 6.88913 , 0.146071 , -0.0907237 , 5.99373 ]), Male allosomes
360 , -202.186 , array([ 0.761489 , -0.589286 , 8.28143 , 0.883072 , -0.860925 , 6.88913 , 0.146071 , -0.16818 , 5.99373 ]), Female allosomes
360 , -97.7435 , array([ 0.761489 , -0.589286 , 8.28143 , 0.883072 , -0.860925 , 6.88913 , 0.146071 , -0.16818 , 5.99373 ]), Male allosomes
390 , -201.865 , array([ 0.761489 , -0.631992 , 8.28143 , 0.883072 , -0.884009 , 6.88913 , 0.146071 , -0.254057 , 5.99373 ]), Female allosomes
390 , -97.5966 , array([ 0.761489 , -0.631992 , 8.28143 , 0.883072 , -0.884009 , 6.88913 , 0.146071 , -0.254057 , 5.99373 ]), Male allosomes
420 , -201.619 , array([ 0.761489 , -0.670873 , 8.28143 , 0.883072 , -0.901055 , 6.88913 , 0.146071 , -0.346009 , 5.99373 ]), Female allosomes
420 , -97.4215 , array([ 0.761489 , -0.670873 , 8.28143 , 0.883072 , -0.901055 , 6.88913 , 0.146071 , -0.346009 , 5.99373 ]), Male allosomes
450 , -201.392 , array([ 0.761489 , -0.705726 , 8.28143 , 0.883072 , -0.916567 , 6.88913 , 0.146071 , -0.429779 , 5.99373 ]), Female allosomes
450 , -97.2673 , array([ 0.761489 , -0.705726 , 8.28143 , 0.883072 , -0.916567 , 6.88913 , 0.146071 , -0.429779 , 5.99373 ]), Male allosomes
480 , -201.188 , array([ 0.761489 , -0.739325 , 8.28143 , 0.883072 , -0.929625 , 6.88913 , 0.146071 , -0.505002 , 5.99373 ]), Female allosomes
480 , -97.1322 , array([ 0.761489 , -0.739325 , 8.28143 , 0.883072 , -0.929625 , 6.88913 , 0.146071 , -0.505002 , 5.99373 ]), Male allosomes
510 , -201.015 , array([ 0.761489 , -0.769161 , 8.28143 , 0.883072 , -0.940325 , 6.88913 , 0.146071 , -0.575796 , 5.99373 ]), Female allosomes
510 , -97.0046 , array([ 0.761489 , -0.769161 , 8.28143 , 0.883072 , -0.940325 , 6.88913 , 0.146071 , -0.575796 , 5.99373 ]), Male allosomes
540 , -200.863 , array([ 0.761489 , -0.795658 , 8.28143 , 0.883072 , -0.949528 , 6.88913 , 0.146071 , -0.63839 , 5.99373 ]), Female allosomes
540 , -96.8946 , array([ 0.761489 , -0.795658 , 8.28143 , 0.883072 , -0.949528 , 6.88913 , 0.146071 , -0.63839 , 5.99373 ]), Male allosomes
570 , -200.725 , array([ 0.761489 , -0.821169 , 8.28143 , 0.883072 , -0.957293 , 6.88913 , 0.146071 , -0.691935 , 5.99373 ]), Female allosomes
570 , -96.8041 , array([ 0.761489 , -0.821169 , 8.28143 , 0.883072 , -0.957293 , 6.88913 , 0.146071 , -0.691935 , 5.99373 ]), Male allosomes
600 , -200.608 , array([ 0.761489 , -0.844097 , 8.28143 , 0.883072 , -0.963638 , 6.88913 , 0.146071 , -0.739852 , 5.99373 ]), Female allosomes
600 , -96.7238 , array([ 0.761489 , -0.844097 , 8.28143 , 0.883072 , -0.963638 , 6.88913 , 0.146071 , -0.739852 , 5.99373 ]), Male allosomes
630 , -200.5 , array([ 0.761489 , -0.86474 , 8.28143 , 0.883072 , -0.96931 , 6.88913 , 0.146071 , -0.779294 , 5.99373 ]), Female allosomes
630 , -96.661 , array([ 0.761489 , -0.86474 , 8.28143 , 0.883072 , -0.96931 , 6.88913 , 0.146071 , -0.779294 , 5.99373 ]), Male allosomes
660 , -200.409 , array([ 0.761489 , -0.882181 , 8.28143 , 0.883072 , -0.974125 , 6.88913 , 0.146071 , -0.813688 , 5.99373 ]), Female allosomes
660 , -96.6067 , array([ 0.761489 , -0.882181 , 8.28143 , 0.883072 , -0.974125 , 6.88913 , 0.146071 , -0.813688 , 5.99373 ]), Male allosomes
690 , -200.33 , array([ 0.761489 , -0.897349 , 8.28143 , 0.883072 , -0.978219 , 6.88913 , 0.146071 , -0.842857 , 5.99373 ]), Female allosomes
690 , -96.5617 , array([ 0.761489 , -0.897349 , 8.28143 , 0.883072 , -0.978219 , 6.88913 , 0.146071 , -0.842857 , 5.99373 ]), Male allosomes
720 , -200.26 , array([ 0.761489 , -0.91161 , 8.28143 , 0.883072 , -0.981609 , 6.88913 , 0.146071 , -0.867797 , 5.99373 ]), Female allosomes
720 , -96.5241 , array([ 0.761489 , -0.91161 , 8.28143 , 0.883072 , -0.981609 , 6.88913 , 0.146071 , -0.867797 , 5.99373 ]), Male allosomes
750 , -200.201 , array([ 0.761489 , -0.923751 , 8.28143 , 0.883072 , -0.984469 , 6.88913 , 0.146071 , -0.889251 , 5.99373 ]), Female allosomes
750 , -96.492 , array([ 0.761489 , -0.923751 , 8.28143 , 0.883072 , -0.984469 , 6.88913 , 0.146071 , -0.889251 , 5.99373 ]), Male allosomes
780 , -200.15 , array([ 0.761489 , -0.934609 , 8.28143 , 0.883072 , -0.98685 , 6.88913 , 0.146071 , -0.907002 , 5.99373 ]), Female allosomes
780 , -96.4662 , array([ 0.761489 , -0.934609 , 8.28143 , 0.883072 , -0.98685 , 6.88913 , 0.146071 , -0.907002 , 5.99373 ]), Male allosomes
810 , -200.105 , array([ 0.761489 , -0.944309 , 8.28143 , 0.883072 , -0.988858 , 6.88913 , 0.146071 , -0.921983 , 5.99373 ]), Female allosomes
810 , -96.4449 , array([ 0.761489 , -0.944309 , 8.28143 , 0.883072 , -0.988858 , 6.88913 , 0.146071 , -0.921983 , 5.99373 ]), Male allosomes
840 , -200.067 , array([ 0.761489 , -0.952547 , 8.28143 , 0.883072 , -0.990577 , 6.88913 , 0.146071 , -0.934415 , 5.99373 ]), Female allosomes
840 , -96.4275 , array([ 0.761489 , -0.952547 , 8.28143 , 0.883072 , -0.990577 , 6.88913 , 0.146071 , -0.934415 , 5.99373 ]), Male allosomes
870 , -200.036 , array([ 0.761489 , -0.959485 , 8.28143 , 0.883072 , -0.991991 , 6.88913 , 0.146071 , -0.945241 , 5.99373 ]), Female allosomes
870 , -96.4123 , array([ 0.761489 , -0.959485 , 8.28143 , 0.883072 , -0.991991 , 6.88913 , 0.146071 , -0.945241 , 5.99373 ]), Male allosomes
900 , -200.007 , array([ 0.761489 , -0.965634 , 8.28143 , 0.883072 , -0.993229 , 6.88913 , 0.146071 , -0.95385 , 5.99373 ]), Female allosomes
900 , -96.4007 , array([ 0.761489 , -0.965634 , 8.28143 , 0.883072 , -0.993229 , 6.88913 , 0.146071 , -0.95385 , 5.99373 ]), Male allosomes
930 , -199.983 , array([ 0.761489 , -0.970915 , 8.28143 , 0.883072 , -0.994299 , 6.88913 , 0.146071 , -0.960913 , 5.99373 ]), Female allosomes
930 , -96.3915 , array([ 0.761489 , -0.970915 , 8.28143 , 0.883072 , -0.994299 , 6.88913 , 0.146071 , -0.960913 , 5.99373 ]), Male allosomes
960 , -199.962 , array([ 0.761489 , -0.975507 , 8.28143 , 0.883072 , -0.995166 , 6.88913 , 0.146071 , -0.966984 , 5.99373 ]), Female allosomes
960 , -96.3835 , array([ 0.761489 , -0.975507 , 8.28143 , 0.883072 , -0.995166 , 6.88913 , 0.146071 , -0.966984 , 5.99373 ]), Male allosomes
990 , -199.946 , array([ 0.761489 , -0.979191 , 8.28143 , 0.883072 , -0.995908 , 6.88913 , 0.146071 , -0.972347 , 5.99373 ]), Female allosomes
990 , -96.3763 , array([ 0.761489 , -0.979191 , 8.28143 , 0.883072 , -0.995908 , 6.88913 , 0.146071 , -0.972347 , 5.99373 ]), Male allosomes
1000 , -199.941 , array([ 0.761489 , -0.980308 , 8.28143 , 0.883072 , -0.996128 , 6.88913 , 0.146071 , -0.973959 , 5.99373 ]), Female allosomes
1000 , -96.3741 , array([ 0.761489 , -0.980308 , 8.28143 , 0.883072 , -0.996128 , 6.88913 , 0.146071 , -0.973959 , 5.99373 ]), Male allosomes
Optimization completed.
-----------------------
Optimization run #2
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -214.27 , array([ 0.761489 , 0.148176 , 8.28143 , 0.883072 , -0.08749 , 6.88913 , 0.146071 , 0.211743 , 5.99373 ]), Female allosomes
30 , -96.3988 , array([ 0.761489 , 0.148176 , 8.28143 , 0.883072 , -0.08749 , 6.88913 , 0.146071 , 0.211743 , 5.99373 ]), Male allosomes
60 , -211.933 , array([ 0.761489 , 0.0997021 , 8.28143 , 0.883072 , -0.223548 , 6.88913 , 0.146071 , 0.218962 , 5.99373 ]), Female allosomes
60 , -96.7678 , array([ 0.761489 , 0.0997021 , 8.28143 , 0.883072 , -0.223548 , 6.88913 , 0.146071 , 0.218962 , 5.99373 ]), Male allosomes
90 , -209.924 , array([ 0.761489 , 0.0479259 , 8.28143 , 0.883072 , -0.350179 , 6.88913 , 0.146071 , 0.216962 , 5.99373 ]), Female allosomes
90 , -97.1168 , array([ 0.761489 , 0.0479259 , 8.28143 , 0.883072 , -0.350179 , 6.88913 , 0.146071 , 0.216962 , 5.99373 ]), Male allosomes
120 , -208.21 , array([ 0.761489 , -0.00551056 , 8.28143 , 0.883072 , -0.464947 , 6.88913 , 0.146071 , 0.21296 , 5.99373 ]), Female allosomes
120 , -97.4474 , array([ 0.761489 , -0.00551056 , 8.28143 , 0.883072 , -0.464947 , 6.88913 , 0.146071 , 0.21296 , 5.99373 ]), Male allosomes
150 , -206.811 , array([ 0.761489 , -0.063352 , 8.28143 , 0.883072 , -0.5633 , 6.88913 , 0.146071 , 0.194717 , 5.99373 ]), Female allosomes
150 , -97.7042 , array([ 0.761489 , -0.063352 , 8.28143 , 0.883072 , -0.5633 , 6.88913 , 0.146071 , 0.194717 , 5.99373 ]), Male allosomes
180 , -205.678 , array([ 0.761489 , -0.126496 , 8.28143 , 0.883072 , -0.644901 , 6.88913 , 0.146071 , 0.161625 , 5.99373 ]), Female allosomes
180 , -97.8756 , array([ 0.761489 , -0.126496 , 8.28143 , 0.883072 , -0.644901 , 6.88913 , 0.146071 , 0.161625 , 5.99373 ]), Male allosomes
210 , -204.759 , array([ 0.761489 , -0.195654 , 8.28143 , 0.883072 , -0.710612 , 6.88913 , 0.146071 , 0.117227 , 5.99373 ]), Female allosomes
210 , -97.9695 , array([ 0.761489 , -0.195654 , 8.28143 , 0.883072 , -0.710612 , 6.88913 , 0.146071 , 0.117227 , 5.99373 ]), Male allosomes
240 , -204.066 , array([ 0.761489 , -0.262399 , 8.28143 , 0.883072 , -0.758588 , 6.88913 , 0.146071 , 0.0612442 , 5.99373 ]), Female allosomes
240 , -97.9461 , array([ 0.761489 , -0.262399 , 8.28143 , 0.883072 , -0.758588 , 6.88913 , 0.146071 , 0.0612442 , 5.99373 ]), Male allosomes
270 , -203.451 , array([ 0.761489 , -0.334432 , 8.28143 , 0.883072 , -0.800285 , 6.88913 , 0.146071 , -0.00368077 , 5.99373 ]), Female allosomes
270 , -97.9114 , array([ 0.761489 , -0.334432 , 8.28143 , 0.883072 , -0.800285 , 6.88913 , 0.146071 , -0.00368077 , 5.99373 ]), Male allosomes
300 , -202.939 , array([ 0.761489 , -0.403534 , 8.28143 , 0.883072 , -0.834034 , 6.88913 , 0.146071 , -0.0760559 , 5.99373 ]), Female allosomes
300 , -97.8324 , array([ 0.761489 , -0.403534 , 8.28143 , 0.883072 , -0.834034 , 6.88913 , 0.146071 , -0.0760559 , 5.99373 ]), Male allosomes
330 , -202.488 , array([ 0.761489 , -0.469274 , 8.28143 , 0.883072 , -0.861541 , 6.88913 , 0.146071 , -0.146645 , 5.99373 ]), Female allosomes
330 , -97.7576 , array([ 0.761489 , -0.469274 , 8.28143 , 0.883072 , -0.861541 , 6.88913 , 0.146071 , -0.146645 , 5.99373 ]), Male allosomes
360 , -202.14 , array([ 0.761489 , -0.52673 , 8.28143 , 0.883072 , -0.883412 , 6.88913 , 0.146071 , -0.232583 , 5.99373 ]), Female allosomes
360 , -97.6103 , array([ 0.761489 , -0.52673 , 8.28143 , 0.883072 , -0.883412 , 6.88913 , 0.146071 , -0.232583 , 5.99373 ]), Male allosomes
390 , -201.843 , array([ 0.761489 , -0.579774 , 8.28143 , 0.883072 , -0.901177 , 6.88913 , 0.146071 , -0.31747 , 5.99373 ]), Female allosomes
390 , -97.4581 , array([ 0.761489 , -0.579774 , 8.28143 , 0.883072 , -0.901177 , 6.88913 , 0.146071 , -0.31747 , 5.99373 ]), Male allosomes
420 , -201.592 , array([ 0.761489 , -0.624828 , 8.28143 , 0.883072 , -0.916354 , 6.88913 , 0.146071 , -0.400736 , 5.99373 ]), Female allosomes
420 , -97.3055 , array([ 0.761489 , -0.624828 , 8.28143 , 0.883072 , -0.916354 , 6.88913 , 0.146071 , -0.400736 , 5.99373 ]), Male allosomes
450 , -201.369 , array([ 0.761489 , -0.669135 , 8.28143 , 0.883072 , -0.928662 , 6.88913 , 0.146071 , -0.479339 , 5.99373 ]), Female allosomes
450 , -97.1622 , array([ 0.761489 , -0.669135 , 8.28143 , 0.883072 , -0.928662 , 6.88913 , 0.146071 , -0.479339 , 5.99373 ]), Male allosomes
480 , -201.173 , array([ 0.761489 , -0.710541 , 8.28143 , 0.883072 , -0.938811 , 6.88913 , 0.146071 , -0.550396 , 5.99373 ]), Female allosomes
480 , -97.0353 , array([ 0.761489 , -0.710541 , 8.28143 , 0.883072 , -0.938811 , 6.88913 , 0.146071 , -0.550396 , 5.99373 ]), Male allosomes
510 , -200.998 , array([ 0.761489 , -0.748064 , 8.28143 , 0.883072 , -0.947664 , 6.88913 , 0.146071 , -0.613903 , 5.99373 ]), Female allosomes
510 , -96.9254 , array([ 0.761489 , -0.748064 , 8.28143 , 0.883072 , -0.947664 , 6.88913 , 0.146071 , -0.613903 , 5.99373 ]), Male allosomes
540 , -200.846 , array([ 0.761489 , -0.779207 , 8.28143 , 0.883072 , -0.955457 , 6.88913 , 0.146071 , -0.67085 , 5.99373 ]), Female allosomes
540 , -96.8283 , array([ 0.761489 , -0.779207 , 8.28143 , 0.883072 , -0.955457 , 6.88913 , 0.146071 , -0.67085 , 5.99373 ]), Male allosomes
570 , -200.711 , array([ 0.761489 , -0.808663 , 8.28143 , 0.883072 , -0.961969 , 6.88913 , 0.146071 , -0.719934 , 5.99373 ]), Female allosomes
570 , -96.7477 , array([ 0.761489 , -0.808663 , 8.28143 , 0.883072 , -0.961969 , 6.88913 , 0.146071 , -0.719934 , 5.99373 ]), Male allosomes
600 , -200.596 , array([ 0.761489 , -0.834103 , 8.28143 , 0.883072 , -0.967391 , 6.88913 , 0.146071 , -0.764249 , 5.99373 ]), Female allosomes
600 , -96.6752 , array([ 0.761489 , -0.834103 , 8.28143 , 0.883072 , -0.967391 , 6.88913 , 0.146071 , -0.764249 , 5.99373 ]), Male allosomes
630 , -200.489 , array([ 0.761489 , -0.858144 , 8.28143 , 0.883072 , -0.972168 , 6.88913 , 0.146071 , -0.799682 , 5.99373 ]), Female allosomes
630 , -96.6211 , array([ 0.761489 , -0.858144 , 8.28143 , 0.883072 , -0.972168 , 6.88913 , 0.146071 , -0.799682 , 5.99373 ]), Male allosomes
660 , -200.396 , array([ 0.761489 , -0.877771 , 8.28143 , 0.883072 , -0.976543 , 6.88913 , 0.146071 , -0.829751 , 5.99373 ]), Female allosomes
660 , -96.5764 , array([ 0.761489 , -0.877771 , 8.28143 , 0.883072 , -0.976543 , 6.88913 , 0.146071 , -0.829751 , 5.99373 ]), Male allosomes
690 , -200.315 , array([ 0.761489 , -0.896534 , 8.28143 , 0.883072 , -0.97997 , 6.88913 , 0.146071 , -0.855346 , 5.99373 ]), Female allosomes
690 , -96.5393 , array([ 0.761489 , -0.896534 , 8.28143 , 0.883072 , -0.97997 , 6.88913 , 0.146071 , -0.855346 , 5.99373 ]), Male allosomes
720 , -200.246 , array([ 0.761489 , -0.911556 , 8.28143 , 0.883072 , -0.983047 , 6.88913 , 0.146071 , -0.877458 , 5.99373 ]), Female allosomes
720 , -96.5074 , array([ 0.761489 , -0.911556 , 8.28143 , 0.883072 , -0.983047 , 6.88913 , 0.146071 , -0.877458 , 5.99373 ]), Male allosomes
750 , -200.189 , array([ 0.761489 , -0.924543 , 8.28143 , 0.883072 , -0.98551 , 6.88913 , 0.146071 , -0.897387 , 5.99373 ]), Female allosomes
750 , -96.4781 , array([ 0.761489 , -0.924543 , 8.28143 , 0.883072 , -0.98551 , 6.88913 , 0.146071 , -0.897387 , 5.99373 ]), Male allosomes
780 , -200.14 , array([ 0.761489 , -0.935721 , 8.28143 , 0.883072 , -0.987637 , 6.88913 , 0.146071 , -0.914075 , 5.99373 ]), Female allosomes
780 , -96.4542 , array([ 0.761489 , -0.935721 , 8.28143 , 0.883072 , -0.987637 , 6.88913 , 0.146071 , -0.914075 , 5.99373 ]), Male allosomes
810 , -200.098 , array([ 0.761489 , -0.945187 , 8.28143 , 0.883072 , -0.989469 , 6.88913 , 0.146071 , -0.928074 , 5.99373 ]), Female allosomes
810 , -96.4345 , array([ 0.761489 , -0.945187 , 8.28143 , 0.883072 , -0.989469 , 6.88913 , 0.146071 , -0.928074 , 5.99373 ]), Male allosomes
840 , -200.06 , array([ 0.761489 , -0.953327 , 8.28143 , 0.883072 , -0.991105 , 6.88913 , 0.146071 , -0.939497 , 5.99373 ]), Female allosomes
840 , -96.4189 , array([ 0.761489 , -0.953327 , 8.28143 , 0.883072 , -0.991105 , 6.88913 , 0.146071 , -0.939497 , 5.99373 ]), Male allosomes
870 , -200.027 , array([ 0.761489 , -0.960527 , 8.28143 , 0.883072 , -0.992515 , 6.88913 , 0.146071 , -0.948486 , 5.99373 ]), Female allosomes
870 , -96.4074 , array([ 0.761489 , -0.960527 , 8.28143 , 0.883072 , -0.992515 , 6.88913 , 0.146071 , -0.948486 , 5.99373 ]), Male allosomes
900 , -200.001 , array([ 0.761489 , -0.966528 , 8.28143 , 0.883072 , -0.993628 , 6.88913 , 0.146071 , -0.956943 , 5.99373 ]), Female allosomes
900 , -96.3959 , array([ 0.761489 , -0.966528 , 8.28143 , 0.883072 , -0.993628 , 6.88913 , 0.146071 , -0.956943 , 5.99373 ]), Male allosomes
930 , -199.977 , array([ 0.761489 , -0.971661 , 8.28143 , 0.883072 , -0.994684 , 6.88913 , 0.146071 , -0.963291 , 5.99373 ]), Female allosomes
930 , -96.3879 , array([ 0.761489 , -0.971661 , 8.28143 , 0.883072 , -0.994684 , 6.88913 , 0.146071 , -0.963291 , 5.99373 ]), Male allosomes
960 , -199.957 , array([ 0.761489 , -0.976118 , 8.28143 , 0.883072 , -0.995501 , 6.88913 , 0.146071 , -0.968979 , 5.99373 ]), Female allosomes
960 , -96.3806 , array([ 0.761489 , -0.976118 , 8.28143 , 0.883072 , -0.995501 , 6.88913 , 0.146071 , -0.968979 , 5.99373 ]), Male allosomes
990 , -199.941 , array([ 0.761489 , -0.979814 , 8.28143 , 0.883072 , -0.996191 , 6.88913 , 0.146071 , -0.973898 , 5.99373 ]), Female allosomes
990 , -96.3741 , array([ 0.761489 , -0.979814 , 8.28143 , 0.883072 , -0.996191 , 6.88913 , 0.146071 , -0.973898 , 5.99373 ]), Male allosomes
1000 , -199.936 , array([ 0.761489 , -0.980884 , 8.28143 , 0.883072 , -0.996414 , 6.88913 , 0.146071 , -0.975307 , 5.99373 ]), Female allosomes
1000 , -96.3723 , array([ 0.761489 , -0.980884 , 8.28143 , 0.883072 , -0.996414 , 6.88913 , 0.146071 , -0.975307 , 5.99373 ]), Male allosomes
Optimization completed.
-----------------------
Optimization run #3
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -212.449 , array([ 0.761489 , -0.0375189 , 8.28143 , 0.883072 , -0.161583 , 6.88913 , 0.146071 , 0.132989 , 5.99373 ]), Female allosomes
30 , -96.4412 , array([ 0.761489 , -0.0375189 , 8.28143 , 0.883072 , -0.161583 , 6.88913 , 0.146071 , 0.132989 , 5.99373 ]), Male allosomes
60 , -210.343 , array([ 0.761489 , -0.0881341 , 8.28143 , 0.883072 , -0.293214 , 6.88913 , 0.146071 , 0.135683 , 5.99373 ]), Female allosomes
60 , -96.796 , array([ 0.761489 , -0.0881341 , 8.28143 , 0.883072 , -0.293214 , 6.88913 , 0.146071 , 0.135683 , 5.99373 ]), Male allosomes
90 , -208.548 , array([ 0.761489 , -0.137774 , 8.28143 , 0.883072 , -0.414332 , 6.88913 , 0.146071 , 0.127515 , 5.99373 ]), Female allosomes
90 , -97.1164 , array([ 0.761489 , -0.137774 , 8.28143 , 0.883072 , -0.414332 , 6.88913 , 0.146071 , 0.127515 , 5.99373 ]), Male allosomes
120 , -207.046 , array([ 0.761489 , -0.191719 , 8.28143 , 0.883072 , -0.520056 , 6.88913 , 0.146071 , 0.112984 , 5.99373 ]), Female allosomes
120 , -97.3925 , array([ 0.761489 , -0.191719 , 8.28143 , 0.883072 , -0.520056 , 6.88913 , 0.146071 , 0.112984 , 5.99373 ]), Male allosomes
150 , -205.797 , array([ 0.761489 , -0.243082 , 8.28143 , 0.883072 , -0.611703 , 6.88913 , 0.146071 , 0.0905672 , 5.99373 ]), Female allosomes
150 , -97.5844 , array([ 0.761489 , -0.243082 , 8.28143 , 0.883072 , -0.611703 , 6.88913 , 0.146071 , 0.0905672 , 5.99373 ]), Male allosomes
180 , -204.898 , array([ 0.761489 , -0.304037 , 8.28143 , 0.883072 , -0.676428 , 6.88913 , 0.146071 , 0.0529003 , 5.99373 ]), Female allosomes
180 , -97.6919 , array([ 0.761489 , -0.304037 , 8.28143 , 0.883072 , -0.676428 , 6.88913 , 0.146071 , 0.0529003 , 5.99373 ]), Male allosomes
210 , -204.064 , array([ 0.761489 , -0.358549 , 8.28143 , 0.883072 , -0.739409 , 6.88913 , 0.146071 , 0.00664666 , 5.99373 ]), Female allosomes
210 , -97.7685 , array([ 0.761489 , -0.358549 , 8.28143 , 0.883072 , -0.739409 , 6.88913 , 0.146071 , 0.00664666 , 5.99373 ]), Male allosomes
240 , -203.414 , array([ 0.761489 , -0.412205 , 8.28143 , 0.883072 , -0.788449 , 6.88913 , 0.146071 , -0.0553743 , 5.99373 ]), Female allosomes
240 , -97.7598 , array([ 0.761489 , -0.412205 , 8.28143 , 0.883072 , -0.788449 , 6.88913 , 0.146071 , -0.0553743 , 5.99373 ]), Male allosomes
270 , -202.897 , array([ 0.761489 , -0.470541 , 8.28143 , 0.883072 , -0.823882 , 6.88913 , 0.146071 , -0.127586 , 5.99373 ]), Female allosomes
270 , -97.702 , array([ 0.761489 , -0.470541 , 8.28143 , 0.883072 , -0.823882 , 6.88913 , 0.146071 , -0.127586 , 5.99373 ]), Male allosomes
300 , -202.459 , array([ 0.761489 , -0.522045 , 8.28143 , 0.883072 , -0.854945 , 6.88913 , 0.146071 , -0.204922 , 5.99373 ]), Female allosomes
300 , -97.6025 , array([ 0.761489 , -0.522045 , 8.28143 , 0.883072 , -0.854945 , 6.88913 , 0.146071 , -0.204922 , 5.99373 ]), Male allosomes
330 , -202.093 , array([ 0.761489 , -0.572472 , 8.28143 , 0.883072 , -0.879505 , 6.88913 , 0.146071 , -0.284098 , 5.99373 ]), Female allosomes
330 , -97.4827 , array([ 0.761489 , -0.572472 , 8.28143 , 0.883072 , -0.879505 , 6.88913 , 0.146071 , -0.284098 , 5.99373 ]), Male allosomes
360 , -201.774 , array([ 0.761489 , -0.620885 , 8.28143 , 0.883072 , -0.899721 , 6.88913 , 0.146071 , -0.359514 , 5.99373 ]), Female allosomes
360 , -97.3637 , array([ 0.761489 , -0.620885 , 8.28143 , 0.883072 , -0.899721 , 6.88913 , 0.146071 , -0.359514 , 5.99373 ]), Male allosomes
390 , -201.522 , array([ 0.761489 , -0.662131 , 8.28143 , 0.883072 , -0.915575 , 6.88913 , 0.146071 , -0.439086 , 5.99373 ]), Female allosomes
390 , -97.2241 , array([ 0.761489 , -0.662131 , 8.28143 , 0.883072 , -0.915575 , 6.88913 , 0.146071 , -0.439086 , 5.99373 ]), Male allosomes
420 , -201.301 , array([ 0.761489 , -0.700127 , 8.28143 , 0.883072 , -0.928893 , 6.88913 , 0.146071 , -0.512703 , 5.99373 ]), Female allosomes
420 , -97.0957 , array([ 0.761489 , -0.700127 , 8.28143 , 0.883072 , -0.928893 , 6.88913 , 0.146071 , -0.512703 , 5.99373 ]), Male allosomes
450 , -201.107 , array([ 0.761489 , -0.736039 , 8.28143 , 0.883072 , -0.939884 , 6.88913 , 0.146071 , -0.57871 , 5.99373 ]), Female allosomes
450 , -96.9827 , array([ 0.761489 , -0.736039 , 8.28143 , 0.883072 , -0.939884 , 6.88913 , 0.146071 , -0.57871 , 5.99373 ]), Male allosomes
480 , -200.931 , array([ 0.761489 , -0.769633 , 8.28143 , 0.883072 , -0.949358 , 6.88913 , 0.146071 , -0.635445 , 5.99373 ]), Female allosomes
480 , -96.8896 , array([ 0.761489 , -0.769633 , 8.28143 , 0.883072 , -0.949358 , 6.88913 , 0.146071 , -0.635445 , 5.99373 ]), Male allosomes
510 , -200.787 , array([ 0.761489 , -0.799432 , 8.28143 , 0.883072 , -0.956703 , 6.88913 , 0.146071 , -0.689895 , 5.99373 ]), Female allosomes
510 , -96.7978 , array([ 0.761489 , -0.799432 , 8.28143 , 0.883072 , -0.956703 , 6.88913 , 0.146071 , -0.689895 , 5.99373 ]), Male allosomes
540 , -200.657 , array([ 0.761489 , -0.825519 , 8.28143 , 0.883072 , -0.963331 , 6.88913 , 0.146071 , -0.736604 , 5.99373 ]), Female allosomes
540 , -96.7219 , array([ 0.761489 , -0.825519 , 8.28143 , 0.883072 , -0.963331 , 6.88913 , 0.146071 , -0.736604 , 5.99373 ]), Male allosomes
570 , -200.544 , array([ 0.761489 , -0.848647 , 8.28143 , 0.883072 , -0.968955 , 6.88913 , 0.146071 , -0.776253 , 5.99373 ]), Female allosomes
570 , -96.6595 , array([ 0.761489 , -0.848647 , 8.28143 , 0.883072 , -0.968955 , 6.88913 , 0.146071 , -0.776253 , 5.99373 ]), Male allosomes
600 , -200.444 , array([ 0.761489 , -0.869255 , 8.28143 , 0.883072 , -0.9738 , 6.88913 , 0.146071 , -0.81017 , 5.99373 ]), Female allosomes
600 , -96.6076 , array([ 0.761489 , -0.869255 , 8.28143 , 0.883072 , -0.9738 , 6.88913 , 0.146071 , -0.81017 , 5.99373 ]), Male allosomes
630 , -200.358 , array([ 0.761489 , -0.888111 , 8.28143 , 0.883072 , -0.977716 , 6.88913 , 0.146071 , -0.839579 , 5.99373 ]), Female allosomes
630 , -96.5632 , array([ 0.761489 , -0.888111 , 8.28143 , 0.883072 , -0.977716 , 6.88913 , 0.146071 , -0.839579 , 5.99373 ]), Male allosomes
660 , -200.284 , array([ 0.761489 , -0.905406 , 8.28143 , 0.883072 , -0.980907 , 6.88913 , 0.146071 , -0.864303 , 5.99373 ]), Female allosomes
660 , -96.527 , array([ 0.761489 , -0.905406 , 8.28143 , 0.883072 , -0.980907 , 6.88913 , 0.146071 , -0.864303 , 5.99373 ]), Male allosomes
690 , -200.219 , array([ 0.761489 , -0.919254 , 8.28143 , 0.883072 , -0.983876 , 6.88913 , 0.146071 , -0.884987 , 5.99373 ]), Female allosomes
690 , -96.4974 , array([ 0.761489 , -0.919254 , 8.28143 , 0.883072 , -0.983876 , 6.88913 , 0.146071 , -0.884987 , 5.99373 ]), Male allosomes
720 , -200.164 , array([ 0.761489 , -0.930839 , 8.28143 , 0.883072 , -0.986406 , 6.88913 , 0.146071 , -0.903001 , 5.99373 ]), Female allosomes
720 , -96.4716 , array([ 0.761489 , -0.930839 , 8.28143 , 0.883072 , -0.986406 , 6.88913 , 0.146071 , -0.903001 , 5.99373 ]), Male allosomes
750 , -200.117 , array([ 0.761489 , -0.941028 , 8.28143 , 0.883072 , -0.988492 , 6.88913 , 0.146071 , -0.918454 , 5.99373 ]), Female allosomes
750 , -96.4496 , array([ 0.761489 , -0.941028 , 8.28143 , 0.883072 , -0.988492 , 6.88913 , 0.146071 , -0.918454 , 5.99373 ]), Male allosomes
780 , -200.076 , array([ 0.761489 , -0.950058 , 8.28143 , 0.883072 , -0.990261 , 6.88913 , 0.146071 , -0.931116 , 5.99373 ]), Female allosomes
780 , -96.4322 , array([ 0.761489 , -0.950058 , 8.28143 , 0.883072 , -0.990261 , 6.88913 , 0.146071 , -0.931116 , 5.99373 ]), Male allosomes
810 , -200.042 , array([ 0.761489 , -0.957753 , 8.28143 , 0.883072 , -0.991768 , 6.88913 , 0.146071 , -0.941694 , 5.99373 ]), Female allosomes
810 , -96.418 , array([ 0.761489 , -0.957753 , 8.28143 , 0.883072 , -0.991768 , 6.88913 , 0.146071 , -0.941694 , 5.99373 ]), Male allosomes
840 , -200.013 , array([ 0.761489 , -0.964067 , 8.28143 , 0.883072 , -0.993014 , 6.88913 , 0.146071 , -0.951266 , 5.99373 ]), Female allosomes
840 , -96.4046 , array([ 0.761489 , -0.964067 , 8.28143 , 0.883072 , -0.993014 , 6.88913 , 0.146071 , -0.951266 , 5.99373 ]), Male allosomes
870 , -199.988 , array([ 0.761489 , -0.969546 , 8.28143 , 0.883072 , -0.994093 , 6.88913 , 0.146071 , -0.959008 , 5.99373 ]), Female allosomes
870 , -96.3941 , array([ 0.761489 , -0.969546 , 8.28143 , 0.883072 , -0.994093 , 6.88913 , 0.146071 , -0.959008 , 5.99373 ]), Male allosomes
900 , -199.967 , array([ 0.761489 , -0.974406 , 8.28143 , 0.883072 , -0.994999 , 6.88913 , 0.146071 , -0.965266 , 5.99373 ]), Female allosomes
900 , -96.386 , array([ 0.761489 , -0.974406 , 8.28143 , 0.883072 , -0.994999 , 6.88913 , 0.146071 , -0.965266 , 5.99373 ]), Male allosomes
930 , -199.949 , array([ 0.761489 , -0.978292 , 8.28143 , 0.883072 , -0.995785 , 6.88913 , 0.146071 , -0.970729 , 5.99373 ]), Female allosomes
930 , -96.3787 , array([ 0.761489 , -0.978292 , 8.28143 , 0.883072 , -0.995785 , 6.88913 , 0.146071 , -0.970729 , 5.99373 ]), Male allosomes
960 , -199.934 , array([ 0.761489 , -0.981514 , 8.28143 , 0.883072 , -0.996433 , 6.88913 , 0.146071 , -0.975489 , 5.99373 ]), Female allosomes
960 , -96.3723 , array([ 0.761489 , -0.981514 , 8.28143 , 0.883072 , -0.996433 , 6.88913 , 0.146071 , -0.975489 , 5.99373 ]), Male allosomes
990 , -199.921 , array([ 0.761489 , -0.984343 , 8.28143 , 0.883072 , -0.99702 , 6.88913 , 0.146071 , -0.978984 , 5.99373 ]), Female allosomes
990 , -96.368 , array([ 0.761489 , -0.984343 , 8.28143 , 0.883072 , -0.99702 , 6.88913 , 0.146071 , -0.978984 , 5.99373 ]), Male allosomes
1000 , -199.917 , array([ 0.761489 , -0.985203 , 8.28143 , 0.883072 , -0.997179 , 6.88913 , 0.146071 , -0.980189 , 5.99373 ]), Female allosomes
1000 , -96.3664 , array([ 0.761489 , -0.985203 , 8.28143 , 0.883072 , -0.997179 , 6.88913 , 0.146071 , -0.980189 , 5.99373 ]), Male allosomes
Optimization completed.
-----------------------
Optimization run #4
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -211.28 , array([ 0.761489 , -0.0827976 , 8.28143 , 0.883072 , -0.282535 , 6.88913 , 0.146071 , -0.181987 , 5.99373 ]), Female allosomes
30 , -96.04 , array([ 0.761489 , -0.0827976 , 8.28143 , 0.883072 , -0.282535 , 6.88913 , 0.146071 , -0.181987 , 5.99373 ]), Male allosomes
60 , -209.319 , array([ 0.761489 , -0.122425 , 8.28143 , 0.883072 , -0.4071 , 6.88913 , 0.146071 , -0.179445 , 5.99373 ]), Female allosomes
60 , -96.3474 , array([ 0.761489 , -0.122425 , 8.28143 , 0.883072 , -0.4071 , 6.88913 , 0.146071 , -0.179445 , 5.99373 ]), Male allosomes
90 , -207.647 , array([ 0.761489 , -0.16752 , 8.28143 , 0.883072 , -0.517582 , 6.88913 , 0.146071 , -0.181974 , 5.99373 ]), Female allosomes
90 , -96.6286 , array([ 0.761489 , -0.16752 , 8.28143 , 0.883072 , -0.517582 , 6.88913 , 0.146071 , -0.181974 , 5.99373 ]), Male allosomes
120 , -206.276 , array([ 0.761489 , -0.216859 , 8.28143 , 0.883072 , -0.611043 , 6.88913 , 0.146071 , -0.193755 , 5.99373 ]), Female allosomes
120 , -96.8575 , array([ 0.761489 , -0.216859 , 8.28143 , 0.883072 , -0.611043 , 6.88913 , 0.146071 , -0.193755 , 5.99373 ]), Male allosomes
150 , -205.177 , array([ 0.761489 , -0.267448 , 8.28143 , 0.883072 , -0.688263 , 6.88913 , 0.146071 , -0.222838 , 5.99373 ]), Female allosomes
150 , -97.0078 , array([ 0.761489 , -0.267448 , 8.28143 , 0.883072 , -0.688263 , 6.88913 , 0.146071 , -0.222838 , 5.99373 ]), Male allosomes
180 , -204.296 , array([ 0.761489 , -0.325207 , 8.28143 , 0.883072 , -0.748933 , 6.88913 , 0.146071 , -0.260563 , 5.99373 ]), Female allosomes
180 , -97.0971 , array([ 0.761489 , -0.325207 , 8.28143 , 0.883072 , -0.748933 , 6.88913 , 0.146071 , -0.260563 , 5.99373 ]), Male allosomes
210 , -203.599 , array([ 0.761489 , -0.379841 , 8.28143 , 0.883072 , -0.797063 , 6.88913 , 0.146071 , -0.309084 , 5.99373 ]), Female allosomes
210 , -97.125 , array([ 0.761489 , -0.379841 , 8.28143 , 0.883072 , -0.797063 , 6.88913 , 0.146071 , -0.309084 , 5.99373 ]), Male allosomes
240 , -203.006 , array([ 0.761489 , -0.437448 , 8.28143 , 0.883072 , -0.836296 , 6.88913 , 0.146071 , -0.358835 , 5.99373 ]), Female allosomes
240 , -97.1274 , array([ 0.761489 , -0.437448 , 8.28143 , 0.883072 , -0.836296 , 6.88913 , 0.146071 , -0.358835 , 5.99373 ]), Male allosomes
270 , -202.522 , array([ 0.761489 , -0.493605 , 8.28143 , 0.883072 , -0.867054 , 6.88913 , 0.146071 , -0.413488 , 5.99373 ]), Female allosomes
270 , -97.0954 , array([ 0.761489 , -0.493605 , 8.28143 , 0.883072 , -0.867054 , 6.88913 , 0.146071 , -0.413488 , 5.99373 ]), Male allosomes
300 , -202.12 , array([ 0.761489 , -0.550269 , 8.28143 , 0.883072 , -0.890867 , 6.88913 , 0.146071 , -0.469564 , 5.99373 ]), Female allosomes
300 , -97.0432 , array([ 0.761489 , -0.550269 , 8.28143 , 0.883072 , -0.890867 , 6.88913 , 0.146071 , -0.469564 , 5.99373 ]), Male allosomes
330 , -201.777 , array([ 0.761489 , -0.606176 , 8.28143 , 0.883072 , -0.90962 , 6.88913 , 0.146071 , -0.521638 , 5.99373 ]), Female allosomes
330 , -96.9889 , array([ 0.761489 , -0.606176 , 8.28143 , 0.883072 , -0.90962 , 6.88913 , 0.146071 , -0.521638 , 5.99373 ]), Male allosomes
360 , -201.504 , array([ 0.761489 , -0.655197 , 8.28143 , 0.883072 , -0.924108 , 6.88913 , 0.146071 , -0.57985 , 5.99373 ]), Female allosomes
360 , -96.9089 , array([ 0.761489 , -0.655197 , 8.28143 , 0.883072 , -0.924108 , 6.88913 , 0.146071 , -0.57985 , 5.99373 ]), Male allosomes
390 , -201.262 , array([ 0.761489 , -0.699992 , 8.28143 , 0.883072 , -0.936497 , 6.88913 , 0.146071 , -0.631259 , 5.99373 ]), Female allosomes
390 , -96.8399 , array([ 0.761489 , -0.699992 , 8.28143 , 0.883072 , -0.936497 , 6.88913 , 0.146071 , -0.631259 , 5.99373 ]), Male allosomes
420 , -201.058 , array([ 0.761489 , -0.738627 , 8.28143 , 0.883072 , -0.946785 , 6.88913 , 0.146071 , -0.679403 , 5.99373 ]), Female allosomes
420 , -96.7725 , array([ 0.761489 , -0.738627 , 8.28143 , 0.883072 , -0.946785 , 6.88913 , 0.146071 , -0.679403 , 5.99373 ]), Male allosomes
450 , -200.884 , array([ 0.761489 , -0.771314 , 8.28143 , 0.883072 , -0.955596 , 6.88913 , 0.146071 , -0.723003 , 5.99373 ]), Female allosomes
450 , -96.7108 , array([ 0.761489 , -0.771314 , 8.28143 , 0.883072 , -0.955596 , 6.88913 , 0.146071 , -0.723003 , 5.99373 ]), Male allosomes
480 , -200.732 , array([ 0.761489 , -0.803133 , 8.28143 , 0.883072 , -0.962583 , 6.88913 , 0.146071 , -0.761496 , 5.99373 ]), Female allosomes
480 , -96.6572 , array([ 0.761489 , -0.803133 , 8.28143 , 0.883072 , -0.962583 , 6.88913 , 0.146071 , -0.761496 , 5.99373 ]), Male allosomes
510 , -200.607 , array([ 0.761489 , -0.83034 , 8.28143 , 0.883072 , -0.968236 , 6.88913 , 0.146071 , -0.797126 , 5.99373 ]), Female allosomes
510 , -96.6059 , array([ 0.761489 , -0.83034 , 8.28143 , 0.883072 , -0.968236 , 6.88913 , 0.146071 , -0.797126 , 5.99373 ]), Male allosomes
540 , -200.493 , array([ 0.761489 , -0.856543 , 8.28143 , 0.883072 , -0.972924 , 6.88913 , 0.146071 , -0.825786 , 5.99373 ]), Female allosomes
540 , -96.5676 , array([ 0.761489 , -0.856543 , 8.28143 , 0.883072 , -0.972924 , 6.88913 , 0.146071 , -0.825786 , 5.99373 ]), Male allosomes
570 , -200.395 , array([ 0.761489 , -0.878416 , 8.28143 , 0.883072 , -0.977057 , 6.88913 , 0.146071 , -0.85083 , 5.99373 ]), Female allosomes
570 , -96.5344 , array([ 0.761489 , -0.878416 , 8.28143 , 0.883072 , -0.977057 , 6.88913 , 0.146071 , -0.85083 , 5.99373 ]), Male allosomes
600 , -200.312 , array([ 0.761489 , -0.896786 , 8.28143 , 0.883072 , -0.980571 , 6.88913 , 0.146071 , -0.873136 , 5.99373 ]), Female allosomes
600 , -96.5045 , array([ 0.761489 , -0.896786 , 8.28143 , 0.883072 , -0.980571 , 6.88913 , 0.146071 , -0.873136 , 5.99373 ]), Male allosomes
630 , -200.244 , array([ 0.761489 , -0.912008 , 8.28143 , 0.883072 , -0.983538 , 6.88913 , 0.146071 , -0.89278 , 5.99373 ]), Female allosomes
630 , -96.4778 , array([ 0.761489 , -0.912008 , 8.28143 , 0.883072 , -0.983538 , 6.88913 , 0.146071 , -0.89278 , 5.99373 ]), Male allosomes
660 , -200.182 , array([ 0.761489 , -0.925926 , 8.28143 , 0.883072 , -0.986078 , 6.88913 , 0.146071 , -0.908247 , 5.99373 ]), Female allosomes
660 , -96.4582 , array([ 0.761489 , -0.925926 , 8.28143 , 0.883072 , -0.986078 , 6.88913 , 0.146071 , -0.908247 , 5.99373 ]), Male allosomes
690 , -200.131 , array([ 0.761489 , -0.937379 , 8.28143 , 0.883072 , -0.98819 , 6.88913 , 0.146071 , -0.922127 , 5.99373 ]), Female allosomes
690 , -96.4401 , array([ 0.761489 , -0.937379 , 8.28143 , 0.883072 , -0.98819 , 6.88913 , 0.146071 , -0.922127 , 5.99373 ]), Male allosomes
720 , -200.088 , array([ 0.761489 , -0.947054 , 8.28143 , 0.883072 , -0.99003 , 6.88913 , 0.146071 , -0.933679 , 5.99373 ]), Female allosomes
720 , -96.4254 , array([ 0.761489 , -0.947054 , 8.28143 , 0.883072 , -0.99003 , 6.88913 , 0.146071 , -0.933679 , 5.99373 ]), Male allosomes
750 , -200.051 , array([ 0.761489 , -0.955288 , 8.28143 , 0.883072 , -0.991587 , 6.88913 , 0.146071 , -0.94395 , 5.99373 ]), Female allosomes
750 , -96.4121 , array([ 0.761489 , -0.955288 , 8.28143 , 0.883072 , -0.991587 , 6.88913 , 0.146071 , -0.94395 , 5.99373 ]), Male allosomes
780 , -200.019 , array([ 0.761489 , -0.96229 , 8.28143 , 0.883072 , -0.992905 , 6.88913 , 0.146071 , -0.952627 , 5.99373 ]), Female allosomes
780 , -96.4009 , array([ 0.761489 , -0.96229 , 8.28143 , 0.883072 , -0.992905 , 6.88913 , 0.146071 , -0.952627 , 5.99373 ]), Male allosomes
810 , -199.993 , array([ 0.761489 , -0.968189 , 8.28143 , 0.883072 , -0.994019 , 6.88913 , 0.146071 , -0.960001 , 5.99373 ]), Female allosomes
810 , -96.3914 , array([ 0.761489 , -0.968189 , 8.28143 , 0.883072 , -0.994019 , 6.88913 , 0.146071 , -0.960001 , 5.99373 ]), Male allosomes
840 , -199.971 , array([ 0.761489 , -0.973181 , 8.28143 , 0.883072 , -0.994955 , 6.88913 , 0.146071 , -0.966266 , 5.99373 ]), Female allosomes
840 , -96.3834 , array([ 0.761489 , -0.973181 , 8.28143 , 0.883072 , -0.994955 , 6.88913 , 0.146071 , -0.966266 , 5.99373 ]), Male allosomes
870 , -199.952 , array([ 0.761489 , -0.977416 , 8.28143 , 0.883072 , -0.995744 , 6.88913 , 0.146071 , -0.971548 , 5.99373 ]), Female allosomes
870 , -96.3767 , array([ 0.761489 , -0.977416 , 8.28143 , 0.883072 , -0.995744 , 6.88913 , 0.146071 , -0.971548 , 5.99373 ]), Male allosomes
900 , -199.936 , array([ 0.761489 , -0.980987 , 8.28143 , 0.883072 , -0.996413 , 6.88913 , 0.146071 , -0.975996 , 5.99373 ]), Female allosomes
900 , -96.371 , array([ 0.761489 , -0.980987 , 8.28143 , 0.883072 , -0.996413 , 6.88913 , 0.146071 , -0.975996 , 5.99373 ]), Male allosomes
930 , -199.922 , array([ 0.761489 , -0.983985 , 8.28143 , 0.883072 , -0.996978 , 6.88913 , 0.146071 , -0.979762 , 5.99373 ]), Female allosomes
930 , -96.3663 , array([ 0.761489 , -0.983985 , 8.28143 , 0.883072 , -0.996978 , 6.88913 , 0.146071 , -0.979762 , 5.99373 ]), Male allosomes
960 , -199.911 , array([ 0.761489 , -0.986512 , 8.28143 , 0.883072 , -0.997453 , 6.88913 , 0.146071 , -0.982953 , 5.99373 ]), Female allosomes
960 , -96.3622 , array([ 0.761489 , -0.986512 , 8.28143 , 0.883072 , -0.997453 , 6.88913 , 0.146071 , -0.982953 , 5.99373 ]), Male allosomes
990 , -199.902 , array([ 0.761489 , -0.988651 , 8.28143 , 0.883072 , -0.997853 , 6.88913 , 0.146071 , -0.985639 , 5.99373 ]), Female allosomes
990 , -96.3588 , array([ 0.761489 , -0.988651 , 8.28143 , 0.883072 , -0.997853 , 6.88913 , 0.146071 , -0.985639 , 5.99373 ]), Male allosomes
1000 , -199.899 , array([ 0.761489 , -0.989281 , 8.28143 , 0.883072 , -0.997973 , 6.88913 , 0.146071 , -0.986434 , 5.99373 ]), Female allosomes
1000 , -96.3578 , array([ 0.761489 , -0.989281 , 8.28143 , 0.883072 , -0.997973 , 6.88913 , 0.146071 , -0.986434 , 5.99373 ]), Male allosomes
Optimization completed.
-----------------------
Optimization run #5
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -215.844 , array([ 0.761489 , -0.0384447 , 8.28143 , 0.883072 , 0.0652884 , 6.88913 , 0.146071 , 0.165736 , 5.99373 ]), Female allosomes
30 , -95.9915 , array([ 0.761489 , -0.0384447 , 8.28143 , 0.883072 , 0.0652884 , 6.88913 , 0.146071 , 0.165736 , 5.99373 ]), Male allosomes
60 , -213.363 , array([ 0.761489 , -0.0911971 , 8.28143 , 0.883072 , -0.07342 , 6.88913 , 0.146071 , 0.177397 , 5.99373 ]), Female allosomes
60 , -96.3466 , array([ 0.761489 , -0.0911971 , 8.28143 , 0.883072 , -0.07342 , 6.88913 , 0.146071 , 0.177397 , 5.99373 ]), Male allosomes
90 , -211.134 , array([ 0.761489 , -0.144072 , 8.28143 , 0.883072 , -0.20911 , 6.88913 , 0.146071 , 0.186119 , 5.99373 ]), Female allosomes
90 , -96.7195 , array([ 0.761489 , -0.144072 , 8.28143 , 0.883072 , -0.20911 , 6.88913 , 0.146071 , 0.186119 , 5.99373 ]), Male allosomes
120 , -209.208 , array([ 0.761489 , -0.192941 , 8.28143 , 0.883072 , -0.336671 , 6.88913 , 0.146071 , 0.189345 , 5.99373 ]), Female allosomes
120 , -97.0822 , array([ 0.761489 , -0.192941 , 8.28143 , 0.883072 , -0.336671 , 6.88913 , 0.146071 , 0.189345 , 5.99373 ]), Male allosomes
150 , -207.586 , array([ 0.761489 , -0.2365 , 8.28143 , 0.883072 , -0.453706 , 6.88913 , 0.146071 , 0.17879 , 5.99373 ]), Female allosomes
150 , -97.3964 , array([ 0.761489 , -0.2365 , 8.28143 , 0.883072 , -0.453706 , 6.88913 , 0.146071 , 0.17879 , 5.99373 ]), Male allosomes
180 , -206.228 , array([ 0.761489 , -0.286446 , 8.28143 , 0.883072 , -0.55499 , 6.88913 , 0.146071 , 0.159534 , 5.99373 ]), Female allosomes
180 , -97.6546 , array([ 0.761489 , -0.286446 , 8.28143 , 0.883072 , -0.55499 , 6.88913 , 0.146071 , 0.159534 , 5.99373 ]), Male allosomes
210 , -205.136 , array([ 0.761489 , -0.344374 , 8.28143 , 0.883072 , -0.637496 , 6.88913 , 0.146071 , 0.124609 , 5.99373 ]), Female allosomes
210 , -97.8215 , array([ 0.761489 , -0.344374 , 8.28143 , 0.883072 , -0.637496 , 6.88913 , 0.146071 , 0.124609 , 5.99373 ]), Male allosomes
240 , -204.278 , array([ 0.761489 , -0.392142 , 8.28143 , 0.883072 , -0.704259 , 6.88913 , 0.146071 , 0.0842843 , 5.99373 ]), Female allosomes
240 , -97.8943 , array([ 0.761489 , -0.392142 , 8.28143 , 0.883072 , -0.704259 , 6.88913 , 0.146071 , 0.0842843 , 5.99373 ]), Male allosomes
270 , -203.57 , array([ 0.761489 , -0.44391 , 8.28143 , 0.883072 , -0.759708 , 6.88913 , 0.146071 , 0.0259366 , 5.99373 ]), Female allosomes
270 , -97.914 , array([ 0.761489 , -0.44391 , 8.28143 , 0.883072 , -0.759708 , 6.88913 , 0.146071 , 0.0259366 , 5.99373 ]), Male allosomes
300 , -203 , array([ 0.761489 , -0.496911 , 8.28143 , 0.883072 , -0.803286 , 6.88913 , 0.146071 , -0.0434763 , 5.99373 ]), Female allosomes
300 , -97.8679 , array([ 0.761489 , -0.496911 , 8.28143 , 0.883072 , -0.803286 , 6.88913 , 0.146071 , -0.0434763 , 5.99373 ]), Male allosomes
330 , -202.542 , array([ 0.761489 , -0.543158 , 8.28143 , 0.883072 , -0.837425 , 6.88913 , 0.146071 , -0.122183 , 5.99373 ]), Female allosomes
330 , -97.7831 , array([ 0.761489 , -0.543158 , 8.28143 , 0.883072 , -0.837425 , 6.88913 , 0.146071 , -0.122183 , 5.99373 ]), Male allosomes
360 , -202.172 , array([ 0.761489 , -0.584853 , 8.28143 , 0.883072 , -0.865464 , 6.88913 , 0.146071 , -0.207128 , 5.99373 ]), Female allosomes
360 , -97.6504 , array([ 0.761489 , -0.584853 , 8.28143 , 0.883072 , -0.865464 , 6.88913 , 0.146071 , -0.207128 , 5.99373 ]), Male allosomes
390 , -201.858 , array([ 0.761489 , -0.630202 , 8.28143 , 0.883072 , -0.886982 , 6.88913 , 0.146071 , -0.292492 , 5.99373 ]), Female allosomes
390 , -97.5042 , array([ 0.761489 , -0.630202 , 8.28143 , 0.883072 , -0.886982 , 6.88913 , 0.146071 , -0.292492 , 5.99373 ]), Male allosomes
420 , -201.597 , array([ 0.761489 , -0.66937 , 8.28143 , 0.883072 , -0.904811 , 6.88913 , 0.146071 , -0.377741 , 5.99373 ]), Female allosomes
420 , -97.3509 , array([ 0.761489 , -0.66937 , 8.28143 , 0.883072 , -0.904811 , 6.88913 , 0.146071 , -0.377741 , 5.99373 ]), Male allosomes
450 , -201.37 , array([ 0.761489 , -0.704875 , 8.28143 , 0.883072 , -0.919812 , 6.88913 , 0.146071 , -0.457497 , 5.99373 ]), Female allosomes
450 , -97.2076 , array([ 0.761489 , -0.704875 , 8.28143 , 0.883072 , -0.919812 , 6.88913 , 0.146071 , -0.457497 , 5.99373 ]), Male allosomes
480 , -201.168 , array([ 0.761489 , -0.739744 , 8.28143 , 0.883072 , -0.932179 , 6.88913 , 0.146071 , -0.529707 , 5.99373 ]), Female allosomes
480 , -97.0801 , array([ 0.761489 , -0.739744 , 8.28143 , 0.883072 , -0.932179 , 6.88913 , 0.146071 , -0.529707 , 5.99373 ]), Male allosomes
510 , -200.992 , array([ 0.761489 , -0.768707 , 8.28143 , 0.883072 , -0.94306 , 6.88913 , 0.146071 , -0.595916 , 5.99373 ]), Female allosomes
510 , -96.9645 , array([ 0.761489 , -0.768707 , 8.28143 , 0.883072 , -0.94306 , 6.88913 , 0.146071 , -0.595916 , 5.99373 ]), Male allosomes
540 , -200.844 , array([ 0.761489 , -0.79586 , 8.28143 , 0.883072 , -0.951666 , 6.88913 , 0.146071 , -0.657017 , 5.99373 ]), Female allosomes
540 , -96.8578 , array([ 0.761489 , -0.79586 , 8.28143 , 0.883072 , -0.951666 , 6.88913 , 0.146071 , -0.657017 , 5.99373 ]), Male allosomes
570 , -200.717 , array([ 0.761489 , -0.820266 , 8.28143 , 0.883072 , -0.958775 , 6.88913 , 0.146071 , -0.711551 , 5.99373 ]), Female allosomes
570 , -96.764 , array([ 0.761489 , -0.820266 , 8.28143 , 0.883072 , -0.958775 , 6.88913 , 0.146071 , -0.711551 , 5.99373 ]), Male allosomes
600 , -200.6 , array([ 0.761489 , -0.844483 , 8.28143 , 0.883072 , -0.964681 , 6.88913 , 0.146071 , -0.755478 , 5.99373 ]), Female allosomes
600 , -96.6927 , array([ 0.761489 , -0.844483 , 8.28143 , 0.883072 , -0.964681 , 6.88913 , 0.146071 , -0.755478 , 5.99373 ]), Male allosomes
630 , -200.498 , array([ 0.761489 , -0.864994 , 8.28143 , 0.883072 , -0.969886 , 6.88913 , 0.146071 , -0.794182 , 5.99373 ]), Female allosomes
630 , -96.6309 , array([ 0.761489 , -0.864994 , 8.28143 , 0.883072 , -0.969886 , 6.88913 , 0.146071 , -0.794182 , 5.99373 ]), Male allosomes
660 , -200.406 , array([ 0.761489 , -0.882812 , 8.28143 , 0.883072 , -0.97457 , 6.88913 , 0.146071 , -0.82674 , 5.99373 ]), Female allosomes
660 , -96.5808 , array([ 0.761489 , -0.882812 , 8.28143 , 0.883072 , -0.97457 , 6.88913 , 0.146071 , -0.82674 , 5.99373 ]), Male allosomes
690 , -200.324 , array([ 0.761489 , -0.899136 , 8.28143 , 0.883072 , -0.978609 , 6.88913 , 0.146071 , -0.853215 , 5.99373 ]), Female allosomes
690 , -96.5421 , array([ 0.761489 , -0.899136 , 8.28143 , 0.883072 , -0.978609 , 6.88913 , 0.146071 , -0.853215 , 5.99373 ]), Male allosomes
720 , -200.254 , array([ 0.761489 , -0.913279 , 8.28143 , 0.883072 , -0.981976 , 6.88913 , 0.146071 , -0.876027 , 5.99373 ]), Female allosomes
720 , -96.5089 , array([ 0.761489 , -0.913279 , 8.28143 , 0.883072 , -0.981976 , 6.88913 , 0.146071 , -0.876027 , 5.99373 ]), Male allosomes
750 , -200.195 , array([ 0.761489 , -0.925418 , 8.28143 , 0.883072 , -0.984804 , 6.88913 , 0.146071 , -0.895621 , 5.99373 ]), Female allosomes
750 , -96.4807 , array([ 0.761489 , -0.925418 , 8.28143 , 0.883072 , -0.984804 , 6.88913 , 0.146071 , -0.895621 , 5.99373 ]), Male allosomes
780 , -200.144 , array([ 0.761489 , -0.936159 , 8.28143 , 0.883072 , -0.987127 , 6.88913 , 0.146071 , -0.91228 , 5.99373 ]), Female allosomes
780 , -96.457 , array([ 0.761489 , -0.936159 , 8.28143 , 0.883072 , -0.987127 , 6.88913 , 0.146071 , -0.91228 , 5.99373 ]), Male allosomes
810 , -200.099 , array([ 0.761489 , -0.945562 , 8.28143 , 0.883072 , -0.989194 , 6.88913 , 0.146071 , -0.925661 , 5.99373 ]), Female allosomes
810 , -96.4389 , array([ 0.761489 , -0.945562 , 8.28143 , 0.883072 , -0.989194 , 6.88913 , 0.146071 , -0.925661 , 5.99373 ]), Male allosomes
840 , -200.064 , array([ 0.761489 , -0.953167 , 8.28143 , 0.883072 , -0.990799 , 6.88913 , 0.146071 , -0.938275 , 5.99373 ]), Female allosomes
840 , -96.4207 , array([ 0.761489 , -0.953167 , 8.28143 , 0.883072 , -0.990799 , 6.88913 , 0.146071 , -0.938275 , 5.99373 ]), Male allosomes
870 , -200.031 , array([ 0.761489 , -0.960368 , 8.28143 , 0.883072 , -0.992216 , 6.88913 , 0.146071 , -0.948 , 5.99373 ]), Female allosomes
870 , -96.4077 , array([ 0.761489 , -0.960368 , 8.28143 , 0.883072 , -0.992216 , 6.88913 , 0.146071 , -0.948 , 5.99373 ]), Male allosomes
900 , -200.003 , array([ 0.761489 , -0.966521 , 8.28143 , 0.883072 , -0.99342 , 6.88913 , 0.146071 , -0.955889 , 5.99373 ]), Female allosomes
900 , -96.3975 , array([ 0.761489 , -0.966521 , 8.28143 , 0.883072 , -0.99342 , 6.88913 , 0.146071 , -0.955889 , 5.99373 ]), Male allosomes
930 , -199.98 , array([ 0.761489 , -0.97178 , 8.28143 , 0.883072 , -0.994442 , 6.88913 , 0.146071 , -0.962679 , 5.99373 ]), Female allosomes
930 , -96.3887 , array([ 0.761489 , -0.97178 , 8.28143 , 0.883072 , -0.994442 , 6.88913 , 0.146071 , -0.962679 , 5.99373 ]), Male allosomes
960 , -199.96 , array([ 0.761489 , -0.976123 , 8.28143 , 0.883072 , -0.995308 , 6.88913 , 0.146071 , -0.968534 , 5.99373 ]), Female allosomes
960 , -96.3811 , array([ 0.761489 , -0.976123 , 8.28143 , 0.883072 , -0.995308 , 6.88913 , 0.146071 , -0.968534 , 5.99373 ]), Male allosomes
990 , -199.942 , array([ 0.761489 , -0.979765 , 8.28143 , 0.883072 , -0.996073 , 6.88913 , 0.146071 , -0.973306 , 5.99373 ]), Female allosomes
990 , -96.375 , array([ 0.761489 , -0.979765 , 8.28143 , 0.883072 , -0.996073 , 6.88913 , 0.146071 , -0.973306 , 5.99373 ]), Male allosomes
1000 , -199.937 , array([ 0.761489 , -0.980863 , 8.28143 , 0.883072 , -0.99629 , 6.88913 , 0.146071 , -0.974805 , 5.99373 ]), Female allosomes
1000 , -96.373 , array([ 0.761489 , -0.980863 , 8.28143 , 0.883072 , -0.99629 , 6.88913 , 0.146071 , -0.974805 , 5.99373 ]), Male allosomes
Optimization completed.
-----------------------
In Step 2: Results from multiple optimization runs with different starting parameters:
-------------------------------------------------------------------------------------------------------------------------------------------------------------
Run | LogLik | REUR | REUR_sex_bias | t1 | RAFR | RAFR_sex_bias | t2 | REUR2 | REUR2_sex_bias | t3
-------------------------------------------------------------------------------------------------------------------------------------------------------------
1 | -296.315 | 0.7615 | -0.9803 | 8.281 | 0.8831 | -0.9961 | 6.889 | 0.1461 | -0.974 | 5.994
2 | -296.308 | 0.7615 | -0.9809 | 8.281 | 0.8831 | -0.9964 | 6.889 | 0.1461 | -0.9753 | 5.994
3 | -296.283 | 0.7615 | -0.9852 | 8.281 | 0.8831 | -0.9972 | 6.889 | 0.1461 | -0.9802 | 5.994
4 | -296.257 | 0.7615 | -0.9893 | 8.281 | 0.8831 | -0.998 | 6.889 | 0.1461 | -0.9864 | 5.994
5 | -296.31 | 0.7615 | -0.9809 | 8.281 | 0.8831 | -0.9963 | 6.889 | 0.1461 | -0.9748 | 5.994
-------------------------------------------------------------------------------------------------------------------------------------------------------------
Selecting best parameters from step 2.
Step 2 used allosomal data only. Final likelihood is evaluated on autosomal + allosomal data at the selected optimal parameters.
--------------------------------------------------------------------------------------------------
Launching re-optimization until convergence is achieved or 5 re-optimizations have been performed.
--------------------------------------------------------------------------------------------------
Re-optimization 1/5: re-optimizing starting from the current optimal parameters (likelihood = -965.051689).
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -668.343 , array([ 0.76093 , -0.989281 , 8.25567 , 0.883314 , -0.997973 , 6.97332 , 0.146155 , -0.986434 , 5.98967 ]), Autosomes
43 , -666.826 , array([ 0.760915 , -0.989281 , 8.25578 , 0.883298 , -0.997973 , 6.97519 , 0.146159 , -0.986434 , 5.99338 ]), Autosomes
Optimization completed.
-----------------------
----------------------------------------------------------------------------------------------------------------
Step 2 : Optimizing allosomal likelihood over parameters : ['REUR_sex_bias', 'RAFR_sex_bias', 'REUR2_sex_bias'].
----------------------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -198.769 , array([ 0.760915 , -0.990654 , 8.25578 , 0.883298 , -0.998233 , 6.97519 , 0.146159 , -0.988495 , 5.99338 ]), Female allosomes
30 , -96.6165 , array([ 0.760915 , -0.990654 , 8.25578 , 0.883298 , -0.998233 , 6.97519 , 0.146159 , -0.988495 , 5.99338 ]), Male allosomes
60 , -198.762 , array([ 0.760915 , -0.991923 , 8.25578 , 0.883298 , -0.998494 , 6.97519 , 0.146159 , -0.990584 , 5.99338 ]), Female allosomes
60 , -96.6136 , array([ 0.760915 , -0.991923 , 8.25578 , 0.883298 , -0.998494 , 6.97519 , 0.146159 , -0.990584 , 5.99338 ]), Male allosomes
90 , -198.756 , array([ 0.760915 , -0.993156 , 8.25578 , 0.883298 , -0.998714 , 6.97519 , 0.146159 , -0.992179 , 5.99338 ]), Female allosomes
90 , -96.6115 , array([ 0.760915 , -0.993156 , 8.25578 , 0.883298 , -0.998714 , 6.97519 , 0.146159 , -0.992179 , 5.99338 ]), Male allosomes
120 , -198.752 , array([ 0.760915 , -0.993878 , 8.25578 , 0.883298 , -0.998856 , 6.97519 , 0.146159 , -0.993078 , 5.99338 ]), Female allosomes
120 , -96.6104 , array([ 0.760915 , -0.993878 , 8.25578 , 0.883298 , -0.998856 , 6.97519 , 0.146159 , -0.993078 , 5.99338 ]), Male allosomes
144 , -198.752 , array([ 0.760915 , -0.993906 , 8.25578 , 0.883298 , -0.998858 , 6.97519 , 0.146159 , -0.993112 , 5.99338 ]), Female allosomes
144 , -96.6103 , array([ 0.760915 , -0.993906 , 8.25578 , 0.883298 , -0.998858 , 6.97519 , 0.146159 , -0.993112 , 5.99338 ]), Male allosomes
Optimization completed.
-----------------------
Change in likelihood from -965.051689 to -962.190994 after re-optimizing.
Re-optimization 2/5: re-optimizing starting from the current optimal parameters (likelihood = -962.190994).
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -666.684 , array([ 0.76063 , -0.993906 , 8.2447 , 0.883298 , -0.998858 , 6.97209 , 0.146162 , -0.993112 , 5.99393 ]), Autosomes
36 , -666.609 , array([ 0.760631 , -0.993906 , 8.24446 , 0.883292 , -0.998858 , 6.97302 , 0.146159 , -0.993112 , 5.99353 ]), Autosomes
Optimization completed.
-----------------------
----------------------------------------------------------------------------------------------------------------
Step 2 : Optimizing allosomal likelihood over parameters : ['REUR_sex_bias', 'RAFR_sex_bias', 'REUR2_sex_bias'].
----------------------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -198.834 , array([ 0.760631 , -0.990647 , 8.24446 , 0.883292 , -0.998236 , 6.97302 , 0.146159 , -0.988491 , 5.99353 ]), Female allosomes
30 , -96.5967 , array([ 0.760631 , -0.990647 , 8.24446 , 0.883292 , -0.998236 , 6.97302 , 0.146159 , -0.988491 , 5.99353 ]), Male allosomes
60 , -198.827 , array([ 0.760631 , -0.992034 , 8.24446 , 0.883292 , -0.998477 , 6.97302 , 0.146159 , -0.990569 , 5.99353 ]), Female allosomes
60 , -96.5938 , array([ 0.760631 , -0.992034 , 8.24446 , 0.883292 , -0.998477 , 6.97302 , 0.146159 , -0.990569 , 5.99353 ]), Male allosomes
90 , -198.82 , array([ 0.760631 , -0.993196 , 8.24446 , 0.883292 , -0.99871 , 6.97302 , 0.146159 , -0.992182 , 5.99353 ]), Female allosomes
90 , -96.5916 , array([ 0.760631 , -0.993196 , 8.24446 , 0.883292 , -0.99871 , 6.97302 , 0.146159 , -0.992182 , 5.99353 ]), Male allosomes
120 , -198.817 , array([ 0.760631 , -0.993935 , 8.24446 , 0.883292 , -0.998857 , 6.97302 , 0.146159 , -0.993084 , 5.99353 ]), Female allosomes
120 , -96.5905 , array([ 0.760631 , -0.993935 , 8.24446 , 0.883292 , -0.998857 , 6.97302 , 0.146159 , -0.993084 , 5.99353 ]), Male allosomes
139 , -198.816 , array([ 0.760631 , -0.993954 , 8.24446 , 0.883292 , -0.998857 , 6.97302 , 0.146159 , -0.993113 , 5.99353 ]), Female allosomes
139 , -96.5904 , array([ 0.760631 , -0.993954 , 8.24446 , 0.883292 , -0.998857 , 6.97302 , 0.146159 , -0.993113 , 5.99353 ]), Male allosomes
Optimization completed.
-----------------------
Change in likelihood from -962.190994 to -962.015782 after re-optimizing.
Re-optimization 3/5: re-optimizing starting from the current optimal parameters (likelihood = -962.015782).
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -666.367 , array([ 0.760365 , -0.993954 , 8.23446 , 0.883292 , -0.998857 , 6.97097 , 0.146158 , -0.993113 , 5.99335 ]), Autosomes
34 , -666.33 , array([ 0.760347 , -0.993954 , 8.23446 , 0.883291 , -0.998857 , 6.97097 , 0.146152 , -0.993113 , 5.99323 ]), Autosomes
Optimization completed.
-----------------------
----------------------------------------------------------------------------------------------------------------
Step 2 : Optimizing allosomal likelihood over parameters : ['REUR_sex_bias', 'RAFR_sex_bias', 'REUR2_sex_bias'].
----------------------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -198.896 , array([ 0.760347 , -0.990757 , 8.23446 , 0.883291 , -0.998225 , 6.97097 , 0.146152 , -0.988419 , 5.99323 ]), Female allosomes
30 , -96.5792 , array([ 0.760347 , -0.990757 , 8.23446 , 0.883291 , -0.998225 , 6.97097 , 0.146152 , -0.988419 , 5.99323 ]), Male allosomes
60 , -198.889 , array([ 0.760347 , -0.9921 , 8.23446 , 0.883291 , -0.99847 , 6.97097 , 0.146152 , -0.990524 , 5.99323 ]), Female allosomes
60 , -96.5762 , array([ 0.760347 , -0.9921 , 8.23446 , 0.883291 , -0.99847 , 6.97097 , 0.146152 , -0.990524 , 5.99323 ]), Male allosomes
90 , -198.883 , array([ 0.760347 , -0.993298 , 8.23446 , 0.883291 , -0.998697 , 6.97097 , 0.146152 , -0.992133 , 5.99323 ]), Female allosomes
90 , -96.574 , array([ 0.760347 , -0.993298 , 8.23446 , 0.883291 , -0.998697 , 6.97097 , 0.146152 , -0.992133 , 5.99323 ]), Male allosomes
120 , -198.947 , array([ 0.760347 , -0.99407 , 8.23446 , 0.883291 , -0.998832 , 6.97097 , 0.146152 , -0.993114 , 5.99323 ]), Female allosomes
120 , -96.5443 , array([ 0.760347 , -0.99407 , 8.23446 , 0.883291 , -0.998832 , 6.97097 , 0.146152 , -0.993114 , 5.99323 ]), Male allosomes
145 , -198.879 , array([ 0.760347 , -0.994072 , 8.23446 , 0.883291 , -0.998832 , 6.97097 , 0.146152 , -0.993111 , 5.99323 ]), Female allosomes
145 , -96.5727 , array([ 0.760347 , -0.994072 , 8.23446 , 0.883291 , -0.998832 , 6.97097 , 0.146152 , -0.993111 , 5.99323 ]), Male allosomes
Optimization completed.
-----------------------
Change in likelihood from -962.015782 to -961.781942 after re-optimizing.
Re-optimization 4/5: re-optimizing starting from the current optimal parameters (likelihood = -961.781942).
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -666.06 , array([ 0.760063 , -0.994072 , 8.22209 , 0.883283 , -0.998832 , 6.96792 , 0.146106 , -0.993111 , 5.99341 ]), Autosomes
38 , -666.06 , array([ 0.760049 , -0.994072 , 8.222 , 0.88328 , -0.998832 , 6.96826 , 0.146115 , -0.993111 , 5.99341 ]), Autosomes
Optimization completed.
-----------------------
----------------------------------------------------------------------------------------------------------------
Step 2 : Optimizing allosomal likelihood over parameters : ['REUR_sex_bias', 'RAFR_sex_bias', 'REUR2_sex_bias'].
----------------------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -198.974 , array([ 0.760049 , -0.990685 , 8.222 , 0.88328 , -0.99823 , 6.96826 , 0.146115 , -0.988504 , 5.99341 ]), Female allosomes
30 , -96.5555 , array([ 0.760049 , -0.990685 , 8.222 , 0.88328 , -0.99823 , 6.96826 , 0.146115 , -0.988504 , 5.99341 ]), Male allosomes
60 , -198.967 , array([ 0.760049 , -0.992042 , 8.222 , 0.88328 , -0.998484 , 6.96826 , 0.146115 , -0.990538 , 5.99341 ]), Female allosomes
60 , -96.5526 , array([ 0.760049 , -0.992042 , 8.222 , 0.88328 , -0.998484 , 6.96826 , 0.146115 , -0.990538 , 5.99341 ]), Male allosomes
90 , -198.961 , array([ 0.760049 , -0.993261 , 8.222 , 0.88328 , -0.998706 , 6.96826 , 0.146115 , -0.992155 , 5.99341 ]), Female allosomes
90 , -96.5504 , array([ 0.760049 , -0.993261 , 8.222 , 0.88328 , -0.998706 , 6.96826 , 0.146115 , -0.992155 , 5.99341 ]), Male allosomes
120 , -198.957 , array([ 0.760049 , -0.99403 , 8.222 , 0.88328 , -0.998845 , 6.96826 , 0.146115 , -0.993108 , 5.99341 ]), Female allosomes
120 , -96.5491 , array([ 0.760049 , -0.99403 , 8.222 , 0.88328 , -0.998845 , 6.96826 , 0.146115 , -0.993108 , 5.99341 ]), Male allosomes
141 , -198.957 , array([ 0.760049 , -0.994037 , 8.222 , 0.88328 , -0.998847 , 6.96826 , 0.146115 , -0.993111 , 5.99341 ]), Female allosomes
141 , -96.5491 , array([ 0.760049 , -0.994037 , 8.222 , 0.88328 , -0.998847 , 6.96826 , 0.146115 , -0.993111 , 5.99341 ]), Male allosomes
Optimization completed.
-----------------------
Change in likelihood from -961.781942 to -961.565994 after re-optimizing.
Re-optimization 5/5: re-optimizing starting from the current optimal parameters (likelihood = -961.565994).
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -665.801 , array([ 0.759764 , -0.994037 , 8.2121 , 0.883302 , -0.998847 , 6.96586 , 0.14606 , -0.993111 , 5.99355 ]), Autosomes
38 , -665.801 , array([ 0.75975 , -0.994037 , 8.2127 , 0.883302 , -0.998847 , 6.96608 , 0.146059 , -0.993111 , 5.99359 ]), Autosomes
Optimization completed.
-----------------------
----------------------------------------------------------------------------------------------------------------
Step 2 : Optimizing allosomal likelihood over parameters : ['REUR_sex_bias', 'RAFR_sex_bias', 'REUR2_sex_bias'].
----------------------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -199.049 , array([ 0.75975 , -0.990678 , 8.2127 , 0.883302 , -0.99823 , 6.96608 , 0.146059 , -0.988513 , 5.99359 ]), Female allosomes
30 , -96.5344 , array([ 0.75975 , -0.990678 , 8.2127 , 0.883302 , -0.99823 , 6.96608 , 0.146059 , -0.988513 , 5.99359 ]), Male allosomes
60 , -199.042 , array([ 0.75975 , -0.992077 , 8.2127 , 0.883302 , -0.998487 , 6.96608 , 0.146059 , -0.990489 , 5.99359 ]), Female allosomes
60 , -96.5316 , array([ 0.75975 , -0.992077 , 8.2127 , 0.883302 , -0.998487 , 6.96608 , 0.146059 , -0.990489 , 5.99359 ]), Male allosomes
90 , -199.036 , array([ 0.75975 , -0.99329 , 8.2127 , 0.883302 , -0.998701 , 6.96608 , 0.146059 , -0.992131 , 5.99359 ]), Female allosomes
90 , -96.5293 , array([ 0.75975 , -0.99329 , 8.2127 , 0.883302 , -0.998701 , 6.96608 , 0.146059 , -0.992131 , 5.99359 ]), Male allosomes
120 , -199.1 , array([ 0.75975 , -0.994032 , 8.2127 , 0.883302 , -0.998837 , 6.96608 , 0.146059 , -0.993111 , 5.99359 ]), Female allosomes
120 , -96.4995 , array([ 0.75975 , -0.994032 , 8.2127 , 0.883302 , -0.998837 , 6.96608 , 0.146059 , -0.993111 , 5.99359 ]), Male allosomes
137 , -199.032 , array([ 0.75975 , -0.994032 , 8.2127 , 0.883302 , -0.998837 , 6.96608 , 0.146059 , -0.993109 , 5.99359 ]), Female allosomes
137 , -96.5279 , array([ 0.75975 , -0.994032 , 8.2127 , 0.883302 , -0.998837 , 6.96608 , 0.146059 , -0.993109 , 5.99359 ]), Male allosomes
Optimization completed.
-----------------------
Change in likelihood from -961.565994 to -961.360865 after re-optimizing.
Convergence not achieved after 5 repetitions. Stopping re-optimization.
Final parameters and corresponding likelihood computed on autosomal + allosomal data:
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
LogLik | REUR | REUR_sex_bias | t1 | RAFR | RAFR_sex_bias | t2 | REUR2 | REUR2_sex_bias | t3 | targetpop_NAT_rate | targetpop_NAT_sex_bias
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
-961.361 | 0.7598 | -0.994 | 8.213 | 0.8833 | -0.9988 | 6.966 | 0.1461 | -0.9931 | 5.994 | 0.2402 | 0.994
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Parameters targetpop_NAT_rate, targetpop_NAT_sex_bias correspond to the dependent ancestry and were not free in the optimization.
The optimal solution has sex-bias parameter(s) REUR_sex_bias, RAFR_sex_bias, REUR2_sex_bias, targetpop_NAT_sex_bias near their ±1 boundary.
The implicit population's sex-bias parameter (targetpop_NAT_sex_bias) is near the +-1 boundary,
but every other source population in the same founder event either also has its sex-bias parameter
at a boundary or is already fixed by value: no alternate implicit population can be chosen.
Re-optimizing with sex-bias parameter(s) REUR_sex_bias, RAFR_sex_bias, REUR2_sex_bias fixed at their boundary value.
-----------------------------------------------------------------------------------------------------
Step 1 : Optimizing autosomal likelihood over parameters ['REUR', 't1', 'RAFR', 't2', 'REUR2', 't3'].
-----------------------------------------------------------------------------------------------------
Iter. Log-likelihood Model parameters Transmission
-------------------------------------------------------------
30 , -665.547 , array([ 0.759465 , -0.999 , 8.20032 , 0.883293 , -0.999 , 6.96313 , 0.146014 , -0.999 , 5.99378 ]), Autosomes
41 , -665.46 , array([ 0.759449 , -0.999 , 8.20018 , 0.883289 , -0.999 , 6.96354 , 0.14602 , -0.999 , 5.99321 ]), Autosomes
Optimization completed.
-----------------------
All sex-bias parameters are fixed (REUR_sex_bias, RAFR_sex_bias, REUR2_sex_bias by user-provided values). Step 2 has no free parameters to optimize and will be skipped.
Final parameters and corresponding likelihood computed on autosomal + allosomal data:
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
LogLik | REUR | REUR_sex_bias | t1 | RAFR | RAFR_sex_bias | t2 | REUR2 | REUR2_sex_bias | t3 | targetpop_NAT_rate | targetpop_NAT_sex_bias
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
-961.058 | 0.7594 | -0.999 | 8.2 | 0.8833 | -0.999 | 6.964 | 0.146 | -0.999 | 5.993 | 0.2406 | 0.999
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Parameters targetpop_NAT_rate, targetpop_NAT_sex_bias correspond to the dependent ancestry and were not free in the optimization.
All sex-bias parameters have been fixed by value. Boundary re-optimization completed.
Predicted autosome proportions: [0.22108777 0.02373156 0.75518067]
Predicted allosome proportions: [0.19749592 0.0316578 0.77084628]
Results saved to : output_three_pulses
{'destination_dir': PosixPath('/home/runner/work/tracts/tracts/docs/source/auto_examples/ASW/output_three_pulses'), 'table_file': PosixPath('/home/runner/work/tracts/tracts/docs/source/auto_examples/ASW/output_three_pulses/ASW_test_output_optimal_parameters.txt')}
import sys
from pathlib import Path
from tracts.driver import run_tracts
# Read files automatically for online documentation
sys.path.append('.')
script_dir = Path.cwd()
driver_filename = script_dir / "ASW_three_pulses.yaml"
run_tracts(
driver_filename=str(driver_filename),
script_dir=str(script_dir),
)
Total running time of the script: (72 minutes 52.534 seconds)